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ABSTRACT 


The  Naval  Postgraduate  School  possesses  software  designed  to  perform  thermal 
analysis  of  electronic  components.  At  the  core  of  this  package  is  a  model  builder  whose 
purpose  is  to  generate  a  thermal  model  for  use  in  steady  state  and  transient  thermal 
analyzers.  The  current  version  of  the  model  builder  requires  excessive  amounts  of  time 
for  data  input  and  model  verification.  This  thesis  describes  the  development  of  a  model 
builder  specifically  designed  to  reduce  the  time  required  to  model  a  printed  circuit  board 
containing  up  to  four  copper  layers. 
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I.  INTRODUCTION 


The  Naval  Postgraduate  School  possesses  software  designed  to  perform  thermal 
analysis  of  electronic  components.  The  software  package  contains  two  routines  de¬ 
signed  to  generate  a  thermal  model  (in  the  form  of  an  ASCII  input  data  file)  to  be  read 
by  the  thermal  analyzer  program.  The  first  routine  is  a  generalized  model  builder  used 
in  developmental  stages,  as  well  as  an  editor  used  to  modify  existing  models.  The  second 
routine  was  designed  to  work  with  models  that  have  a  specific  geometric  configuration. 

In  order  to  accurately  model  electrical  components,  the  structure  must  be  subdivided 
into  small,  equal  sized  subvolumes.  The  centroid  of  each  subvolume  is  referred  to  as  a 
node,  and  due  to  the  assumed  isothermal  nature  of  each  subvolume,  a  node  is  considered 
as  representative  of  the  total  subvolume. 

Producing  an  accurate  thermal  model,  requires  the  design  engineer  to  deal  with 
enormous  amounts  of  equations  and  temperatures,  each  describing  an  indiviual  node. 
Additionally,  each  node  is  connected  to  adjacent  nodes  by  thermal  conductances.  As  the 
desired  accuracy  increases,  the  number  of  nodes  and  equations  also  increases.  Modeling 
the  electronic  component  without  the  aid  of  a  computerized  model  builder  is  a  task  re¬ 
quiring  inordinate  amounts  of  time. 

Current  versions  of  the  thermal  analyzer  contain  provisions  for  the  generation  of  the 
node  equations  and  node  interrelationships;  however,  data  input  is  still  a  manual,  time- 
consuming  process.  At  present  there  exists  a  need  for  a  program  that  will  generate  a 
data-input  file  for  the  thermal  analyzer  that  is  both  generated  through  a  menu-driven 
routine,  and  allows  the  design  engineer  enough  flexibility  to  model  the  electronic  com¬ 
ponent  to  suit  his  or  her  needs. 

This  thesis  describes  the  development  of  a  model  builder  designed  specifically  to  re¬ 
duce  the  time  required  to  model  the  copper  and  epoxy  layers  of  a  printed  circuit  board. 
A  typical  printed  circuit  board  configuration  is  shown  in  Figure  1.  The  printed  circuit 
board  may  contain  up  to  four  copper  layers.  Additional  features  of  the  model  builder 
covered  are: 

1.  The  capability  of  working  in  SI  or  English  units. 

2.  The  choice  of  a  total  of  sixteen  aspect  ratios. 

3.  The  provision  for  up  to  740  nodes. 

4.  The  ability  to  specify  the  percent  of  copper  coverage  for  each  layer. 


1 


5.  The  ability  to  specify  dimensional  parameters  for  each  individual  copper  layer. 

6.  The  ability  to  input  heat  dissipation  using  several  methods. 

7.  The  provision  for  six  ambient  temperatures. 

8.  The  automatic  calculation  of  conductance  values  based  on  user  input. 
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II.  A  JUSTIFICATION  FOR  THERMAL  ANALYSIS 


A  printed  circuit  board  is  a  conglomeration  of  organic  and  inorganic  materials  with 
external  and  internal  wiring,  which  allows  electronic  components  to  be  physically  sup¬ 
ported,  and  electrically  connected.  [Ref.  1] 

Over  the  past  several  decades,  printed  circuit  board  technology  has  developed  sub¬ 
stantially.  Early  printed  circuits  were  fabricated  by  printing  a  pattern  of  polymer  resist 
on  a  copper  plane  and  then  chemically  etching.  Holes  drilled  in  the  laminate  held  the 
component  leads  that  were  soldered  to  the  copper-printed  patterns.  The  technology  has 
progressed  in  developing  the  sophistication  and  uses  of  the  printed  circuit  board  inter¬ 
connections. 

Today,  the  basic  functions  of  the  printed  circuit  board  are  the  same;  the  intercon¬ 
necting  copper  signal  lines  join  two  I/O  leads  from  two  different  components.  The 
components  may  be  resistors,  inductors,  capacitors,  or  semiconductor  chips.  When 
applying  multichip  technologies,  there  may  be  hundreds  of  components  attached  to  the 
printed  circuit  board.  The  ever-increasing  level  of  complexity  of  printed  circuit  boards 
has  forced  primitive  boards  using  a  few  yards  of  printed  wiring  in  me  1960s,  to  evolve 
into  sophisticated,  multi-layered  structures  requiring  kilometers  of  printed  wiring.  This 
increase  in  the  circuit  board  level  of  sophistication  can  be  attibuted  to  the  integration 
of  semiconductors  and  an  increased  need  for  I/O  capabilities.  [Ref.  2] 

The  most  complex  printed  circuit  boards  contain  kilometers  of  copper  intercon¬ 
nection,  roughly  50  to  100  microns  wide,  and  half  as  thick.  These  boards  distribute 
KlVjm^  of  power  internally  in  very  densely  packed  layers  of  copper.  The  drive  for  higher 
performance  means  that  there  is  a  greater  requirement  for  power-handling  and  cooling 
capabilities.  [Ref  2] 

It  is  the  responsibility  of  the  designers  to  ensure  that  cooling  on  the  printed  circuit 
board  is  adequate  under  all  possible  load  conditions  in  order  to  allow-  proper  perform¬ 
ance  of  the  individual  components  and  the  board  as  a  whole.  Therefore,  it  is  imperative 
that  the  designers  understand  and  are  able  to  predict  the  temperature  distribution  on 
multilayered  structures  prior  to  prototype  production.  The  overriding  reasons  for  per¬ 
forming  a  precise  thermal  analysis  are  to  increase  component  reliability,  ensure  proper 
material  selection,  reduce  the  possibility  of  catastrophic  thermal  failure,  and  guarantee 
electrical  performance. 
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A.  RELIABILITY 

Reliability  is  defined  as  the  probability  that  a  component  is  functioning  as  designed, 
while  failure  is  defined  as  the  probability  that  a  component  is  not  functioning  as  de¬ 
signed.  There  is  a  predictable  relationship  between  the  operating  temperature  of  elec¬ 
tronic  components  and  reliability.  The  materials  used  in  the  fabrication  of  these 
components  have  thermal  limitations,  and  should  these  thermal  limitations  be  exceeded, 
the  physical  and  chemical  properties  of  the  material  are  affected,  and  the  device  fails. 
[Ref  1] 

For  a  large  number  of  components,  a  typical  plot  of  failure  rate  as  a  function  of  time 
is  shown  in  Figure  2.  Failures  at  short  times  are  called  early  fails,  or  infant  mortality, 
while  failures  at  long  times  are  called  wearout  fails  because  they  result  from  usage.  At 
all  times,  failures  can  occur  from  intrinsic  mechanisms,  or  from  random  overstress.  [Ref 

3] 

Provided  the  device  has  been  adequately  designed,  early  failures  can  arise  as  a  result 
of  manufacturing  defects.  Defects  that  occur  early  on,  or  the  "burn-in"  period,  are 
considered  to  be  the  result  of  poor  or  inadequate  quality  control  mechanisms  in  the 
manufacturing  process.  [Ref  2] 

Of  greatest  concern  are  the  failures  that  occur  during  the  useful  life  of  the  device 
because  the  probability  of  failure  during  this  period  should  be  nearly  zero.  Should  a 
device  fail  during  its  useful  life,  the  probable  cause  of  failure  would  be  due  to  a  variety 
of  external  factors,  and  are  unpredictable.  [Ref  2] 

As  time  and  usage  progress,  the  terminal  period  of  wearout  is  encountered,  usually 
well  past  the  the  end  of  system  life  (EOL),  in  which  the  failure  rate  increases.  Con¬ 
ductor  electromigration  is  a  typical  example  of  a  wearout  mechanism,  in  which  the 
electron  flow  itself  causes  irreversible  mass  flow,  which  causes  the  formation  of  voids  and 
consequent  conductor  failure.  Printed  circuit  boards  seldom  have  significant  failure  rates 
when  they  are  produced  with  sufficient  characterization  and  control.  When  failures  do 
occur,  they  can  normally  be  attributed  to  either  manufacturing  defects  or  lack  of  integ¬ 
rity  to  the  design  specifications.  The  most  notable  concerns  are  the  resin  and  laminate 
effects  that  can  give  rise  to  insulation  integrity.  As  printed  circuit  boards  become  more 
densely  packaged  and  continue  to  spread  into  more  unconventional  environments  such 
as  homes,  automobiles  or  marine  engines,  there  is  a  need  to  establish  improved  materials 
and  process  controls  for  even  better  reliability.  Thermal  analysis  of  the  printed  circuit 
board  becomes  one  of  the  fundamental  facets  of  the  design  process  in  order  to  better 
characterize  the  printed  circuit  board.  [Ref  2] 
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The  infant  mortality, • or  early  fail 
portion  of  the  curve  is  the  result 
of  failures  during  the  "burn-in" 
period.  The  high  failure  rate  shown 
on  the  right  side  of  the  curve  is  due 
to  usage,  or  "wear-out." 


Figure  2.  Failure  rate  as  a  function  of  time  for  typical  component. 
Source:  Reference  3. 
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B.  MATERIAL 


The  fabrication  of  printed  circuit  boards  results  in  the  joining  of  different  materials. 
The  materials  selected  can  have  a  significant  impact  in  thermal  properties  of  the  printed 
circuit  board.  In  the  largest  circuit  boards  hundreds  of  amps  may  be  switched  at  once. 
As  packaging  densities  increase  thermal,  mechanical,  electrical,  and  chemical  coupling 
becomes  very  strong.  In  view  of  the  complexity  of  electrical  structures  today,  computer 
modeling  of  total  thermal  response  of  a  printed  circuit  board  is  a  requirement  conducive 
to  understanding  one  of  many  interactions.  [Ref  2] 

C.  THERMAL  FATIGUE  AND  CATASTROPHIC  THERMAL  FAILURE 

Printed  circuit  boards  are  comprised  of  dissimilar  materials  that  expand  at  different 
rates  of  heating.  Table  1  shows  the  thermal  coefficients  of  expansion  of  materials  com¬ 
monly  used  in  the  fabrication  of  printed  circuit  boards.  The  differential  expansion  of 
mismatch  must  be  accommodated  by  the  various  elements  on  the  board.  The  increasing 
levels  of  packaging  densities  and  board  complexity  dictate  the  need  for  designing  a 
thermal  environment  that  can  accommodate  the  diverse  components  that  are  in  close 
proximity  to  each  other.  [Ref  4] 


Table  1.  MATERIALS  USED  IN  PCB  PRODUCTION 


Chemical  Name 

Coefficient  of  Thermal  Expansion  lO-’/'C 

Polyethylene 

650 

Bismaleimide 

500 

Polyamide  imide 

360 

Polyaryiate 

625 

Peek 

400 

Polyether  amide 

500 

Polyimide 

500 

Polytetrafluorethylene 

700-1200 

Epoxy  glass  cloth 

170  (x,y);  600  (z) 

PBZT 

90  (x,y);  200  (z) 

Source:  Reference  2. 
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As  previously  menlioned  it  is  very  important  that  the  board  designers  have  an 
understanding  of  the  operating  environment  in  which  the  board  will  be  operating  in  or¬ 
der  to  incorporate  into  the  design  the  tolerances  that  will  allow  the  product  to  operate 
reliably. 

Catastrophic  thermal  failure  is  defined  as  an  immediate,  thermally  induced,  total  loss 
of  electronic  function  in  a  specified  component.  This  type  of  failure  comes  as  a  result 
of  excessive  temperature,  or  a  thermal  fracture.  Catastrophic  failure  comes  about  as  a 
result  of  many  factors  including  the  operating  environment,  equipment  histoiy,  me¬ 
chanical  loading,  and  operational  modes  of  the  component.  Although  it  is  difficult  to 
predict  the  temperature  at  which  thermal  failure  may  occur,  it  is  possible  to  establish 
with  the  aid  of  thermal  analysis  the  boundaries  at  which  the  board  can  be  expected  to 
operate  reliably  and  within  its  useful  operating  life  cycle.  [Ref.  1,2] 
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III.  HEAT  TRANSFER 


The  degradation  of  the  heat  flow  capabilities  of  a  printed  circuit  board  can  lead  to 
reliability  problems  due  to  excessive  operating  temperatures.  It  is  imperative  that  de¬ 
signers  incorporate  into  the  board  the  capability  to  maintain  temperatures  within  upper 
operational  limits  while  operating  in  all  possible  environments  in  which  the  board  will 
be  exposed.  [Ref  2] 

Heat  transfer  is  defmed  as  all  energy  flows  that  arise  as  a  result  of  temperature  dif¬ 
ferences.  Because  the  components  mounted  on  printed  circuit  boards  and,  indeed,  the 
printed  circuit  boards  themselves  are  not  one  hundred  percent  efficient,  heat  is  gener¬ 
ated.  The  primary  modes  of  heat  transfer  are  conduction  and  convection.  Conductive 
modes  include  mechanical  thermal  contact  and  solid  thermal  interfaces  between  materi¬ 
als,  such  as  copper,  solder,  or  epoxy  layers.  Convective  modes  include  natural  and  air- 
forced  cooling,  and  forced  liquid  cooling.  Radiation  is  also  a  factor;  however,  it  is  not 
as  significant  as  conduction  and  convection  at  the  temperatures  in  which  printed  circuit 
boards  operate.  (Ref  5] 

A.  CONDUCTION 

Conduction  is  the  transfer  by  molecular  motion  of  heat  between  one  part  of  a  body 
to  another  part  of  the  same  body  or  one  body  and  another  in  physical  contact.  [Ref  1] 

For  the  case  of  conduction,  the  heat  flow  equation  is  the  basis  for  understanding  this 
behavior 


A  7^ 

Q  —  heat  flow  along  the  thermal  gradient, 

2 

A  =  the  area  through  which  the  heat  is  flowing,  m 

W 

k  =  the  thermal  conductivity  of  the  material,  - 
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AL 

A.V 


=  change  in  temperature  per  unit  length, 


Rearranging  Equation  (1)  leads  to 


^th  — 


AT 

Q 


which  is  the  thermal  resistance 


R,h  in 


W 


(2) 


(3) 


Equation  (3)  demonstrates  that  thermal  resistance  will  increase  with  an  increase  in 
path  length  for  heat  flow,  AX,  with  a  decrease  in  area  of  heat  (low,  A,  or  change  in 
conductivity  k  to  a.  lower  value. 

1.  General  Equation  of  Heat  Conduction 
The  general  equation  of  heat  conduction  is 


^  ^  ^  ^  ^  ^  ‘T'  ^  nr* 

C  ,,  C I  ,  ,  C  ,,  C  /  ..  ,  C  ,,  cl  s  .  ^  cl 

-r-  {k  -T—  )  +  —  (A  -T—  )  +  -T-  (A  )  +  <?,  =  pC 
cx  cx  oy  cy  cz  cz  dt 


(4) 


where 


kg 

p  =  density,  — j- 


m 


C  =  specific  heat. 


J 

kg^C 


T  =  temperature,  °C 

x,y,  and  z  =  cartesian  coordmates,  m 


t  =  time,  sec 


k  =  thermal  conductivity, 


W 


m-'‘C 


q  =  internal  heat  generation, 


W 


m 
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Assuming  k,  C  and  p  are  independent  of  temperature,  direction,  and  time,  the  resulting 
equatica  is 


ill  4.  ilL  .  ilL  . -1  =  J_  iL 

CX  cv  cz  ^ 


where 

a  =  thermal  diffusivitv 


k 

’  pC 


_  m 
'  ’  see 


Several  variations  of  the  general  conduction  equation  exist.  Fourier's  equation, 
which  contains  no  heat  sources  is 

,  d^T  1  cT 

CX  CJ  cz 

Another  variation,  known  as  Poisson's  equation,  solves  a  system  in  which  tem¬ 
perature  is  not  time  dependent 

■ff +  -Tf +  +  f  =  (7) 

c.T  cy  cz 

The  last  variant  of  the  general  equation  of  conduction  is  intended  for  a  system 
operating  in  a  steady-state  condition,  and  does  not  contain  any  heat  sources.  This 
equation,  known  as  La  Place's  equation  is 


+  +  =  V'r=0  (8) 

cx^  cy-  cz^ 

2.  Single  Plane  Slab  Modeling 

For  illustrative  purposes,  a  single  plane  slab,  with  T,  and  being  the  face 
temperatures,  is  considered.  If  Equation  (8)  is  limited  to  only  one  coordinate  then 


If  Equation  (9)  is  integrated  twice,  and  boundaries  are  established,  then  the  temperature 
distribution  across  the  slab  can  be  expressed  as  follows 
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(10) 


T=T,-j-{T,-T,) 

If  Equation  (10)  is  substituted  into  Equation  (1),  then  a  solution  for  the  heat  flow  across 
the  single  plane  slab  can  be  obtained 

(11) 

It  is  known  that  Ohm's  law  relates  the  voltage  drop  across  a  resistor  with  the 
current  flow'  through  the  resistor,  V=1R.  Ohm's  law  is  analogous  to  Equation  (1), 
where  the  current  flowing  through  a  resistor  is  equivalent  to  heat  flow,  and  the  electrical 
resistance  is  equivalent  to  the  thermal  resistance  expressed  in  Equation  (2). 

Ohm's  law  can  be  expressed  as 

where  the  analogy  between  Ohm's  law  and  Equation  (1)  is 
Current  /<*Heat  flow  q 
Potential  F«>Temperature  difference  AT 
Resistance  /^oThermal  Resistance  R 

It  can  be  observed  that  for  the  heat  flow  across  a  simple  plane  slab  described 
by  Equation  (11),  the  thermal  resistance  is  described  by  Equations  (2)  and  (3). 

B.  CONVECTION 

Convection  is  defined  as  the  process  by  which  thermal  energy  is  transferred  to  or 
from  a  solid  by  a  fluid  flowing  past  it.  If  the  movement  of  fluid  is  a  result  of  a  temper¬ 
ature  differential,  then  the  process  is  called  free  or  natural  convection.  When  natural 
convection  is  present,  the  movement  of  fluid  can  be  accelerated  by  increasing  the  tem¬ 
perature  differential.  When  a  pressure  differential  is  introduced  to  force  the  movement 
of  the  fluid  by  using  a  pump  or  a  fan,  the  process  is  called  forced  convection.  [Ref 
1.5] 

Newion's  law  of  cooling  states  that  the  heat  flow  through  a  body  is  proportional  to 
the  normal  area  and  the  temperature  difference  between  the  body  and  the  surrounding 
fluid  [Ref  2).  To  make  it  an  equality,  a  proportionality  constant,  h  ,  is  introduced.  The 
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proportionality  constant  is  the  surface  heat  transfer  coefficient.  Newton's  law  of  cool¬ 
ing  can  be  expressed  as 

g  =  hA{To~T^)  (13) 

where 

/i  is  a  proportionality  factor  known  as  the  surface  heat  transfer  coefficient 

Comparing  Newton's  law  of  cooling  with  Fourier's  law  leads  to  an  expression  re¬ 
lating  the  surface  heat  transfer  coefficient  to  thermal  conductivity,  the  surface  fluid 
temperature  differences,  and  the  wall  temperature  grtdient  of  the  fluid 


h  = 


q 

At^T 


(14) 


Consequently,  correlating  heat  transfer  coefficients  must  be  based  on  the  depend¬ 
ence  of  h  on  the  thermal  conductivity  of  the  fluid  and  on  the  ratio  of  the  w’all  temper¬ 
ature  gradient  to  the  temperature  difference.  [Ref  1] 

1.  Electrothermal  Analog 

For  the  case  of  convective  heat  transfer  the  thermal  resistance  is  represented  by 


(15) 


Consequently,  the  total  thermal  resistance  is  now  defined  for  a  single  slab  and  convective 
heat  transfer  on  its  two  faces  as  follows 


(16) 


The  heat  transfer  equation  can  now  be  represented  as 


AT  (T,  -  Fj) 


(17) 


C.  RADIATION 

The  third  mechanism  for  heat  transfer  is  radiation  in  the  form  of  electromagnetic 
waves.  1  he  rate  at  which  a  body  radiates  thermal  energy  is  proportional  to  the  area 
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of  the  body  and  to  the  fourth  power  of  the  absolute  temperature  (Ref  5].  This  result, 
found  empirically  by  Josef  Stefan  in  1S79,  is  written  as 

q  =  eaAT*  (18) 

where 

9  =  power  radiated,  If 
A  =  area, 

e  =  emissivity  of  the  body,  a  value  between  0  and  1 

-8  If’ 

a  =  Stefan-Boltzmann  constant,  5.6703  x  10  — - r 

m  •  k 

When  radiation  falls  on  an  opaque  body,  part  of  the  radiation  is  reflected  and  part 
is  absorbed.  Light-colored  bodies  reflect  most  of  the  radiation,  whereas  dark  bodies 
absorb  most  of  it.  (Ref  5] 

Materials  employed  in  the  manufacture  of  electronic  components  are  classified  as 
gray.  Gray  bodies  arc  diffusely-reflecting  opaque  surfaces.  These  surfaces  reflect  equal 
amounts  of  energy  over  the  thermal  radiation  spectrum  in  all  directions.  (Ref  2] 

1.  Transformation  of  the  General  Radiation  Equation 

The  use  of  the  thermal  radiation  equation  in  analytical  studies  is  a  difficult  task 
due  to  the  fourth  power  relationship  with  temperatures.  The  complexity  of  the  calcu¬ 
lations  involved  with  radiation  dictate  the  need  for  the  aid  of  a  computer  (Ref  2].  The 
general  equation  for  radiation  interchange  is 

q  =  nF,F,A{r,-T^)  (19) 

where,  as  before,  a  is  the  Stefan-Boltzmann  constant  and 

Fa  =  shape  factor  accounting  for  source  and  receiver  arrangement. 

2 

A  =  area,  m 

Fg  =  emissivity  factor  accounting  for  properties  of  the  source  and  receiver. 

=  temperature  of  the  source,  ^  K 
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T,  =  temperature  of  the  receiver,  °K 

It  must  be  observed  that  for  the  radiation  case  absolute  temperature  is  the  norm. 

Radiation  based  heat  transfer  is  represented  by  transforming  the  general  radi¬ 
ation  equation.  Equation  (19>,  into  a  form  compatible  with  Fourier's  law.  Linearization 
of  the  general  radiation  equation  is  the  method  used  to  produce  the  desired  result.  The 
difference  between  the  two  fourth  powers  can  be  reduced  to: 

(7^-7^)  =  (r?  +  r^)(7f-r,') 


=  {T^+Th{T,+  T,)iT,-T,)  (20) 

Inserting  this  into  Equation  (19)  results  in 

q  =  oFMT]  +  Tl){T,  +  T,){T,  -  7))  (21) 

A  radiative  heat  transfer  coefficient  may  be  defined  as 

h,  =  aF,F,if;  +  Tj){T,+  T,)  (22a) 

or 

h,  =  aF,F,{f,  +  T]T,  -f  +  Ti)  (226) 

Substitution  of  6,  into  Equation  (21)  shows  that  radiative  heat  transfer  can  now  be 
treated  similarly  to  convection  at  the  boundary  (Ref  6j.  Thermal  resistance  in  the  case 
of  radiation  heat  transfer  can  now  be  denoted  as 

^  =  (23) 

Now  that  all  three  methods  of  heat  flow  have  been  discussed,  a  graph  depicting 
the  electrothermal  equivalent  is  shown  in  Figure  3.  Here  the  plane  slab  has  both  radi¬ 
ation  and  convection  on  its  wo  faces,  and  the  nonlinearity  of  the  radiation  coefficient, 
h, ,  leads  to  a  more  detailed  analysis  procedure. 

When  considering  the  calculation  of  heat  transfer  by  radiation,  it  is  usuaUy 
necessary  to  approximate  real  mateiial  behavior  with  gray-body  idealization.  There  is 
also  a  great  deal  of  data  required  in  order  to  model  radiation  for  a  real  body  not  only 
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Source:  Reference  8. 
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because  the  surfaces  are  non-black,  but  aho  because  the  non- uniformity  of  the  body's 
surfaces  creates  greater  complexities  in  the  collection  and  processing  of  data.  [Ref  7] 
Basic  problems,  such  as  the  one  shown  in  Figure  3,  defy  a  simple  solution  be¬ 
cause  of  h,.  Therefore,  a  computer-aided  approach  becomes  necessary  in  solving  even 
elementar)’  problems.  .Models  e.xist  for  aiding  the  user  in  solving  radiative  heat  transfer 
problems;  however,  their  development  is  not  within  the  boundaries  of  this  thesis. 
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IV.  FINITE  DIFFERENCE  ANALYSIS 


r-'inite  diilerence  methods  represent  the  most  appropriate  approach  in  determining 
the  temperature  distribution  within  an  electronic  component.  By  using  a  finite  differ¬ 
ence  approach,  complex  problems  involving  diflicult  geometries,  non  linearities,  or 
complicated  boundary  conditions  can  be  successfully  approached.  The  purpose  of  this 
section  is  to  explain  the  fundamental  concepts  behind  finite  difference  analysis  and  its 
use  in  solving  differential  equations.  (Ref  6] 

A.  FUNDAMENTALS  OF  FINITE  DIFFERENCE  ANALYSIS 

The  general  equation  for  heat  transfer  by  conduction.  Equation  (4),  can  be  ex¬ 
pressed  in  a  three-dimensional  variant  as 


ex  ex  cy 


^  '7- 

cv  cz 


cT  .  ,  8T 


(24) 


where  all  the  equation  components  have  been  previously  defined.  In  order  to  formulate 
a  problem  using  Equation  (24),  it  is  first  necessarv'  to  approximate  the  first  and  second 
derivatives. 

I.  First  and  Second  Derivative  Approximation 

The  derivative  of  a  function  at  a  point  can  be  formulated  by  a  finite  difference 
approximation.  It  is  first  necessarx’  to  determine  the  slope  of  a  line  tangent  to  the  point 
of  interest  on  the  plot  of  temperature  as  a  function  of  x,  y,  or  z.  (Figure  4).  In  choosing 
X  as  the  coordinate  system,  it  is  observed  that  the  slope  of  the  line  (m)  is  defined  as 


lim 

AX-.0 


./lA-p  -f  A.r)  -.A-Vq) 
Ax 


(25) 


The  term  of  Equation  (25),  ]if{x^  -I-  Ax)  — /(xo)]/Ajr  is  referred  to  as  the  difference  quo¬ 
tient,  and  is  the  ratio  of  the  change  in  the  value  of  the  function  at  x^  and  -f  Ax  to  the 
change  in  x.  The  limit  of  the  difference  quotient  is  called  the  derivative  of  the  function 
at -Xo.  [Ref  8] 


/’(xq)  =  lim 

A,'c-*0 


/(xq  -f  Ax)  -JIxq) 
Ax 


(26) 
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Figure  4.  First  and  second  derivative  approximation. 
Source:  Reference  8. 


19 


By  letting  J[x)  =  T{x),  the  derivative  for  this  particular  curve  at  the  point  jrj  can  be  de¬ 
rived 


f'{xo)  =  hm 

Ajc-»0 


T{xo  +  Ax)  -  TI-Xq 
Ajc 


As  Ajc  ->  0  Equation  (27)  becomes 
.  dT  T{xo  +  Ax)  -  Tix.) 

- s - 


Thus,  for  Ax  finite,  but  veiy  small 
AT 


/'(^o)  = 


Ax 


(27) 


(28) 


(29) 


The  second  derivative,  using  the  definition  of  a  derivative,  will  be  the  ratio  of  the 
change  in  value  of  the  first  derivative  to  the  change  in  Ax  [Ref.  8].  Again,  looking  at 
Figure  4,  the  second  derivative  is 


/''(xo)  =  lim 

Ax-*0 


fixp  +  Ax)  -/'(xq) 
Ax 


(30) 


/"(^o) 


d^T 

dx^ 


T{xq  ■¥  2Ar)  —  2r(xp  +  Ax)  +  T{xq) 
Ax^ 


(31) 


For  both  Equations  (28)  and  (31),  the  approximations  are  valid  as  long  as  Ax  is  small. 

Taylor's  theorem  states  that  a  function  can  be  approximated  by  a  polynomial 
of  the  form 

Ax)  =Aa)  +f'ia){x  -a)  +  (x  -  «)'  -h  ...  (x  -  o)"  (32) 

l  n. 


where  the  polynomial  is  for  the  function  expanded  about  x  =  a  [Ref.  8].  The  Taylor  se¬ 
ries  for  the  function  T(x)  at  T(xo  -I-  2Ar)  is 


r(xo  -f  2Ax)  =  T(xq  +  Ax)  -f 


dT(X[)  +  Ax) 
dx 


Ax  + 


1  d^T{x,  +  Ax) 

2  dx^ 


Ax"-t- 


(33) 


i 

6 


d  T(xq  -I-  Ax)  3 
—  Ax 


dx^ 


j  d'’  7'(xq  -K  Ax) 

d7^ 


Ax” 
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and  for  the  function  T{x)  at  T{xo)  : 


7'(^o)  =  7'(^o  +  - 


dT{xQ  +  A.v)  I  d  T{Xq  +  Ajc)  , 

-  Ax +  — - - Ax  - 


dx 


dx 


1  d"^  T{xq  +  Ax) 
^  d? 


Ax^  +  ...  + 


1  d'‘  T{xq  +  Ax) 
«!  dx'' 


Ax'' 


(34) 


If  Ax-  is  kept  small,  then  the  terms  above  second  order  become  negligible.  When 
Equations  (33)  and  (34)  are  added,  an  expression  for  the  second  derivative  can  be  ob¬ 
tained 


d  T{xq  -f  Ax)  , 

r(xo  +  2Ax)  +  r(xo)  =  2r(xo  +  Ac)  + - - -  Ax^ 

dx 

d^T{xQ  +  Ax)  T(xq  2Ax)  —  27'(xo  -f  Ax)  -f-  T{xq) 


(35) 


(36) 


As  previously  stated,  the  configuration  to  be  analyzed  can  be  subdivided  into 
small,  finite  subvolumes  considered  to  be  isothermal.  The  centroid  of  each  subvolume 
is  called  a  node,  and  the  node  is  representative  of  the  entire  subvolume.  Nodes  may  are 
connected  to  their  adjacent  nodes  through  thermal  resistances  and  the  nodal  analysis  for 
the  node  temperatures  can  best  be  solved  by  means  of  a  computer  aided  model  [Ref  8|. 

The  first  law  of  thermodynamics  states  that  energy  can  neither  be  created  nor 
destroyed,  but  can  be  transformed  from  one  form  to  another.  An  energy  balance  may, 
therefore,  be  formed  on  the  typical  node  shown  in  Figure  5.  For  the  sake  of  simplicity, 
the  environment  is  composed  by  a  single  node,  node  5  (Figure  6).  Close  observation  of 
node  5  shows  that  it  is  connected  to  nodes  2,  4,  6,  8,  14,  and  100.  If  any  energy  is  di¬ 
rectly  applied  to  or  removed  from  this  node,  it  would  become  part  of  the  equation.  The 
node  equation  for  node  5,  with  rate  of  heat  input  q„  becomes 

+  <76  +  <?8  +  <?14  +  ?100  =  0  (37) 


where  each  one  of  the  q  values  with  a  numerical  subscript  represents  the  rate  of  heat 
flow  from  node  5  to  the  node  indicated  by  the  subscript.  [Ref  1,6,8] 
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»  -I  Reastance 
_ •  Nodes 


Figure  5.  Graphical  representation  of  a  single  node. 
Source:  Reference  6. 


100 


Figure  6.  Node  arrengement  on  t)i>iad  printed  circuit  board. 
Source:  Reference  8. 
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Equation  (37)  does  not  display  the  node  temperatures;  however,  q  —  ATI R  and 
Equation  (37)  can  be  combined  to  produce 


r 'T'  'T'  \ 

L  Ay 

(38a) 

k 4  kAvAz 

0‘-  ^  (n  7-.) 

(38^1) 

(38c) 

kA  .  kAxAz  T- , 

L  Ay 

im 

lA  kAxAv 

q,,-  ^2  (Fs-r.J 

(38c) 

kA  kAxAv 

?,oo-  1  at-  (rj-r.oo) 

(387) 

2 


where  Fioo  is  on  the  top  face  of  node  5  and,  Ax,  Ay,  Az  =  the  distance  between  nodes  in 
the  x,y,  and  z  directions  respectively.  Assuming  Ajt  =  Ay  =  Az  -  1,  then  the  node 
equation  becomes 

k{Ts  -  T,)  +  k{T,  -  T,)  +  k{T,-  T,)  +  (39) 


k{T,  -  T,)  +  kiT,  -  r.J  +  2k{Ts  -  Tioo)  =  ?,• 


or 


-T,-T,  +  IT, +  2r,„o  (40) 

where  the  terms  on  the  right  hand  side  of  Equations  (39)  and  (40)  are  known  values. 
[Ref  1,6,81 

A  node  equation  is  written  for  each  node  in  the  configuration,  and  in  Figure 
8,  which  contains  18  nodes,  18  node  equations  would  be  required  to  obtain  the  18  un¬ 
known  node  temperatures.  A  larger  problem  having  several  hundred  nodes  would  re¬ 
quire  prohibitive  amounts  of  time  to  solve  or  even  enter  data  for  computer-aided  thermal 
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analysis.  That  is  why  a  computer-aided  model  builder  capable  of  constructing  the 
needed  input  data  for  thermal  analysis,  is  so  desirable.  The  thrust  of  this  thesis  is  not 
to  solve  the  equations,  but  to  allow  the  user  to  model  a  printed  circuit  board  to  his  or 
her  specifications  with  minimal  effort,  and  then  produce  a  data  file  compatible  for  use 
in  the  thermal  analvzer. 
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V.  THE  MODEL  BUILDER 


As  previously  indicated,  the  model  builder  currently  in  use  in  the  thermal  analysis 
software  package  requires  time-consuming  manual  data  entiy'.  This  thesis  proposes  that 
a  model  builder,  PCB,  can  be  incorporated  into  the  software  package.  PCB  is  a  menu- 
driven,  user-friendly  program  which  will  assist  the  user  in  developing  a  thermal  model 
of  a  printed  circuit  board  with  pre-defmed  geometries.  PCB  will  generate  a  properly 
formatted  ASCII  output  data  file  for  use  w'ith  the  thermal  analyzer.  The  current  version 
of  the  thermal  analyzer  is  designed  to  handle  a  maximum  of  300  nodes  [Ref  9].  An 
updated  version  capable  of  handling  the  amount  of  nodes  generated  by  the  model  builder 
is  under  development;  however,  abridged  versions  of  the  model  builder  have  been  suc¬ 
cessfully  tested  with  the  current  version  of  the  model  builder.  This  chapter  addresses 
many  aspects  of  the  model  builder,  including  a  terminal  session  demonstrating  the  design 
process  used  in  modeling  a  circuit  board,  and  the  interface  between  the  output  data  file 
and  the  thermal  analyzer.  Also  addressed  is  the  methodology  used  in  validating  the  PCB 
model  builder  output  data  file.  Finally,  limitations  and  possible  improvements  to  the 
model  builder  are  discussed. 

A.  FEATURES 

PCB  is  a  menu  driven  interactive  program  providing  the  user  with  a  variety  of  sys¬ 
tematic  choices  progressing  in  a  logical  sequence  that  eventually  will  lead  to  the  devel¬ 
opment  of  the  model.  Upon  entering  the  program,  the  user  will  be  provided  with  the 
option  of  viewing  an  optional  overview.  First-time  users  should  consider  reading  the 
overview.  The  program  is  designed  to  operate  with  upper-case  letters,  therefore,  it  is 
suggested  that  the  caps-lock  key  is  pressed  before  starting  the  program. 

Following  the  overview  option  the  user  will  be  asked  to  choose  either  SI  or  English 
units.  For  this  choice  as  well  all  others,  the  program  provides  the  user  with  the  op¬ 
portunity  to  correct  any  erroneous  entries. 

The  third  section  queries  the  user  for  the  number  of  copper  layers  desired  in  the 
design.  The  allowable  number  ranges  from  one  to  four  layers.  For  every  choice  of 
number  of  copper  layers,  there  are  four  possible  board  configurations  the  user  can  select. 
If  one  copper  layer  is  selected,  the  program  will  display  the  following  choices  of  aspect 
ratios: 

1.  10  BY  36  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES. 
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2.  10  BY  35  PROVIDES  A  1:3.5  RATIO  WITH  700  NODES. 

3.  12  BY  25  PROVIDES  A  1. 2:2.5  RATIO  WITH  600  NODES. 

4.  15  BY  20  PROVIDES  A  1.5:2.0  RATIO  WITH  600  NODES. 

When  two  copper  layers  are  selected,  the  user  will  be  provided  with  the  following 
choices: 

1.  10  BY  18  PROVIDES  A  1:1.8  RATIO  WITH  720  NODES. 

2.  9  BY  20  PROVIDES  A  1:2.11  RATIO  WITH  720  NODES. 

3.  8  BY  20  PROVIDES  A  2:5  RATIO  WITH  640  NODES. 

4.  10  BY  15  PROVIDES  A  1:1.5  RATIO  WITH  600  NODES. 

If  three  copper  layers  are  selected,  the  user  will  have  the  following  choices; 

1.  5  BY  24  PROVIDES  A  1:4.8  RATIO  WITH  720  NODES. 

2.  10  BY  10  PROVIDES  A  1:1  RATIO  WITH  600  NODES. 

3.  12  BY  10  PROVIDES  A  1.2:1  RATIO  WITH  720  NODES. 

4.  8  BY  15  PROVIDES  A  1:1.875  RATIO  WITH  720  NODES. 

By  selecting  four  copper  layers,  the  program  will  provide  the  user  the  following 
options: 

1.  5  BY  18  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES. 

2.  9  BY  10  PROVIDES  A  1:1.1 1  RATIO  WITH  720  NODES. 

3.  8  BY  10  PROVIDES  A  1:1.25  RATIO  WITH  640  NODES. 

4.  13  BY  6  PROVIDES  A  1:1  RATIO  WITH  624  NODES. 

After  selecting  the  desired  aspect  ratio,  the  user  must  specify  the  epoxy  character¬ 
istics.  The  program  will  prompt  the  user  to  enter  the  epoxy  layers  length,  width, 
thickness  and  thermal  conductivity.  An  option  available  to  the  user  will  be  the  choice 
of  a  default  epoxy  thickness  value,  0.0625  in.,  (0.15875  cm.). 

Once  the  epoxy  characteristics  have  been  entered,  the  program  will  then  display  the 
corresponding  length  and  thickness  for  the  copper  layers  selected.  The  length  and  width 
for  the  copper  layers  will  be  the  same  as  those  of  the  epoxy  layers.  The  user  will  then 
be  provided  with  the  option  of  entering  the  thickness  of  each  selected  copper  layer  either 
by  length  or  by  weight.  If  input  by  length  is  selected,  the  units  should  be  consistent  with 
the  system  of  measurement  selected;  however,  if  entry  of  thickness  by  weight  is  selected, 
then  the  entries  should  be  in  ounces.  A  copper  layer  weighing  one  oz.  is  equivalent  to 
a  copper  layer  0.0014  in.  thick  with  a  surface  area  of  144  in^  .  The  program  adjusts 
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thelayer  thickness  to  the  surface  area  provided  automatically.  Entry  of  thickness  by 
weight  is  a  standard  practice  used  in  industry.  |Ref.  2] 

After  entering  the  copper  layer  thickness  the  program  will  query  the  user  for  the 
thermal  conductivity  of  copper.  After  this  entry,  the  program  will  redisplay  all  infor¬ 
mation  pertaining  to  the  copper  layers  and  will  require  the  user  to  confirm  all  entries  in 
order  to  proceed. 

The  following  section  requests  the  input  of  initial  and  ambient  temperatures.  PCB 
has  six  ambient  temperatures. 

After  the  temperatures  have  been  entered  the  program  requires  the  user  to  input  the 
external  heat  sources.  PCB  provides  for  heat  input  into  the  upper  copper  layer.  There 
are  four  methods  of  external  heat  input  from  which  to  choose.  The  first  alternative  al¬ 
lows  for  a  total  rate  of  heat  applied  to  the  upper  surface.  An  entry  for  this  choice  would 
be  divided  by  the  number  of  nodes  and  distributed  uniformly.  The  second  choice  pro¬ 
vides  for  the  entry  of  average  heat  per  unit  area.  The  third  alternative  gives  the  user  the 
ability  to  enter  heat  in  specifically  designated  nodes.  The  last  option  provides  for  no 
heat  input  and  was  developed  to  test  compatibility  with  the  thermal  analyzer. 

The  following  selection  is  the  percent  copper  coverage  for  each  copper  layer.  The 
program  will  prompt  the  user  to  enter  the  percentage  of  copper  coverage  for  each  layer 
in  the  printed  circuit  board. 

Finally,  the  program  requires  a  name  and  a  title  for  the  file  that  will  be  created.  This 
file  will  be  in  the  proper  format  for  use  with  the  thermal  analyzer. 

B.  THE  THERMAL  ANALYZER  INPUT  DATA  FILE 

The  model  builder  generates  an  ASCII  data  file  from  the  physical  characteristics  of 
the  printed  circuit  board  provided  by  the  user.  In  order  for  the  input  data  file  to  be  ac¬ 
ceptable  to  the  thermal  analyzer,  it  must  be  in  a  specific  format  which  is  compatible  to 
the  thermal  analyzer  [Ref  9].  Because  each  value  and  position  of  the  output  data  file 
has  a  meaning  to  the  thermal  analyzer,  and  is  not  readily  identifiable  to  the  user,  it  is 
beneficial  to  describe  each  line  and  data  set  of  the  output  data  file  and  their  relationship 
to  PCB.  Figure  7  shows  a  partial  output  data  file. 

Line  one  is  the  title  line.  It  may  be  left  blank  or  may  contain  up  to  79  alphanumeric 
characters.  The  user-selected  title  appears  at  the  top  of  the  data  file.  [Ref.  6] 

Line  two  is  the  problem  data  line.  It  has  nine  entries  of  which  two  are  under  user 
control,  the  number  of  nodes  under  consideration  and  the  unit  type.  One  entr>,  the 
number  of  constant  temperatures,  is  preset  at  six  for  this  specific  model.  The  remaining 
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Figure  7.  Sample  PCB  partial  output  data  file. 


entries  have  applications  to  models  associated  with  heaters,  unique  exponents,  secondarx’ 
heat  input,  temperature  coefficients  and  curves,  and  nodes  controlling  fast  heat.  These 
entries  are  not  applicable  to  this  model  and  are  preset  to  zero.  [Ref  6] 

Line  three  places  a  zero  at  three  points  and  is  beyond  the  user's  control.  Therefore, 
no  further  discussion  is  required.  [Ref  6] 

Line  four  is  the  problem  capability  line.  This  line  defines  the  maximum  values  for 
the  entries  in  line  two.  The  first  entiy  is  750,  which  is  the  number  of  nodes  for  which 
the  analysis  is  dimensioned.  The  number  750  is  significant  because  the  first  constant 
temperature  will  be  assigned  to  node  number  751.  The  second  entry  is  50  which  re¬ 
presents  the  largest  possible  number  of  constant  temperatures  in  accordance  with  the 
analyzer  dimension  statement.  The  third  entry  is  set  to  6  and  does  not  change.  This 
entry  is  related  to  heaters  and  is  not  applicable  to  the  model.  The  balance  of  the  entries 
in  line  four  represent  a  listing  of  data  sets  that  are  required  for  the  particular  analysis 
at  hand.  PCB  uses  three  data  sets  that  will  be  discussed  in  what  follows.  [Ref  6] 

The  fifth  line  contains  five  values  that  relate  to  the  accuracy  level  that  the  thermal 
analyzer  will  achieve.  These  entries  are  preset.  The  first  value  provides  the  level  of  ac¬ 
curacy  between  iterations.  The  accuracy  level  number  is  critical  because  too  small  a 
tolerance  will  cause  the  computer  to  run  for  excessive  amounts  of  time,  and  too  large  a 
number  will  provide  inaccurate  results.  The  second  value  is  the  damping  factor  used 
between  iterations  in  order  to  prevent  temperature  oscillations  between  iterations.  The 
third  number  provides  the  maximum  number  of  iterations.  If  erroneous  data  is  entered, 
the  computer  will  not  run  for  excessive  amounts  of  time.  The  fourth  value  is  the  con¬ 
vergence  factor  which  adjusts  the  damping  in  order  to  close  the  critical  value.  The  fifth 
entrv'  is  the  initial  temperature  at  which  the  iterative  process  begins.  This  value  is  sup¬ 
plied  by  the  user.  [Ref  6] 

Line  six  contains  the  temperature  dependent  coefficients,  and  is  not  used  in  this 
model.  Line  seven  contains  up  to  50  constant  temperature  inputs.  This  model  has  six 
ambient  temperatures  entered  by  the  user.  [Ref  6] 

The  following  lines  contain  all  pertinent  information  concerning  the  n-node 
equations.  Each  node  requires  two  lines  of  data.  Even-numbered  lines  are  used  for 
specifying  the  nodes  that  interact  with  the  node  in  question  and  the  modes  by  which  this 
interaction  takes  place.  For  example,  a  line  of  the  form; 

6  7551  21  7521  111  7511  2411 
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is  for  node  number  one  of  aspect  ratio  selection  one.  The  first  entr>’  indicates  the  num¬ 
ber  of  connections  to  that  node.  The  second  entr\'  is  indicating  that  node  755,  an  am¬ 
bient  temperature,  is  connected  to  the  node  in  question  and  the  1  indicates  that  the 
connection  is  conductive.  The  same  procedure  applies  to  the  rest  of  the  values  on  the 
line.  If  the  entry  is  9991  then  the  heat  input  is  external.  [Ref  6] 

Other  data  sets  that  appear  relate  to  unique  exponents,  secondary  heat,  and 
temperature-dependent  heat  input  curves.  These  data  sets  are  not  used  in  this  model. 
[Ref  6] 

C.  PCB  MODEL  BUILDER  SAMPLE  PROBLEM 

The  following  terminal  session  is  a  typical  example  of  the  printed  circuit  board 
model  builder,  PCB,  and  its  capabilities.  Figure  8  depicts  the  printed  circuit  board 
modeled  in  the  terminal  session. 

1.  Printed  Circuit  Board  Specifications 

1.  Unit  system:  British. 

2.  Copper  layers  in  printed  circuit  board:  three. 

3.  Aspect  ratio  desired:  12  by  10,  with  720  nodes. 

4.  Epoxy  layer  length:  7.2  in. 

5.  Epoxy  layer  width:  4.2  in. 

6.  Epoxy  layer  thickness  to  default  value:  Yes  (.0625  in.) 

I.  Epoxy  layer  thermal  conductivity:  0.087  Btu  hr/'F. 

8.  Specify  copper  layer  thickness  by  weight  or  length:  length. 

9.  Thickness  for  the  three  copper  layers,  respectively:  .20,  .10,  .11  in. 

10.  Copper  thermal  conductivity:  243.000  Btu'hr  "f. 

II.  Initial  board  temperature:  87  T. 

12.  Upper  surface  ambient  temperature:  86  'f. 

13.  Lower  surface  ambient  temperature:  85  °F. 

14.  Right  surface  ambient  temperature:  84  °F. 

15.  Left  surface  ambient  temperature:  86  °F. 

16.  Front  surface  ambient  temperature:  86  °f. 

17.  Rear  surface  ambient  temperature:  86  °F. 

18.  Heat  input  to  the  upper  copper  layer:  300.23  Btu  hr  °F. 

19.  Percent  copper  coverage  for  copper  layers,  respectively:  87.3.  76.2,  65.7  %. 

20.  Name  of  output  data  file:  LAYERS. 
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21.  Header  of  the  output  data  file:  THIS  IS  THE  OUTPUT  OF  THE  THERMAL 
A^’ALYZER. 

2.  Terminal  Session 

Using  the  board  specifications  and  Figure  8,  the  terminal  session  for  the  PCB 
program  is  as  follows; 


THIS  PROGRAM  WAS  WRITTEN  TO  INTEGRATE  WITH  EXISTING 
THERMAL  ANALYSIS  SOFTWARE  AND  TO  REDUCE  THE  AMOUNT 
OF  TIME  REQUIRED  FOR  DATA  ENTRY. 

WOULD  YOU  LIKE  AND  C'/ERVIEW  OF  THE  PROGRAM  PRIOR  TO 
BEGINNING?  ENTER  Y  FOR  YES  AND  N  FOR  NO:  N 


************************************************* 

************************************************* 


**  ** 

**  PLEASE  SELECT  UPPER  CASE  LETTERS  ** 

**  PRIOR  TO  BEGINNING  ** 

**  ** 


************************************************* 

************************************************* 


PRIOR  TO  ENTERING  DATA  INTO  THIS  PROGRAM  ENSURE  THAT 
YOU  HAVE  A  DRAWING  OF  YOUR  DESIGN  AND  ALL  PERTINENT 
DATA. 


PRESS  <ENTER>  TO  CONTINUE 
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THIS  PROGRAM  IS  CAPABLE  OF  OPERATIONS  IN  EITHER  SI  OR 
ENGLISH  UNITS.  AFTER  THE  SELECTION  OF  THE  UNITS,  ALL 
ENTRIES  MUST  BE  COMPATIBLE.  PLEASE  MAKE  YOUR  SELECTION. 

S  FOR  SI  NOTATION 
E  FOR  ENGLISH  NOTATION;  E 


YOU  HAVE  SELECTED  ENGLISH  NOTATION. 


IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FOR  NO:  Y 


THE  PCB  CAN  HAVE  UP  TO  FOUR  COPPER  LAYERS 

HOW  MANY  COPPER  LAYERS  DO  YOU  DESIRE? 

PICK  A  NUMBER  1  TO  4:  3 


YOU  SELECTED  3  COPPER  LAYER (S)  FOR  THE  PCB 

IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FOR  NO;  Y 
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YOU  SELECTED  3  COPPER  LAYERS,  GIVING  YOU  THE  FOLLOWING  ALTERNATIVES 


1.  5  BY  24  PROVIDES  A  1:4.8  RATIO  WITH  720  NODES. 

2.  10  BY  10  PROVIDES  A  1:1  RATIO  WITH  600  NODES. 

3.  12  BY  10  PROVIDES  A  1.2:1  RATIO  WITH  720  NODES. 

4.  8  BY  15  PROVIDES  A  1:1.875  RATIO  WITH  720  NODES. 

PLEASE  SELECT  A  NUMBER  1  THROUGH  4:  3 


YOU  SELECTED  NUMBER  3  OF  THE  FOLLOWING  4  ALTERNATIVES: 


1.  5  BY  24  PROVIDES  A  1:4.8  RATIO  WITH  720  NODES. 

2.  10  BY  10  PROVIDES  A  1:1  RATIO  WITH  600  NODES. 

3.  12  BY  10  PROVIDES  A  1.2:1  RATIO  WITH  720  NODES. 

4.  8  BY  15  PROVIDES  A  1:1.875  RATIO  WITH  720  NODES. 


IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FOR  NO:  Y 
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******************** *******CPOXY  CHARACTERISTICS************************ 
************************************************************************ 

ALL  ENTRIES  ARE  IN  ENGLISH  NOTATION 


ENTER  EPOXY  LAYER  LENGTH  (in):  7.2 
ENTER  EPOXY  LAYER  WIDTH  (in):  4.2 

THE  EPOXY  LAYER  THICKNESS  WILL  DEFAULT  TO  0.0625  in  (0.15875  cm) 
DO  YOU  WANT  TO  CHANGE  THE  EPOXY  THICKNESS?  (Y  OR  N) :  N 

ENTER  EPOXY  LAYER  THERMAL  CONDUCTIVITY  (Btu/hr/F  ):  0.087 


YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  EPOXY  LAYER: 


1.) 

LENGTH: 

7.2000 

in 

2.) 

WIDTH: 

4.2000 

in 

3.) 

THICKNESS: 

.0625 

in 

4.) 

k: 

.0870 

Btu/hr/F 

DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  YES  AND  N  FOR  NO:  N 
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*****************COPPER  LAYER  CHARACTERISTICS*************************** 

ALL  ENTRIES  ARE  IN  ENGLISH  NOTATION 


THE  COPPER  LAYER  LENGTH  IS  THE  SAME  AS  THE  EPOXY  LAYER;  7.2000  in 

THE  COPPER  LAYER  WIDTH  IS  THE  SAME  AS  THE  EPOXY  LAYER:  4.2000  in 

YOU  SELECTED  3  COPPER  LAYER (S)  FOR  THE  PCB. 

YOU  NOW  WILL  BE  ASKED  TO  ENTER  THE  COPPER  LAYER  THICKNESS 

DO  YOU  WANT  TO  SPECIFY  THICKNESS  BY  LENGTH  OR  WEIGHT? 

ENTER  L  IF  YOU  WANT  TO  ENTER  LENGTH,  W  IF  YOU  WANT  TO  ENTER  WEIGHT:  L 


ENTER  THE  THICKNESS  FOR  LAYER  1  (in):  .2 

ENTER  THE  THICKNESS  FOR  LAYER  2  (in):  .1 

ENTER  THE  THICKNESS  FOR  LAYER  3  (in):  .11 

ENTER  COPPER  LAYER  THERMAL  CONDUCTIVITY  (Btu/hr/F  ) :  243 


YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  COPPER  LAYER (S) . 


1.) 

LENGTH : 

2.) 

WIDTH; 

3. ) 

THICKNESS 

LAYER 

1 

THICKNESS 

LAYER 

2 

THICKNESS 

LAYER 

3 

4.) 

k: 

7.2000 

4.2000 
.2000 
.1000 
.1100 

243.0000 


in 

in 

in 

in 

in 

Btu/hr/F 
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************************************************************************ 
******************ambient  temperature  inputs**************************** 
************************************************************************ 


ENTER  THE  INITIAL  BOARD  TEMPERATURE  (F) :  87 
ENTER  THE  UPPER  SURFACE  AMBIENT  TEMPERATURE  (F) :  86 
ENTER  THE  LOWER  SURFACE  AMBIENT  TEMPERATURE  (F) :  85 
ENTER  THE  RIGHT  SURFACE  AMBIENT  TEMPERATURE  (F) :  84 
ENTER  THE  LEFT  SURFACE  AMBIENT  TEMPERATURE  (F) :  86 
ENTER  THE  FRONT  SURFACE  AMBIENT  TEMPERATURE  (F) :  86 
ENTER  THE  REAR  SURFACE  AMBIENT  TEMPERATURE  (F) :  86 


YOU  HAVE  MADE  THE  FOLLOWING  AMBIENT  TEMPERATURE  ENTRIES: 


1.  INITIAL  BOARD  TEMPERATURE:  87.000  F 

2.  UPPER  AMBIENT  TEMPERATURE:  86.000  F 

3.  LOWER  AMBIENT  TEMPERATURE:  85.000  F 

4.  RIGHT  AMBIENT  TEMPERATURE:  84.000  F 

5.  LEFT  AMBIENT  TEMPERATURE:  86.000  F 

6.  FRONT  AMBIENT  TEMPERATURE:  86.000  F 

7.  REAR  AMBIENT  TEMPERATURE:  86.000  F 


DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  YES  AND  N  FOR  NO:  N 


DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  YES  AND  N  FOR  NO:  N 
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************************************************************************ 
**************** **************H£AT  inputs******************************* 
************************************************************************ 


HEAT  INPUT  TO  THE  PCB  OCCURS  ONLY  ON  THE  UPPER  COPPER 
LAYER.  HEAT  INPUT  IS  ACCOMPLISHED  BY  ONE  OF  THE 
FOLLOWING  METHODS: 

1.  ENTER  AS  TOTAL  HEAT  APPLIED  TO  THE  PCB 

2.  ENTER  AS  AVERAGE  HEAT  PER  UNIT  AREA 

3.  ENTER  HEAT  NODE  BY  NODE 

4.  NO  HEAT  INPUT 

PLEASE  SELECT  NUMBER  1  THROUGH  4:  1 


YOU  HAVE  SELECTED  NUMER  1  OF  FOUR  ALTERNATIVES. 
IS  THIS  THE  DESIRED  SELECTION?  (Y  OR  N) :  Y 


YOU  HAVE  SELECTED  TO  INPUT  HEAT  AS  A  TOTAL  HEAT 
APPLIED  TO  THE  SURFACE. 

ENTER  TOTAL  HEAT  APPLIED  TO  THE  SURFACE  (Btu/hr) :  300.23 
IS  THIS  THE  CORRECT  ENTRY?  (Y  OR  N) :  Y 
TOTAL  HEAT  PER  NODE  IS;  2.5019  Btu/hr 

<PRESS  ENTER  TO  CONTINUE> 


ENTER 

THE 

PERCENT 

COPPER 

COVERAGE 

(i.e. 

67) 

FOR 

LAYER 

1: 

87.3 

ENTER 

THE 

PERCENT 

COPPER 

COVERAGE 

(i.e. 

67) 

FOR 

LAYER 

2: 

76.2 

ENTER 

THE 

PERCENT 

COPPER 

COVERAGE 

(i.e. 

67) 

FOR 

LAYER 

3: 

65.7 
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YOU  HAVE  SELECTED;  87.30%  COVERAGE,  LAYER  1 

76.20%  COVERAGE,  LAYER  2 
65.70%  COVERAGE,  LAYER  3. 

IS  THIS  YOUR  DESIRED  ENTRY?  (Y  OR  N) ;  y 


YOU  SELECTED  LAYERS  FOR  YOUR  DATA  FILE  NAME 


IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FOR  NO:  Y 


THIS  PROGRAM  CREATES  AN  OUTPUT  DATA  FILE  FOR  ENTRY  INTO  THE 
EXISTING  THERMAL  ANALYZER,  FURTHERMORE,  THIS  PROGRAM  DOES 
NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE.  THEREFORE 
THE  USER  WILL  NAME  THE  DATA  FILE  FOR  EACH  RUN  OF  THIS 
PROGRAM.  THE  FILE  NAME  IS  LIMITED  TO  SIX  CHARACTERS,  AND 
SHOULD  NOT  HAVE  ANY  SPACES. 


PLEASE  ENTER  THE  DESIRED  DATA  FILE  NAME;  LAYERS 


YOU  SELECi’ED  LAYERS  FOR  YOUR  DATA  FILE  NAME 

IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FOR  NO;  Y 
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ENTER  THE  DESIRED  TITLE  TO  BE  PLACED  ON  LINE 
NUMBER  ONE  OF  THE  OUTPUT  DATA  FILE: 


THIS  IS  THE  OUTPUT  OF  THE  THERMAL  ANALYZER 


DO  YOU  WISH  TO  CHANGE  THE  TITLE  OF  YOUR  OUTPUT  DATA  FILE? 
ENTER  Y  FOR  YES  AND  N  FOR  NO:  N 


THE  OUTPUT  DATA  HAS  BEEN  PLACED  IN  A  FILE  NAMED  LAYERS 


<PRESS  ENTER  TO  CONTINUE> 


3.  The  output  data  file 

Figure  9  shows  a  partial  listing  of  the  model  builder  output  computed  in  the 
terminal  session.  A  detailed  description  of  the  output  data  file  is  provided  in  section  B 
of  this  chapter. 

4.  Thermal  analyzer  output 

The  file  obtained  in  the  model  builder  produces  an  output  for  720  nodes.  In 
order  to  use  this  data  file  with  the  thermal  analyzer  the  data  file  must  first  be  shortened. 
In  this  particular  case,  the  output  file,  LAYERS,  shortened  to  eight  nodes  was  entered 
to  the  thermal  analysis  program,  TASS.  Three  files  are  produced  by  the  thermal  analysis 
program,  TASS.  The  first  two  output  files  consist  of  detailed  and  summarized  thermal 
analysis  results,  respectively.  The  third  file  is  an  error  message  summary  in  the  event 
there  is  an  error  in  the  input  data  file.  Figure  10  shows  the  detailed  output  file  of  the 
eight  node  thermal  analysis. 

D.  PCB  MODEL  BUILDER  VALIDATION 

The  output  of  the  model  builder  was  validated  by  following  a  three-step  process. 
The  first  step  in  the  validation  process  concerned  format.  In  order  to  proceed  with  any 
other  validation  steps  it  was  first  necessary  to  ensure  that  all  data  fields  were  in  accord¬ 
ance  with  the  specifications  set  forth  in  the  user's  manual  of  the  thermal  analyzer  (Ref 
9].  The  thermal  analyzer  input  data  file  consists  of  five  lines  and  as  many  as  seven  input 
data  sets.  The  thermal  analyzer's  manual  model  builder  program,  THANSS,  was  in¬ 
strumental  in  verifying  the  accuracy  of  the  model  builder  output  data  file  format. 

The  second  step  entailed  content.  The  output  of  the  model  builder  was  checked  for 
accuracy  and  completeness.  Thermal  resistances  in  the  output  uata  files  were  checked 
and  compared  with  manually  calculated  benchmark  models.  A  benchmark  model  was 
developed  for  each  copper  layer  configuration.  Other  information  such  as  node  numbers 
and  node  relationships  were  also  compared  with  the  benchmark  model  and  verified  for 
accuracy.  To  aid  in  the  verification  process,  the  manual  model  builder  program, 
THANSS,  part  of  the  thermal  analyzer  software  was  used  to  verify  node  relationships 
[Ref  9).  Other  items  checked  for  correctness  were  specification  codes  relating  to  units 
(British  or  SI),  temperature  inputs,  and  thermal  analyzer  default  values. 

The  third  and  final  step  in  the  validation  process  required  running  the  thermal  ana¬ 
lyzer  program  with  the  PCB  output  data  file  as  the  program  input.  Abridged  versions 
of  the  model  builder  were  successfully  executed  with  the  current  version  of  the  thermal 
analyzer  (maximum  node  capability  is  300  nodes).  The  results  of  the  thermal  analysis 
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THIS  IS  THE  1 

OUTPUT  OF  THE 

THERMAL  ANALYZER 

720  6 

0 

0  0 

0 

0  0 

1 

0  0 

0 

750  50 

6 

2  4 

6 

0  0 

0 

.0500000  . 

6666670 

12 

.8000000 

87.00000 

S6.000 

86.000 

86.000  84. 

000  86. 

000  85. 

.000 

7  7551 

21 

7521 

131 

7511 

1211 

242.445 

67.280 

41.249 

20.625 

.701 

.701 

2.50 

7  11 

31 

7521 

141 

7511 

1221 

B7.2B0 

87.280 

41.249 

20.625 

.701 

.701 

2.50 

7  21 

41 

7521 

151 

7511 

1231 

87.280 

87.280 

41.249 

20.625 

.701 
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171 

7511 

1251 
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1261 

87.280 

67.280 

41.249 
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7511 

1301 
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71 
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211 

81 
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.701 
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91 
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1411 
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20.625 

.701 

.701 

2.50 

7  211 

231 

101 
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7511 

1421 

87.280 

87.280 

20.625 

20.625 

.701 

.701 

2.50 

7  221 

241 

111 
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7511 

1431 

87.280 

67.280 

20.625 

20.625 

.701 

.701 

2.50 

7  231 

251 

121 
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7511 

1441 

87.280 

87.280 

20.625 

20.625 

.701 

.701 

2.50 

7  241 

261 

131 

371 

7511 

1451 

87.280 

87.280 

20.625 

20.625 

.701 

.701 

2.50 

7  251 

271 

141 

381 

7511 

1461 

87.280 

87.280 

20.625 

20.625 

.701 

.701 

2.50 

7  261 

281 

151 

391 

7511 

1471 

Figure  9.  Partial  output  data  fiie  of  PCB  terminal  session. 


were  compared  with  the  results  obtained  using  the  analyzer's  manual  model  builder, 
THANSS,  and  checked  for  discrepancies. 
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Figure  10.  Output  data  file  of  TASS  thermal  tntlyxcr. 
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E.  MODEL  BUILDER  LIMITATIONS  AND  POSSIBLE  IMPROVEMENIS 
Although  the  model  builder  provides  the  user  with  a  significant  advantage  in  terms 
of  time  and  efTori  necessary  to  model  a  printed  circuit  board,  there  are  at  present  certain 
limitations  to  the  program.  There  is  also  significant  potential  for  improvements  that 
could  enhance  the  capabilities  of  the  current  version  of  the  model  builder. 

1.  PCB  Model  Builder  Limitations 

1.  The  model  builder  is  restricted  to  a  maximum  ol'four  copper  layers  due  to  con¬ 
ventional  personal  computer  operating  system  limitations,  Current  versions  of 
MS-DOS  ®  do  not  allow  the  user  to  access  extended  memory.  Access  to  extended 
memory  would  allow  for  a  larger  copper  layer  capability  for  the  model. 

2.  The  model  builder  will  not  run  when  using  .MS-DOS  ®,  version  4.01,  When  run¬ 
ning  PCI!  on  a  computer  operating  with  MS-DOS  ®,  version  4,01,  the  computer 
will  display  an  "insulficicnt  memoiy"  me.ssagc.  Writing  more  compact  code  for  the 
model  builder  could  overcome  this  problem.  The  program  will  run  on  any  other 
version  of  MS-DOS  ®. 

3.  The  model  builder  doc.s  not  account  for  convective  and  radiative  heat  transfer, 

4.  The  model  builder  only  accepts  heat  inputs  from  the  top  surface  of  the  printed 
circuit  board. 

2.  Pos.sible  Improvements  to  tiie  Thermal  Model  Builder. 

1.  The  present  version  of  the  model  builder  requires  extensive  amounts  of  FORTILAN 
code.  An  updated  model  builder  could  be  designed  using  a  more  compact  and  ef¬ 
ficient  higher  level  language,  such  as  C  +  +,  Smaller,  more  clTicicnt  source  code 
would  allow  for  incorporating  more  copper  layers  into  the  model  builder. 

2.  The  menu  driven  queries  used  in  the  thermal  analyzer  could  be  improved  by  in¬ 
corporating  as.sembly  language  routines  enabling  the  use  of  a  mouse.  Attempts 
were  m.idc  to  introduce  the  use  of  mouse  driven  menus;  however,  memory  handling 
limitations  in  the  Microsoft  ®  linker,  version  3.55,  prevented  the  use  of  a  mouse. 

3.  Allowing  for  the  injection  of  heat  not  only  from  the  top  surface  of  the  printed  cir¬ 
cuit  board,  but  also  from  other  external  surfaces  (to  include  the  heat  being  trans¬ 
ferred  through  the  wedges  supporting  the  printed  circuit  board),  is  another 
potential  improvement. 

4.  Adding  graphics  capabilities  to  the  model  builder  would  enable  the  user  to  model 
the  printed  circuit  board  without  the  need  for  preliminary  drawings  or  sketches; 
however,  adding  this  feature  would  cause  the  model  builder  executable  code  to  be¬ 
come  excessively  large. 
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VI.  CONCLUSIONS 


The  purpose  of  this  thesis  is  to  develop  a  printed  circuit  board  thermal  analysis 
model  builder  tluat  will  c/Icctivcly  interact  with  thermal  analysis  software  used  by  the 
Naval  Postgraduate  School.  [Ref.  9J 

The  current  version  of  the  model  builder  enables  the  user  to  choose  from  sixteen 
possible  printed  circuit  board  configurations,  ranging  from  one  to  four  copper  layers. 
Other  features  of  the  model  builder  include  the  ability  to  set  the  physical  characteristics 
of  the  board  including  dimensions,  thermal  conductivities,  percent  copper  coverage, 
board  temperatures,  and  heat  input  modes. 

The  main  goal  of  the  model  builder  is  to  enable  the  user  to  model  a  printed  circuit 
board  with  minimum  elTort,  and  produce  a  thermal  analysis  input  data  file.  Using  the 
model  builder,  PCB,  relieves  the  user  of  the  tedious,  time-consuming  manual  data  entry 
required  of  the  thermal  analyzer's  current  model  builder,  THANSS.  [Ref.  9] 

The  expansion  possibilities  for  the  model  builder  are  significant.  Adding  more  cop¬ 
per  layers  and  nodes  would  enable  the  user  to  model  more  sophisticated  printed  circuit 
boards.  Other  potential  features  include  the  addition  of  graphics  capabilities  and  the 
introduction  of  user-friendly  peripherals,  such  as  mouse  drivers  or  digitazing  pads.  As 
the  level  of  complexity  of  future  model  builders  increases,  it  will  become  necessary'  to 
import  the  source  code  from  FORTRAN  compilers  limited  to  640  Kilobytes  of  memory', 
to  32  bit  FORTRAN  compilers  allowing  use  of  all  RAM.  The  source  code  for  PCB 
contains  over  5000  lines  of  code,  effectively  reaching  the  memory  limits  of  both  The 
compiler  and  personal  computer. 

The  possible  applications  for  the  model  builder  are  substantial.  As  the  level  of 
complexity  of  printed  circuit  boards  increases,  there  is  a  real  need  for  a  tool  that  allows 
the  designer  to  efficiently  perform  thermal  analysis  of  printed  circuit  boards  during  the 
design  process.  The  model  builder  employed  must  allow  the  designer  the  necessary 
flexibility  needed,  to  model  a  board  that  meets  the  desired  design  criteria.  The  PCB 
model  builder  is  the  first  step  in  creating  a  tool  that  enables  the  designer  to  effectively 
build  a  model  that  meets  specified  requirements. 
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APPENDIX  PROGRAM  LISTING  FOR  PCB  MODEL  BUILDER 


$ LARGE 


TITLE; 

AUTHOR: 

DATE: 

COMPILER; 

LINKER; 


MODEL  BUILDER  --  MAIN  PROGRAM 
LT  STEVE  GLASER 
09  JUL  1991 

MICROSOFT  VERSION  4.  01 
MICROSOFT  VERSION  3.55 


DEFINE  REAL  VARIABLES 


REAL  EL.EW ,DELE ,UKE ,CL,CW,T1L1 ,T2L1 ,T2L2 ,T3L1 ,T3L2 ,T3L3 ,T4L1 ,T4L2 , 
+T4L3 , T4L4 . TIW 1 , T2W 1 , T2W2 , T3W1 ,T3W2 , T3W3 , T4W 1 , T4W2 , T4W3 , T4W4 , ARE A , U 
+N , VOL 1 , V0L2 , V0L3 , V0L4 , UKC , I BT , UPRT , LWRT , RT , LT , FT , BT , THE AT , THPN , AHE 
+AT , NHEAT , DELX , DELY 

COMMON  I BT, UPRT, LWRT, RT.LT, FT, BT,UT, ANSWER 
INTEGER  NWIDE,NDEEP,NPL,NUM,NUMA,CH,H 


DEFINE  ALL  ONE  CHARACTER  VARIABLES 


CHARACTER* 1  SELECT , ANS , ANSL , OVR , ANSWER , ANSN , ANSA , LAY 1 , LAY2 , LAY3 , LA 
+Y4 , ANSE , ANSEL , ANSEW , ANSET , ANSEK , ANSQ , WORL , ANSE 1 , ANSTl , ANSK 1 , ANSQl , 
+SPEVAL , ANSTEM , ANSTEB , ANSTEU , ANSTEL , ANSTER , ANSLFT , ANSFRT , ANSBCK 


CHARACTER  VARIABLES  OF  MORE  THAN  ONE  POSITION 


CHARACTER  UK*11 ,UT*1 ,UH*13 ,UAH*6 ,UL*2 ,UW*2 
DEFINE  MATRICES 

REAL  HEAT(100,100),COEF(1000,12) 

INTEGER  IH( 1000) , JH( 1000) ,NCON( 1000,12) 
PROVIDE  THE  USER  WITH  AN  INTRODUCTORY  STATEMENT. 

88888  CALL  CLS 
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WRITE(*,7001) 

7001  FORMATC ///////,' 

+'  THIS  PROGRAM  WAS  WRITTEN  TO  INTEGRATE  WITH  EXISTING 
+'  THERMAL  ANALYSIS  SOFTWARE  AND  TO  REDUCE  THE  AMOUNT 
+  ’  OF  TIME  REQUIRED  FOR  DATA  ENTRY. 

+'  WOULD  YOU  LIKE  AND  OVERVIEW  OF  THE  PROGRAM  PRIOR  TO 

+'  BEGINNING?  ENTER  Y  FOR  YES  AND  N  FOR  NO:  ',2X,  ) 

READ(*,7002)  OVR 

7002  FORMAT(Al) 


VARIABLE,  CONSTANT,  AND  STRING  DEFINITION 
PHYSICAL  CHARACTERISTICS 

CL, EL  -  EPOXY  AND  COPPER  LENGTHS 

CW,EW  -  EPOXY  AND  COPPER  WIDTHS 

UKE.UKC  -  EPOXY  AND  COPPER  THERMAL  CONDUCTIVITIES 

DELE,  TlLl  THROUGH  T4L4  -  THICKNESS  OF  EPOXY  AND  COPPER  LAYERS 

DELX  -  CL/NDEEP 

DELY  -  CW/NWIDE 

NPL  -  NUMBER  OF  NODES  PER  LAYER 
NW'IDE  -  NUMBER  OF  NODES  WIDE 
NDEEP  -  NUMBER  OF  NODES  DEEP 
UL  -  UNITS  OF  LENGTH  (SI  OR  ENGLISH) 

UK  -  UNITS  OF  THERMAL  CONDUCTIVITY  (SI  OR  ENGLISH) 

INITIAL  AND  AMBIENT  TEMPERATURES 

IT  -  INITIAL  BOARD  TEMPERATURE 
LT  -  LEFT  SIDE  AMBIENT  TEMPERATURE 
RT  -  RIGHT  SIDE  AMBIENT  TEMPERATURE 
FT  -  FRONT  AMBIENT  TEMPERATURE 
BT  -  BACK  AMBIENT  TEMPERATURE 
UPRT  -  UPPER  AMBIENT  TEMPERATURE 
LWRT  -  LOWER  AMBIENT  TEMPERATURE 

UT  -  UNITS  OF  TEMPERATURE  (CENTIGRADE  OR  FAHRENHEIT) 

HEAT  INPUT 

THEAT  -  TOTAL  INJECTED  HEAT 

THPN  -  TOTAL  HEAT  PER  NODE 

AHEAT  -  AVERAGE  HEAT  OVER  A  GIVEN  SURFACE 

NHEAT  -  NODE  PER  HEAT  INJECTED  NODE  PER  NODE 

UH,  UAH  -  UNITS  OF  HEAT  (SI  OR  ENGLISH) 

NUM,  NUMA,  CH,  H  -  DUMMY  VARIABLES 

TOTNOD,  NN,  NC,  N  -  VARIABLES  USED  TO  ALLOW  NODAL  HEAT  INPUT 
IH,  JH  -  VECTORS  USED  TO  CORRELATE  NODE  NUMBER  WITH 
MATRIX  POSITION 

HEAT  -  MATRIX  USED  TO  HOLD  HEAT  INPUTS 


COEFFICIENT  DEFINITIONS  Y  IMPLIES  WIDTH,  AND  X  IMPLIES  DEPTH 

Z  IMPLIES  HEIGHT 

C  IMPLIES  COPPER,  E  IMPLIES  EPOXY 
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EYLR  -  PROVIDES  COEFFICIENT  FOR  LEFT  OR  RIGHT  EDGE  NODES 
TO  THE  EXTERNAL  NODE,  EYLR  IMPLIES  EPOXY. 

CYLR  -  SAME  AS  ABOVE,  COPPER  LAYER 

EXFB  -  PROVIDES  COEFFICIENT  FOR  FRONT  AND  BACK  EDGE  NODES 
TO  THE  EXTERNAL  NODE. 

CXFB  -  SAME  AS  ABOVE  FOR  COPPER  LAYER 

EYY  -  INTERNAL  COEFFICIENT  IN  THE  Y  DIRECTION,  EPOXY 

CYY  -  INTERNAL  COEFFICIENT  IN  THE  Y  DIRECTION,  COPPER 

EXX  -  INTERNAL  COEFFICIENT  IN  THE  X  DIRECTION,  EPOXY 

CXX  -  INTERNAL  COEFFICIENT  IN  THE  X  DIRECTION,  COPPER 

SZE  -  COEFFICIENT  FOR  EPOXY  TO  UPPER  EXTERNAL  USE 
EZC  -  COEFFICIENT  FOR  COPPER  TO  EPOXY 
PCC  -  PERCENT  COPPER  COVERAGE 

ALL  ONE  CHARACTER  STRINGS  ARE  SIMPLE  'YES’  OR  'NO'  TYPE 
VARIABLES. 

CALL  INTRODUCTION 


IF  (OVR.  EQ.  'Y' )  THEN 
CALL  INTRO 

ELSEIF  (OVR.EQ. 'N')  THEN 

GOTO  899 

ELSE 

GOTO  88888 
ENDIF 

SUPPLY  A  LIST  OF  ACCEPTABLE  LAYER  LEVELS 
BEGIN  PROGRAM 


899 

800 


C 

801 


CALL  CLS 
VRITE(*,800) 

FORMAT( III,'  ********************************************* 


+**** ' , / , 

+ '  **  **'[/’ 
+'  **  PLEASE  SELECT  UPPER  CASE  LETTERS  **'\l\ 

+'  **  PRIOR  TO  BEGINNING  **',/, 

+ '  ** 


+ '  ************************************************* ' , / , 

+ '  ************************************************* ' ) 


WRITE(*,801) 
FORMATC//, ' 


+' 

PRIOR  TO 

+' 

YOU  HAVE 

+' 

DATA. 

+' 

ENTERING  DATA  INTO  THIS  PROGRAM  ENSURE  THAT  ',/, 
A  DRAWING  OF  YOUR  DESIGN  AND  ALL  PERTINENT  ' ,/, 

'  ,///, 

PRESS  <ENTER>  TO  CONTINUE  ' ,2X,  ) 
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READ(*, 802) ANSWER 
802  FORMAT(Al) 

C 

C 

C 

CALL  CLS 

357  WRITE(*,367) 

367  FORMATC ///////, '  THIS  PROGRAM  IS  CAPABLE  OF  OPERATIONS  IN  El 

+THER  SI  OR  ' 

+'  ENGLISH  UNITS.  AFTER  THE  SELECTION  OF  THE  UNITS,  ALL 

+  '  ENTRIES  MUST  BE  COMPATIBLE.  PLEASE  MAKE  YOUR  SELECTION.  '/ 

+/,'  S  FOR  SI  NOTATION’,/, 

+'  E  FOR  ENGLISH  NOTATION:  ’ ,2X,  ) 

READ(*,368)  ANSN 

368  FORMATC Al) 

C  CHECK  FOR  CORRECT  UNIT  SELECTION 
C 

3333  IFCANSN.EQ. 'S' )  THEN 
WRITE(*,369) 

369  FORMATC////,'  YOU  HAVE  SELECTED  SI  NOTATION.') 

ELSEIF  CANSN.  EQ.  'E')  THEN 

WRITEC*,370) 

370  FORMATC////,'  YOU  HAVE  SELECTED  ENGLISH  NOTATION.') 

ELSE 

CALL  CLS 
GOTO  357 
END  IF 
C 

3711  WRITEC*,3710) 

3710  FORMATC//,'  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  A 

+ND  N  FOR  N  '  ,  ) 

READC*,372)  ANSA 
372  FORMATC Al) 

C 

IFCANSA.  EQ.  'Y')  THEN 
GOTO  378 

ELSEIFCANSA.  EQ.  'N')  THEN 
CALL  CLS 
GOTO  357 
ELSE 

CALL  CLS 
GOTO  3333 
378  ENDIF 
C 
C 
C 
C 

7  CALL  CLS 
SPEVAL='A' 

WRITEC*,301) 

301  FORMATC//////,'  THE  PCB  CAN  HAVE  UP  TO  FOUR  COPPER  LAYERS') 
WRITEC*,302) 

302  FORMATC//,'  HOW  MANY  COPPER  LAYERS  DO  YOU  DESIRE?’) 
WRITEC*,303) 

303  FORMATC/,'  PICK  A  NUMBER  1  TO  4;  ' ,2X,  ) 
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RE AD(*, 304) SELECT 
304  FORMAT(Al) 

C 

C 

CALL  CLS 

IF  ( SELECT.  NE.  ' 1 ' . AND. SELECT. NE.  ' 2 ' . AND. SELECT.  NE.  ’ 3 ' .  AND. 
+SELECT. NE.  '4')  THEN 

GOTO  7 


ELSE 

10  WRITE(*,305)SELECT 

305  FORMATC //////,'  YOU  SELECTED  '.Al,'  COPPER  LAYER(S)  FOR  THE  PCB 

+'.//) 

C 

ENDIF 

C 

WRITE(*,306) 

306  FORMATC  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FO 

+R  NO:  ' ,2X,  ) 

READ(*,307)ANS 

307  FORMATC Al) 


IFCANS.EQ.  'N')  THEN 
GOTO  7 
ELSE 

CONTINUE 

ENDIF 

IFCANS.EQ.  'Y')  THEN 
GOTO  9 
ELSE 

CALL  CLS 
GOTO  10 
ENDIF 


CALL  CLS 


IFC SELECT.  EQ.  'l')  THEN 
GOTO  123 
ELSE 

CONTINUE 

ENDIF 

IFC SELECT.  EQ.  ’2')THEN 
GOTO  124 
ELSE 

CONTINUE 

ENDIF 

IFC SELECT. EQ.  ’3')THEN 
GOTO  125 
ELSE 

CONTINUE 

ENDIF 
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IF( SELECT.  EQ.  '4')THEN 
GOTO  126 
ELSE 

CONTINUE 
ENDIF 
GOTO  7 
C 

c 

11  CALL  CLS 

123  WRITE(*,308) 

308  FORMATC //////, '  YOU  SELECTED  1  COPPER  LAYER,  GIVING  YOU  THE  FOLLO 
+WING  ALTERNATIVES:  ' ) 

WRITE(*,309) 

309  FORMATC//,'  1.  10  BY  36  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES. 

2.  10  BY  35  PROVIDES  A  1:3.5  RATIO  WITH  700  NODES.',/, 

+  '  3.  12  BY  25  PROVIDES  A  1. 2: 2. 5  RATIO  WITH  600  NODES,  ',/, 

+'  4.  15  BY  20  PROVIDES  A  1. 5: 2. 0  RATIO  WITH  600  NODES. ',//, 

+'  PLEASE  SELECT  A  NUMBER  1  THROUGH  4:  ' ,2X,  ) 

READ(*,310)LAY1 

310  FORMATC Al) 

IF  (LAYl.NT.  'I'.AND.  LAYl.NE.  '  2* ,  AND.  LAYl.  NE. 

+  '3'.AND.LAY1.NE.  '4')  THEN 
GOTO  11 
ELSE 

3334  CALL  CLS 

WRITE(*,311)LAY1 

311  FORMATC//////,'  YOU  SELECTED  NUMBER  'Al'  OF  THE  FOLLOWING  4  ALTER 
+NATIVES:  ',///, 

+'  1.  10  BY  36  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES. ',/, 

+'  2.  10  BY  35  PROVIDES  A  1:3.5  RATIO  WITH  700  NODES. ' ,/, 

+'  3.  12  BY  25  PROVIDES  A  1. 2: 2. 5  RATIO  WITH  600  NODES. ',/, 

+  '  4.  15  BY  20  PROVIDES  A  1. 5:2.0  RATIO  WITH  600  NODES.  ',///) 

ENDIF 

13  WRITEC*,312) 

312  FORMATC '  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FO 

+R  NO:  ' ,2X,  ) 

READC*,313)ANSL 

313  FORMATC Al) 

C 

C 

IFCANSL.  EQ. 'N')  THEN 
GOTO  11 
ELSE 

CONTINUE 
ENT)  IF 

IFCANSL. EQ. 'Y' )  THEN 
GOTO  12 
ELSE 

GOTO  3334 

12  ENDIF 
GOTO  1234 

14  CALL  CLS 

124  WRITEC*,314) 

314  FORMATC//////,'  YOU  SELECTED  2  COPPER  LAYERS,  GIVING  YOU  THE  POLL 
+OWING  ALTERNATIVES:  ' ) 
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TOITE(*,315) 

315  FORMATC//,'  1.  10  BY  18  PROVIDES  A  1:1.8  RATIO  WITH  720  NODES. 

+  2.  9  BY  20  PROVIDES  A  1:2.11  RATIO  WITH  720  NODES.',/, 

+'  3.  8  BY  20  PROVIDES  A  2: 5  RATIO  WITH  640  NODES. ' ,/, 

+’  4.  10  BY  15  PROVIDES  A  1: 1. 5  RATIO  WITH  600  NODES. ' ,//, 

+'  PLEASE  SELECT  A  NUMBER  1  THROUGH  4:  ' ,2X,  ) 

READ(*,316)LAY2 

316  FORMATC Al) 

IF  ( LAY2.  NE.  ' 1 ' .  AND.  LAY2.  NE.  ' 2 ' .  AND.  LAY2.  NE. 

+  '3'.  AND.LAY2.NE.  '4')  THEN 
GOTO  14 
ELSE 

3335  CALL  CLS 

WRITE(*,317)LAY2 

317  FORMATC//////,'  YOU  SELECTED  NUMBER  'Al'  OF  THE  FOLLOWING  4  ALTER 
+NATIVES: ',///, 

+  '  1.  10  BY  18  PROVIDES  A  1:1.8  RATIO  WITH  720  NODES.',/, 

+'  2.  9  BY  20  PROVIDES  A  1:2.11  RATIO  WITH  720  NODES. ',/, 

+'  3.  8  BY  20  PROVIDES  A  2:5  RATIO  WITH  640  NODES. ',/, 

+  '  4.  10  BY  15  PROVIDES  A  1:1.5  RATIO  WITH  600  NODES.',///) 

ENDIF 

15  TOITEC*,318) 

318  FORMATC '  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FO 
+R  NO;  ' ,2X,  ) 

READC*,319)ANSL 

319  FORMATC Al) 

C 

C 

IFCANSL.EQ.  'N')  THEN 
GOTO  14 
ELSE 

CONTINUE 

ENDIF 

IFCANSL.EQ.  'Y')  THEN 
GOTO  16 
ELSE 

GOTO  3335 

16  ENDIF 
GOTO  1235 

17  CALL  CLS 
125  WRITEC*,320) 

320  FORMATC//////,'  YOU  SELECTED  3  COPPER  LAYERS,  GIVING  YOU  THE  POLL 
+OWING  ALTERNATIVES:  ' ) 

WRITEC*,321) 

321  FORMATC//,'  1.  5  BY  24  PROVIDES  A  1:4.8  RATIO  WITH  720  NODES. 

+',/,'  2.  10  BY  10  PROVIDES  A  1:1  RATIO  WITH  600  NODES.',/, 

+'  3.  12  BY  10  PROVIDES  A  1.2:1  RATIO  WITH  720  NODES. ',/, 

+'  4.  8  BY  15  PROVIDES  A  1:1.875  RATIO  WITH  720  NODES.',//, 

+'  PLEASE  SELECT  A  NUMBER  1  THROUGH  4:  ' ,2X,  ) 

READC*,322)LAY3 

322  FORMATC Al) 

IF  CLAY3.NE.  '  1 ' .  AND.  LAY3.  NE.  '  2' .  AND.  LAY3.  NE. 

+  '3'.AND.  LAY3.NE.  '4')  THEN 
GOTO  17 
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ELSE 

3336  CALL  CLS 
WRITE(*,323)LAY3 

323  FORMATC //////,'  YOU  SELECTED  NUMBER  'Al'  OF  THE  FOLLOWING  4  ALTER 
+NATIVES: ',///, 

+'  1.  5  BY  24  PROVIDES  A  1:4.8  RATIO  WITH  720  NODES. 

+  '  2.  10  BY  10  PROVIDES  A  1;  1  RATIO  WITH  600  NODES. 

+  '  3.  12  BY  10  PROVIDES  A  1.2:1  RATIO  WITH  720  NODES. 

+  '  4.  8  BY  15  PROVIDES  A  1:1.875  RATIO  WITH  720  NODES.  \///) 

ENDIF 

19  WRITE(*,324) 

324  FORMATC*  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FO 

+R  NO:  ' ,2X,  ) 

READ(*,325)ANSL 

325  FORMATC Al) 

C 

C 

IFCANSL.EQ. 'N* )  THEN 
GOTO  17 
ELSE 

CONTINUE 

ENDIF 

IFCANSL.  EQ.  'Y*  )  THEN 
GOTO  18 
ELSE 

GOTO  3336 
18  ENDIF 

GOTO  1236 

20  CALL  CLS 
126  WRITEC*,326) 

326  FORMATC//////.'  YOU  SELECTED  4  COPPER  LAYERS,  GIVING  YOU  THE  FOLL 
+OWING  ALTERNATIVES:  ' ) 

WRITEC*,327) 

327  FORMATC//,'  1.  5  BY  18  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES. 
+',/,'  2.  9  BY  10  PROVIDES  A  1:1.11  RATIO  WITH  720  NODES.',/, 

+'  3.  8  BY  10  PROVIDES  A  1:1.25  RATIO  WITH  640  NODES. ',/, 

+'  4.  13  BY  6  PROVIDES  A  1: 1  RATIO  WITH  624  NODES. ’ ,//, 

+'  PLEASE  SELECT  A  NUMBER  1  THROUGH  4:  ' ,2X,  ) 

READC*,328)LAY4 

328  FORMATC Al) 

IF  C  LAY4.  NE.  ' 1 ' .  AND. LAY4.  NE. ' 2 ’ . AND.  LAY4. NE. 

+  '3'.AND.  LAY4.NE.  '4')  THEN 
GOTO  20 
ELSE 

3337  CALL  CLS 
WRITEC*,329)LAY4 

329  FORMATC//////,'  YOU  SELECTED  NUMBER  'Al'  OF  THE  FOLLOWING  4  ALTER 
+NATIVES: ',///, 

+'  1.  5  BY  18  PROVIDES  A  1:3.6  RATIO  WITH  720  NODES.'./, 

+  '  2.  9  BY  10  PROVIDES  A  1:1.11  RATIO  WITH  720  NODES.  ',/, 

+'  3.  8  BY  10  PROVIDES  A  1:1.25  RATIO  WITH  640  NODES. ',/, 

+'  4.  13  BY  6  PROVIDES  A  1:1  RATIO  WITH  624  NODES. ',///) 

ENDIF 

21  WRITEC*,330) 

330  FORMATC '  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N  FO 
+R  NO;  ' ,2X,  ) 
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READ(*,331)ANSL 
331  FORMAT(Al) 


IF(ANSL.  EQ.  'N')  THEN 
GOTO  20 
ELSE 

CONTINUE 
ENDIF 

IFCANSL.  EQ.  'Y')  THEN 
GOTO  22 
ELSE 

GOTO  3337 

22  ENDIF 

GOTO  1237 

1234  CONTINUE 

IF(LAY1.  EQ.  '  1'  )THEN 
NWIDE=10 
NDEEP=36 
NPL=360 

ELSEIFCLAYl.  EQ.  '2')THEN 
NV^TDE=10 
NDEEP=35 
NPL=350 

ELSEIFCLAYl.  EQ.  '3')THEN 
NWIDE=12 
NDEEP=25 
NPL=300 

ELSEIFCLAYl.  EQ.  '4')THEN 
NWIDE=15 
NDEEP=20 
NPL=300 

ENDIF 

GOTO  1241 

1235  CONTINUE 

IFCLAY2.EQ.  'l')THEN 

NWIDE=10 

NDEEP=18 

NPL=180 

ELSEIFCLAY2.EQ. '2')THEN 
NWIDE=9 
NDEEP=20 
NPL=180 

ELSEIFC  LAY2. EQ.  ' 3 ' )THEN 
NV'IDE=8 
NDEEP=20 
NPL=160 

ELSEIFC LAY2.  EQ.  '4' )THEN 
NWIDE=10 
NDEEP=15 
NPL=150 

ENDIF 

GOTO  1241 

1236  CONTINUE 

IFCLAY3.  EQ.  ' l' )THEN 


NVIDE=5 

NDEEP=24 

NPL=120 

ELSEIF( LAY3.  EQ.  ' 2 ' )THEN 
NVIDE=10 
NDEEP=10 
NPL=100 

ELSEIF(LAY3.  EQ.  '3')THEN 
NWIDE=12 
NDEEP=10 
NPL=120 

ELSEIF(LAY3.EQ. '4')THEN 
NWIDE=8 
NDEEP=15 
NPL=120 

ENDIF 

GOTO  1241 
1237  CONTINUE 

IF(LAY4.  EQ.  ' 1* )THEN 
NWIDE=5 
NDEEP=18 
NPL=90 

ELSEIF(LAY4.  EQ.  '2')THEN 
NWIDE=9 
NDEEP=10 
NPL=90 

ELSEIF( LAY4.  EQ.  ' 3 ' )THEN 
NWIDE=8 
NDEEP=10 
NPL=80 

ELSEIF(LAY4.EQ.  ’4’ )THEN 
NWIDE=13 
NDEEP=6 
NPL=78 

ENDIF 

GOTO  1241 
C 

C  EPOXY  CHARACTERISTICS 

C 

1241  CALL  CLS 

WRITE(*,8000) 

8000  FORMAT(  //  /  J  ' 

+■*******************  '  j  /  j 

+ ' ****************************epoxy  characteristics**************** 

.f '  •k'h-Mtir4!itirifitirfciritir'feie’hieirititisit1riritieirlrkicic-Hic'/cirirfcitirirHirit-feititieie'kititieirkirk'/rkirkir-/c'k 

+******** ' , / , ) 

C 

C  PROVIDE  CORRECT  UNIT  ABBREVIATIONS 

C 

IFCANSN.  EQ.  'S')  THEN 
WRITE(*,5400) 

5400  FORMATC  all  ENTRIES  ARE  IN  SI  NOTATION.  ',/) 

UL=’cm' 

UK='Watts/ciD/C' 

UT='C' 
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UV=' gm' 

ELSEIFCANSN.EQ. 'E' )  THEN 
WRITE(*,5401) 

5401  FORMATC  all  ENTRIES  ARE  IN  ENGLISH  NOTATION',/) 
UL=' in' 

UK='Btu/hr/F' 

UT='F' 

W='oz' 

ENDIF 


C 

C 

c 

WRITE(*,5402)  UL 

5402  FORMATC/,'  ENTER  EPOXY  LAYER  LENGTH  (' ,A2 ,' ):  ' ,2X,  ) 

READ  *,EL 

C 

WRITE(*,5403)  UL 

5403  FORMATC/,'  ENTER  EPOXY  LAYER  WIDTH  C',A2,');  ' ,2X,  ) 

READ  *,EW 

4337  WRITEC*,6403) 

6403  FORMATC//,'  THE  EPOXY  LAYER  THICKNESS  WILL  DEFAULT  TO  0.0625  in  C 
+0. 15875  cm)',/,'  DO  YOU  WANT  TO  CHANGE  THE  EPOXY  THICKNESS?  CY  OR 
+  N):  ' ,2X,  ) 

READC*,5799)  ANSQ 
5799  FORMATCAl) 

IFCANSQ. EQ. 'Y' )  THEN 
GOTO  5798 

ELSEIFCANSQ.EQ.  'N' . AND.  ANSN.  EQ.  'S')  THEN 
DELE=0.  15875 
GOTO  4338 

ELSEIFCANSQ.EQ.  'N' .  AND.  ANSN.  EQ.  'E')  THEN 
DELE=0.  0625 
GOTO  4338 
ELSE 

GOTO  4337 
ENDIF 

5798  WRITEC*,5404)  UL 

5404  FORMATC/,'  ENTER  EPOXY  LAYER  THICKNESS  C ' ,A2, ' ):  ' ,2X,  ) 

READ  *,DELE 

C 

4338  WRITEC*,5405)  UK 

5405  FORMATC/,'  ENTER  EPOXY  LAYER  THERMAL  CONDUCTIVITY  C'.AIO,'):  ',2 

+X,  ) 

READ  *,UKE 
5459  CALL  CLS 
C 

C  MAKE  CHANGES  OR  CORRECTIONS  TO  EPOXY  ENTRIES 
C 


WRITEC*,5406) 

5406  FORMATC///,' 

+  LAYER:  '  , / , ) 

5407  WRITEC*, 89343) 

89343  FORMATC/,' 
WRITEC*, 89344) 

89344  FORMATC/,' 


YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  EPOXY 
EL,UL 

1. )  LENGTH:  '  ,F9.  4 , 1X,A2 ,  ) 

EW,UL 

2.  )  WIDTH:  ' ,F9.4,1X,A2,  ) 
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WRITE(*, 89345)  DELE.UL 

89345  FORMATC/,'  3.  )  THICKNESS:  ' ,F9. 4, 1X,A2,  ) 

WRITE(*, 89346)  UKE.UK 

89346  FORMATC/,'  4.  )  k:  ' ,F9.  4, lX.All ,/// ,  ) 

C 

VRITE(*,5408) 

5408  FORMATC/,'  DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  YES 

+  AND  N  FOR  NO:  ' ,2X,  ) 

READC*,5409)ANSE 

5409  FORMATC Al) 

C 

C 

IFCANSE.  EQ.  'Y')  THEN 

5441  CALL  CLS 
WRITEC*,5410)  EL,UL 

5410  FORMATC////,'  THE  CURRENT  ENTRY  FOR  LENGTH  IS:  ’ ,F9. 4, 1X,A2) 

WRITEC*,5411) 

5411  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  LENGTH?  CY  OR  N): 

+  '  2X  ) 

READC*,5412)ANSEL 

5412  FORMATC Al) 

PRINT  * 

IFC ANSEL.  EQ.  'Y')  THEN 
WRITEC*,5422)  UL 

5422  FORMATC/,'  ENTER  THE  EPOXY  LENGTH  C A2 ,') :  ',2X,  ) 

READ  *,EL 

ELSEIFC ANSEL.  EQ.  'N')  THEN 

GOTO  5440 

ELSE 

GOTO  5441 
ENDIF 
C 

5440  CALL  CLS 

WRITEC*,5442)  EW,UL 

5442  FORMATC////,'  THE  CURRENT  ENTRY  FOR  WIDTH  IS:  ' ,F9. 4 , 1X,A2) 

WRITEC*,5443) 

5443  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  WIDTH?  CY  OR  N): 

+  '  2X  ) 

READC*,5444)ANSEW 

5444  FORMATC Al) 

PRINT  * 

IFCANSEW.  EQ.  'Y' )  THEN 
WRITEC*,5445)  UL 

5445  FORMATC/,'  ENTER  THE  EPOXY  WIDTH  C ' ,A2, ' ):  ' ,2X,  ) 

READ  *,EW 

ELSEIFCANSEW.  EQ.  'N')  THEN 

GOTO  5446 

ELSE 

GOTO  5440 
ENDIF 
C 

5446  CALL  CLS 
WRITEC*,5447)  DELE,UL 

5447  FORMATC////,'  THE  CURRENT  ENTRY  FOR  THICKNESS  IS:  ' ,F9.  4, 1X,A2) 

WRITEC*,5448) 

5448  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  THICKNESS?  CY  OR  N 


T 
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+):  '  2X,  ) 

READ(*,5449)ANSET 

5449  FORMAT(Al) 

PRINT  * 

IF(ANSET. EQ. 'Y')  THEN 
TOITE(*,5450)  UL 

5450  FORMATC/,'  ENTER  THE  EPOXY  THICKNESS  (' ,A2 :  ’,2X,  ) 

READ  *,DELE 

ELSEIF(ANSET. EQ.  'N')  THEN 

GOTO  5451 

ELSE 

GOTO  5446 
ENDIF 
C 

5451  CALL  CLS 

WRITE(*,5452)  UKE.UK 

5452  FORMATC////,'  THE  CURRENT  ENTRY  FOR  THERMAL  CONDUCTIVITY  IS: 

+F9. 4,1X,A10) 

WRITE(*,5453) 

5453  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  THERMAL  CONDUCTIVIT 

+Y?  CY  OR  N):  ' ,2X,  ) 

READC*,5454)ANSEK 

5454  FORMATC Al) 

PRINT  * 

IFCANSEK.  EQ.  'Y')  THEN 
WRITEC*,5455)  UK 

5455  FORMATC/,'  ENTER  THE  THERMAL  CONDUCTIVITY  C ' ,A10, ' ).  ’,2X,  ) 

READ  *,UKE 

ELSEIFCANSEK.EQ. 'N' )  THEN 

GOTO  5456 

ELSE 

GOTO  5451 
ENDIF 

5456  CALL  CLS 
WRITEC*,5457) 

5457  FORMATC///,'  YOU  HAVE  MADE  THE  FOLLOWING  CORRECTIONS  TO  THE 

+  EPOXY  ENTRIES: ' ,//) 

GOTO  5407 
C 
C 

ELSEIFCANSE. EQ. 'N' )  THEN 

GOTO  5458 

ELSE 

GOTO  5459 
ENDIF 
C 
C 
C 

C  CALL  SUBROUTINE  COPPER 
C 
C 
C 

5458  CALL  COPPERC EW, EL, ANSN, UK, SELECT, TlLl ,T2L1 ,T2L2 ,T3L1 ,T3L2 ,T3L3 ,T4L 
+1,T4L2,T4L3,T4L4,SPEVAL,UKC,CL,CW) 

C 

C 
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nnoono  nonoon  nnnonrjooooonnnoorj  onooo  onoonowonono 


CALL  SUBROUTINE  PCBSl 


CALL  PCBSl 

CALL  SUBROUTINE  PCBS2 


CALL  PC B S 2 ( THE AT , THPN , AHE AT , NHEAT , NPL , NW IDE , NDEEP , HEAT , IH , JH , ANSN , 
+EL,EW) 

CALL  SUBROUTINE  PCBS3 


CALL  PCBS3(EW, EL, CW, CL, UKE, DELE, UKC, SELECT, T1L1,T2L1,T2L2,T3L1,T3L 
+2 , T3L3 , T4L1 , T4L2 , T4L3 , T4L4 , NWIDE , NDEEP , NPL , IH , JH , HEAT , COEF , IBT , UPR 
+T , LWRT , RT , LT , FT , BT , ANSN) 


SUBROUTINE  INTRO 


TITLE:  MODEL  BUILDER 

SUBROUT:  INTRO 

AUTHOR:  LT  STEVE  GLASER 

DATE:  09  JUL  1991 

COMPILER:  MICROSOFT  VERSION  4.  01 

LINKER:  MICROSOFT  VERSION  3.55 


THIS  subroutine:  provides  the  user  with  an  overview  of  the  thermal 

ANALYZER  MODEL  BUILDER. 


COMMON  I BT.UPRT, LWRT, RT,LT, FT, BT,UT, ANSWER 

DEFINE  ONE  VARIABLE  CHARACTER  VALUES 
CHARACTER*!  ANSWER 

PROVIDE  THE  USER  WITH  A  PROGRAM  OVERVIEW 
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7002  FORMAT  (Al) 

CALL  CLS 
UT^ITE(*,7003) 

7003  FORMAT(/////A '  *********  OVERVIEW  ********** 

+'  THIS  PROGRAM  PERFORMS  A  NODAL  ANALYSIS  OF  A  PRINTED 
+'  CIRCUIT  BOARD  CONTAINING  UP  TO  FOUR  COPPER  LAYERS 

+'  (WITH  EPOXY  LAYERS  IN  BETWEEN).  THE  OUTPUT  CONSISTS 
+'  OF  UP  TO  720  COEFFICIENTS  THAT  CONTRIBUTE  TO  THE 
+'  DETERMINATION  OF  THE  TEMPERATUTIE  DISTRIBUTION 
+'  OF  THE  PRINTED  CIRCUIT  BOARD  (PCB)  WHEN  FED  INTO  THE 
+'  THERMAL  ANALYZER. 

+'  THE  FOLLOWING  IS  AN  OUTLINE  OF  THE  MAJOR 
+'  SECTIONS  OF  THIS  PROGRAM  AND  WHAT  ENTRIES  ARE 
+'  REQUIRED  OF  THE  USER. 

+  ' 

+'  PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 
+  ' 

+'  PRESS  <ENTER>  TO  CONTINUE 

RE AD(*, 7 002) ANSWER 
C 

CALL  CLS 
C 

WRITE(*,7004) 

7004  FOR.MAT  (//////,'  *********  OVERV I EW  ********** 

+'  A.  DATA  OUTPUT  FILE: 

+  ' 

+'  THIS  PROGRAM  GENERATES  AN  OUTPUT  DATA  FILE  WHICH  IS 
+'  TO  BE  THE  INPUT  TO  THE  THERMAL  ANALYZER. 

+'  DURING  THE  COURSE  OF  THE  PROGRAM,  THE  USER  WILL  BE 
+'  ASKED  TO  PROVIDE  A  NAME  FOR  THE  OUTPUT  FILE. 

+  ' 

+'  WHEN  PROMPTED  PLEASE  ENTER  THE  NAME  OF  THE  OUTPUT 
+'  DATA  FILE.  THE  DATA  FILE  NAME  SHOULD  BE  NO  LONGER 
+'  THAN  SIX  LETTERS,  AND  MAY  NOT  HAVE  ANY  SPACES. 

+  ' 

+'  PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 
+  ' 

+'  PRESS  <ENTER>  TO  CONTINUE 

READ(*,7002)ANSWER 
C 


t 


/, 

,/ 

J 

/. 

/, 

/, 

/, 

/. 

,/ 

,/ 

J 

! 


7005 


CALL  CLS 
WRITE(*,7005) 

FORMATC //////, '  *********  OVERVIEW  ********** 

+'  B.  STRUCTURE  PHYSICAL  CHARACTERISTICS 
+  ' 

+'  1.  YOU  WILL  BE  ASKED  TO  SELECT  UNIT  TYPE 

+'  (SI  OR  ENGLISH). 

+’  2.  THE  PRINTED  CIRCUIT  BOARD  IS  MAINLY  COMPOSED  OF  ' 

+'  ALTERNATING  COPPER  AND  EPOXY  LAYERS.  THE  PROGRAM' 

+'  IS  DESIGNED  TO  PROVIDE  THE  USER  WITH  A  RANGE  ' 

+'  OF  1  TO  4  COPPER  LAYERS.  IT  IS  ASSUMED  THAT  EACH  ' 

+'  COPPER  LAYER  LIES  BETWEEN  EPOXY  LAYERS. 

+’  DEPENDING  ON  THE  NUMBER  OF  COPPER  LAYERS  SELECTED, 

+'  THE  PROGRAM  WILL  ALLOW  THE  USER  TO  CHOOSE  FROM 


'  ./ 
'  ,/ 
'  ,/ 
'  ,/ 

'  ,/ 


// 


/, 


) 


// 


2X,  ) 


// 
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7006 


7007 


7008 


FOUR  ALTERNATIVE  NODAL  ASPECT  RATIOS. 

PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 


+  ' 

+  ' 

+  ' 

+  ' 

+'  PRESS  <ENTER>  TO  CONTINUE 

READ(*, 7002) ANSWER 
CALL  CLS 
WRITE(*,7006) 

FORMATC //////, '  *********  OVERVIEW  ********** 


,2X,  ) 


+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 


B. 

3. 


STRUCTURE  PHYSICAL  CHARACTERISTICS  (CONTINUED) 


AFTER  SELECTING  THE  DESIRED  NUMBER  OF  COPPER 
LAYERS  AND  NODAL  ASPECT  RATIO,  THE  PROGRAM  WILL 
THEN  ASK  THE  USER  TO  PROVIDE  LAYER  CHARACTERISTICS. 
EPOXY  AND  COPPER  LENGTH,  THICKNESS,  AND  WIDTH, 

AS  WELL  AS  CONSTANTS  CONSTITUTE  THE  QUERIES.  ' , 
THE  PROGRAM  WILL  ALSO  ASK  THE  USER  TO  PROVIDE  ',/ 
THE  PERCENT  COVERAGE  FOR  EACH  COPPER  LAYER.  ' 


’.// 

>  /  j 
/, 


PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 


1. 


PRESS  <ENTER>  TO  CONTINUE 
READ(*, 7002) ANSWER 
CALL  CLS 
WRITE(*,7007) 

FORMATC //////,'  *********  OVERVIEW  ********** 

+'  C.  INTITIAL  AND  AMBIENT  TEMPERATURES 
+  ' 

+  ' 

+  ’ 

+  ' 

+  ' 

+  ' 

+  ' 

+  ’ 

+  ' 

+  ' 

+  ' 


/. 


,2X,  ) 


AFTER  SELECTING  THE  DESIRED  NUMBER  OF  COPPER 
LAYERS  AND  NODAL  ASPECT  RATIO,  THE  PROGRAM  WILL 
THEN  ASK  THE  USER  TO  PROVIDE  TEMPERATURES  FOR  THE 
BOARD. 


PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 


PRESS  <ENTER>  TO  CONTINUE 
READ(*,7002)ANSWER 
CALL  CLS 
WRITE(*,7008) 

FORMATC //////,'  *********  OVERVIEW  ********** 


+  ' 
+  ' 
+  ' 
+  ' 
+  ' 
+  ’ 
+  ' 
+  ' 
+  ’ 
+  ' 
+  ' 
+  ' 
+  ' 
+  ' 


D.  HEAT  INPUT. 

1.  HEAT  INJECTION  OCCURS  ONLY  ON  THE  UPPER  COPPER 
LAYER.  THIS  PROGRAM  SUPPLIES  THE  USER  FOR 
ALTERNATIVE  METHODS  FOR  ENTERING  HEAT. 


A.  )  TOTAL  HEAT  OVER  SURFACE. 

B.  )  AVERAGE  HEAT  PER  UNIT  AREA 

C.  )  INPUT  HEAT  NODE  BY  NODE 


D. )  NO  HEAT  INPUT 
PLEASE  NOTE:  ENTRIES  MUST  BE  IN  UPPER  CASE  LETTERS 
****THIS  CONCLUDES  THE  PROGRAM  OVERVIEW**** 


',// 

/, 

>1, 

’,2X,  ) 


/. 


,//. 
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+'  PRESS  <ENTER>  TO  CONTINUE  ',2X,  ) 

READ(*,7002)ANSWER 
CALL  CLS 
END 

SUBROUTINE  PCBSl 

Q  A***********!/,-******************************************************** 

C 

C  TITLE:  MODEL  BUILDER 

C  SUBROUT:  PCBSl 

C  AUTHOR:  LT  STEVE  GLASER 

C  DATE:  09  MAY  1991 

C  COMPILER:  MICROSOFT  VERSION  4. 01 

C  LINKER:  MICROSOFT  VERSION  3.55 

C 

C  THIS  SUBROUTINE  HANDLES  THE  TEMPERATURE  INPUTS  TO  THE  PRINTED  CIRCUIT 
C  BOARD. 

C 

C 

C 

C 

COMMON  IBT,UPRT,LWRT,RT,LT,FT,BT,UT,ANSWER 
C 

C  DEFINE  REAL  VARIABLES 

C 

C 

REAL  IBT,UPRT,LVRT,RT,LT,FT,BT 
C 

C  DEFINE  ONE  CHARACTER  VARIABLES 

C 

CHARACTER*!  ANSTEM , ANSTEB , ANSTEU , ANSTEL , ANSTER , ANSLFT , ANSFRT , ANSBC 
+K,UT 
C 

31180  CALL  CLS 

WRITE(*, 60000) 

60000  FORMAT( III,' ****************************************************** 

+*********yf*  ■>******* '  ^ 

+ '  *******************^0 TEMPERATURE  INPUTS******************** 
+******** ' , / , 

4- ' ***************************************************************** 
4.*****Vf** ',//,) 

c 

C 

VTRITEC*, 60001)  UT 

60001  FORMATC/,'  ENTER  THE  INITIAL  BOARD  TEMPERATURE  (',A1,’):  ',2X, 

+) 

READ  *,IBT 
C 
C 

WRITE(*, 60002)  UT 

60002  FORMATC/,'  ENTER  THE  UPPER  SURFACE  AMBIENT  TEMPERATURE  (',A1,') 

+:  ’.2X,  ) 

READ  *,UPRT 
C 
C 
C 
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WRITEC*, 60003)  UT 

60003  FORMATC/,'  ENTER  THE  LOWER  SURFACE  AMBIENT  TEMPERATURE  (',A1,') 
+:  '.2X,  ) 

READ  *,LWRT 
C 
C 

WRITEC*, 60004)  UT 

60004  FORMATC/,'  ENTER  THE  RIGHT  SURFACE  AMBIENT  TEMPERATURE  (',A1,') 
+:  ',2X,  ) 

READ  *,RT 
C 
C 
C 

WRITEC*, 60005)  UT 

60005  FORMATC/,'  ENTER  THE  LEFT  SURFACE  AMBIENT  TEMPERATURE  (',A1,'): 
+  '.2X,  ) 

READ  *,LT 
C 
C 


WRITEC*, 60006)  UT 

60006  FORMATC/,'  ENTER  THE  FRONT  SURFACE  AMBIENT  TEMPERATURE  (',A1, 
+:  ',2X,  ) 

READ  *,FT 
C 
C 
C 


') 


WRITEC*, 60007)  UT 

60007  FORMATC/,'  ENTER  THE  REAR  SURFACE  AMBIENT  TEMPERATURE  C',A1,'): 
+  ’,2X,  ) 

READ  *,BT 


C 

C 

C  REVIEW  THE  TEMPERATURE  ENTRIES 
C 

60052  CALL  CLS 

WRITEC*, 60008)  IBT,UT 

60008  FORMATC////,'  YOU  HAVE  MADE  THE  FOLLOWING  AMBIENT  TEMPERATURE  E 
+NTRIES:  ',////,'  1.  INITIAL  BOARD  TEMPERATURE:  ' , 1X,F9.  3, 1X,A1) 

WRITEC *,60009)  UPRT,UT 

2.  UPPER  AMBIENT  TEMPERATURE: 


60009  FORMATC' 

WRITEC *,600 10)  LWRT,UT 

60010  FORMATC '  3.  LOWER  AMBIENT  TEMPERATURE 

WRITEC*, 60011)  RT,UT 

60011  FORMATC'  4.  RIGHT  AMBIENT  TEMPERATURE 
WRITEC *,600 12)  LT,UT 

5.  LEFT  AMBIENT  TEMPERATURE: 
WRITEC*, 60013)  FT,UT 

60013  FORMATC '  6.  FRONT  AMBIENT  TEMPERATURE: 

WRITEC*, 60014)  BT,UT 

60014  FORMATC'  7.  REAR  AMBIENT  TEMPERATURE: 

C 

C 

C 


60012  FORMATC' 


,1X,F9. 3,1X,A1) 
,1X,F9.  3,1X,A1) 
,1X,F9.  3,1X,A1) 
,1X,F9.  3,1X,A1) 
,1X,F9.  3,1X,A1) 
,1X,F9.3,1X,A1,//) 


WRITEC *,600 15) 

60015  FORMATC/,'  DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  YES  A 
+ND  N  FOR  NO:  ' ,2X,  ) 
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READ(*,60016)ANSTEM 

60016  FORMAT(Al) 

C 

C 

C  MAKE  CORRECTIONS  BOARD  TEMPERATURES 

C 

C 

IF(ANSTEM.  EQ.  'Y' )  THEN 

60022  CALL  CLS 

WRITEC*, 60017)  IBT.UT 

60017  FORMATC///,'  THE  INITIAL  BOARD  TEMPERATURE  IS;  ' ,F9. 3 , lX.Al) 
WRITE(*, 60018) 

60018  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  (Y  OR  N):  ',2 
+X,  ) 

READC*, 60019)  ANSTEB 

60019  FORMAT(Al) 

IF( ANSTEB.  EQ.  'Y' )  THEN 
WRITEC*, 60020)  UT 

60020  FORMATC/,'  ENTER  THE  NEW  VALUE  C ' ,A1 ,' ):  ' ,2X,  ) 

READ  *,IBT 

ELSE IFC ANSTEB.  EQ.  'N')  THEN 
GOTO  60021 
ELSE 

GOTO  60022 
END  IF 
C 

60021  CALL  CLS 

WRITEC*, 70022)  UPRT.UT 

70022  FORMATC///,'  THE  UPPER  AMBIENT  TEMPERATURE  IS:  ' ,F9. 3 , 1X,A1) 
WRITEC*, 60023) 

60023  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ’,2 
+X,  ) 

READC*, 60024)  ANSTEU 

60024  FORMATC Al) 

IFC ANSTEU.  EQ.  'Y')  THEN 
WRITEC*, 60025)  UT 

60025  FORMATC/,'  ENTER  THE  NEW  VALUE  C ' ,A1 ,' ):  ' ,2X,  ) 

READ  *,UPRT 

ELSEIFC ANSTEU.  EQ.  'N')  THEN 
GOTO  60026 
ELSE 

GOTO  60021 
END  IF 
C 

c 

60026  CALL  CLS 

WRITEC*, 60027)  LWRT,UT 

60027  FORMATC///,'  THE  LOWER  AMBIENT  TEMPERATURE  IS:  ' ,F9.  3 , 1X,A1) 
WRITEC*, 60028) 

60028  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ',2 
+X,  ) 

READC*, 60029)  ANSTEL 

60029  FORMATC Al) 

IFC ANSTEL.  EQ.  'Y')  THEN 
WRITEC*, 60030)  UT 

60030  FORMATC/,'  ENTER  THE  NEW  VALUE  C ' ,A1 ,' );  ' ,2X,  ) 
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READ  *,LWRT 

ELSEIF(ANSTEL.  EQ. 'N')  THEN 
GOTO  60031 
ELSE 

GOTO  60026 
ENDIF 
C 
C 

c 

60031  CALL  CLS 

WRITEC*, 60032)  RT.UT 

60032  FORMATC///,'  THE  RIGHT  AMBIENT  TEMPERATURE  IS:  ' ,F9. 3 , 1X,A1) 
WRITEC*, 60033) 

60033  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ',2 
+X,  ) 

READC*, 60034)  ANSTER 

60034  FORMATC Al) 

IFC ANSTER.  EQ.  'Y')  THEN 
WRITEC*, 60035)  UT 

60035  FORMATC/,'  ENTER  THE  NEW  VALUE  C ' ,A1 ,' ):  ' ,2X,  ) 

READ  *,RT 

ELSEIFC ANSTER.  EQ. 'N')  THEN 
GOTO  60036 
ELSE 

GOTO  60031 
ENDIF 
C 
C 

c 

c 

60036  CALL  CLS 

WRITEC*, 60037)  LT.UT 

60037  FORMATC///,'  THF.  LEFT  AMBIENT  TEMPERATURE  IS:  ' ,F9. 3, lX.Al) 
WRITEC*, 60038) 

60038  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ',2 
+X,  ) 

READC*, 60039)  ANSLFT 

60039  FORMATC Al) 

IFC ANSLFT.  EQ.  'Y' )  THEN 
WRITEC *,60040)  UT 

60040  FORMATC/,'  ENTER  THE  NEW  VALUE  C ’ ,A1 ,' ):  ' ,2X,  ) 

READ  *,LT 

ELSEIFC ANSLFT.  EQ.  'N')  THEN 
GOTO  60041 
ELSE 

GOTO  60036 
ENDIF 
C 
C 
C 
C 

60041  CALL  CLS 

WRITEC *,60042)  FT,UT 

60042  FORMATC///,'  THE  FRONT  AMBIENT  TEMPERATURE  IS:  ' ,F9. 3, 1X,A1) 
WRITEC *,60043) 

60043  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ',2 
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+x,  ) 

READ(*, 60044)  ANSFRT 

60044  FORMAT(Al) 

IF( ANSFRT.  EQ.  'Y' )  THEN 
WRITEC*, 60045)  UT 

60045  FORMATC/,'  ENTER  THE  NEW  VALUE  (' ,A1 ,') :  ' ,2X,  ) 
READ  *,FT 

ELSEIFC ANSFRT. EQ. 'N')  THEN 
GOTO  60046 
ELSE 

GOTO  60041 
END  IF 


C 

C 

C 

C 

60046  CALL  CLS 

WRITEC*, 60047)  BT.UT 

60047  FORMATC///,'  THE  REAR  AMBIENT  TEMPERATURE  IS;  ' ,F9.  3, 1X,A1) 
WRITEC*, 60048) 

60048  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THIS  VALUE?  CY  OR  N):  ',2 
+X,  ) 

READC*, 60049)  ANSBCK 

60049  FORMATC Al) 

IFC ANSBCK. EQ. 'Y')  THEN 
WRITEC*, 60050)  UT 

60050  FORMATC/,'  ENTER  THE  NEW  VALUE  C ' ,A1 ,' ):  ',2X,  ) 

READ  *,BT 

ELSEIFC ANSBCK.  EQ.  'N' )  THEN 
CONTINUE 
ELSE 

GOTO  60046 
END  IF 

ELSEIFC ANSTEM.  EQ. 'N' )  THEN 
GOTO  60051 
ELSE 

GOTO  60052 
END  IF 

60051  END 

SUBROUTINE  PCBS2CTHEAT,THPN,AHEAT,NT{EAT,NPL,NWIDE,NDEEP,HEAT,IH,JH 
+,ANSN,EL,EW) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


TITLE: 

SUBROUT: 

AUTHOR: 

DATE: 

COMPILER: 

LINKER: 


MODEL  BUILDER 

PCBS2  --  HEAT  INPUT  SUBROUTINT; 
LT  STEVE  GLASER 
09  MAY  1991 

MICROSOFT  VERSION  4.  01 
MICROSOFT  VERSION  3. 55 


THIS  SUBROUTINE  HANDLES  THE  HEAT  INPUTS  TO  THE  PRINTED  CIRCUIT  BOARD 


DEFINE  REAL  VARIABLES 


REAL  THEAT,THPN,AHEAT,NHEAT,SL,SW 
C 
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c 

C  DEFINE  INTEGERS 

C 

INTEGER  NWIDE ,NDEEP,NPL,NUM,NUMA,CH,H,TOTNOD,NC,NN, I , J, IM 
C 

c 

c  DEFINE  ONE  CHARACTER  VARIABLES 

C 

CHARACTER* 1  SE  LH , ANSH , ANSN , ATH , ANSHA , AHN , DUMMY 
C 

C  DEFINE  CHARACTER  VARIABLES 

C 

CHARACTER  UH*13,UAH*6 
C 

C  DEFINT  REAL  MATRICES 

C 

REAL  HEAT(100,100),COEFF(740,9) 

C 

C  DEFINE  INTEGER  VARIABLES 

C 

INTEGER  IH( 1000), JH( 1000), UMMY( 1000) 

C 

C  PROVIDE  A  CORRELATION  BETWEEN  NODE  NUMBERS  AND  MATRIX  LOCATION 

C 

C 

51  NUM=1 

DO  60  I=1,NPL/NWIDE 
DO  61  J=1,N\^’IDE 
JH(NUM)=J 
NUM=NUM+1 

61  CONTINUE 

60  CONTINUE 

C 

NUMA=1 

CH=0 

H=1 

DO  62  1=1, NPL 
IH(NUMA)=H 
CH=CH+1 

IF(CH.EQ. NWIDE)  THEN 
H=H+1 
CH=0 
ELSE 

CONTINUE 
END  IF 

NUMA=NUMA+1 

62  CONTINUE 
C 

C 

C 

1241  CALL  CLS 

WRITE(*,8000) 

8000  FORMAT (  /  //  ,  ' 

+***:fr*y<************* '  ,  /  , 

+ '  ***yr************yr**************}{EAT  INPUTS*********************** 
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+ ' **★★*****★******★★★**★*******■*★*****★★*******★****★★************* 
+******** ' , / / , 

+'  heat’ INPUT  TO  THE  PCB  OCCURS  ONLY  ON  THE  UPPER  COPPER',/, 

+'  LAYER.  HEAT  INPUT  IS  ACCOMPLISHED  BY  ONE  OF  THE  ’,/, 

+'  FOLLOWING  METHODS:  ',//, 

+'  1.  ENTER  AS  TOTAL  HEAT  APPLIED  TO  THE  PCB’ ,/, 

+'  2.  ENTER  AS  AVERAGE  HEAT  PER  UNIT  AREA  ’,/, 

+'  3.  ENTER  HEAT  NODE  BY  NODE  ',/, 

+'  4.  NO  HEAT  INPUT' ,//, 

+'  PLEASE  SELECT  NUMBER  1  THROUGH  4:  ',2X,  ) 

C 

READ(*,5799)  SELH 
5799  FORMATCAl) 

IF(  SELH.  EQ.  '  1 ' .  OR.  SELH.  EQ.  '  2 ' .  OR.  SELH.  EQ.  ’  3 ' .  OR.  SELH.  EQ.  '  4 '  )  THEN 
444  WRITE(*,291)  SELH 

291  FORMATC//,'  YOU  HAVE  SELECTED  NUMER  ',A1,'  OF  FOUR  ALTERNATIVES 
+.  ',/,'  IS  THIS  THE  DESIRED  SELECTION?  (Y  OR  N):  ' ,2X,  ) 

ELSE 

GOTO  1241 
END  IF 

READ(*,5409)ANSH 
5409  FORMATCAl) 

C 

C 

IFCANSH. EQ. 'Y')  THEN 
5441  GOTO  63 

ELSEIFCANSH. EQ. 'N' )  THEN 
GOTO  51 
ELSE 

CALL  CLS 
WRITE(*,8000) 

GOTO  444 
63  ENDIF 

C 
C 

C  DETERMINE  UNIT  FOR  HEAT  INPUT 
C 

IF  (CANSN.EQ.  'E').AND.(SELH.EQ.  '2'))  THEN 
UH  =  'Btu/(hr*in  2)' 

UAH  =  'Btu/hr' 

ELSEIF  (CANSN.  EQ.  'E').AND.  (SELH.  EQ. 'I'.OR.  SELH.  EQ.  '3'))  THEN 
UH  =  'Btu/hr' 

ELSEIF  (CANSN.  EQ.  ’S').  AND. (SELH.  EQ. '2'))  THEN 
UH  =  'Watts/(cm  2)' 

UAH  =  'Watts' 

ELSEIF  ((ANSN.EQ.  'S').AND.(SELH.EQ. 'I'.OR.SELH.EQ. '3'))  THEN 
UH  =  'Watts' 

ENDIF 

C 

C 

C  ALLOW  FOR  RE -SELECTION  OF  HEAT  INPUT  METHOD  OR  CONTINUE  WITH 
C  INITIAL  SELECTION. 

C 

C  CHOICE  in 

C 

C 
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CALL  CLS 

IF  (SELH.  EQ.  'I'.AND.  ANSN.  EQ.  'E')  THEN 

989  WRITE(*,988)  UH 

988  FORMATC///,'  YOU  HAVE  SELECTED  TO  INPUT  HEAT  AS  A  TOTAL  HEAT 

+  APPLIED  TO  THE  SLWACE. 

+'  ENTER  TOTAL  HEAT  APPLIED  TO  THE  SURFACE  (',A6,’):  ' ,2X,  ) 

READ  *,THEAT 

990  WRITE(*,991) 

991  FORMATC/,'  IS  THIS  THE  CORRECT  ENTRY?  (Y  OR  N):  ' ,2X,  ) 

READ(*,992)ATH 

992  FORMATC Al) 

C 

C  MAKE  ENTRY  AND  ALLOW  FOR  CORRECTIONS 
C 

IF  CATH.  EQ.  'Y')  THEN 
THPN=THEAT/NPL 
WRITEC*,993)  THPN.UH 

993  FORMATC/,'  TOTAL  HEAT  PER  NODE  IS:  ' ,F9. 4 , 1X,A13) 

WRITEC* , 9323) 

9323  FORMATC//,'  <PRESS  ENTER  TO  CONTINUE>') 

READC*,9324)DUMMY 

9324  FORMATC Al) 

C 

C 

C  FILL  HEAT  MATRIX  WITH  DESIRED  VALUES 

C 

C 

DO  994  I=1,NDEEP 
DO  995  J=1,NWIDE 
HEATCI,J)=THPN 

995  CONTINUE 

994  CONTINUE 
C 

C 

ELSEIFCATH.  EQ.  'N' )  THEN 
CALL  CLS 
GOTO  989 
ELSE 

CALL  CLS 

WRITEC*,996)  THEAT,UH 

996  FORMATC////,'  TOTAL  HEAT  APPLIED  TO  THE  SURFACE  IS:  ' ,F9 

+.4,1X,A13,  ) 

GOTO  990 

END  IF 
C 
C 

C  CHOICE  #2 

C 

C 

ELSEIFCSELH.  EQ.  '2')  THEN 
998  WRITEC*,997)  UH 

997  FORMATC////,'  YOU  HAVE  SELECTED  TO  ENTER  THE  AVERAGE  HEA 

+T  OVER  THE',/,'  UPPER  PCB  SURFACE.',//,'  ENTER  THE  DESIR 

+ED  HEAT  INPUT:  C',A13,'):  ' ,2X,  ) 

READ  *,AHEAT 
*C 
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C  MAKE  ENTRY  AND  ALLOW  FOR  CORRECTION 
C 

1000  WRITE(*,999) 

999  FORMAT(/,'  IS  THIS  THE  CORRECT  ENTRY?  (Y  OR  N):  ' ,2X,  ) 

READ(*,1001)ANSHA 

1001  FORMAT(Al) 

IF  (ANSHA. EQ. 'Y' )  THEN 
THPN=AHEAT*EL*EW/NPL 
WRITE(*,1002)  THPN,UAH 

1002  FORMATC///,'  TOTAL  HEAT  PER  NODE  IS:  ' ,F9. 4, 1X,A6) 

WRITE(*,4323) 

4323  FORMATC//,'  <PRESS  ENTER  TO  CONTINUE> ' ) 

READ(*,4324)DUMMY 

4324  FORMATC Al) 

C 

C  FILL  HEAT  MATRIX  WITH  DESIRED  VALUES 

C 

C 

DO  1003  I=1,NDEEP 
DO  1004  J=1,NWIDE 
HEATCI,J)=THPN 
1004  CONTINUE 

1003  CONTINUE 
C 

C 

ELSEIFC ANSHA.  EQ.  'N')  THEN 
CALL  CLS 
GOTO  998 
ELSE 

CALL  CLS 

WRITEC*,1006)  AHEAT.UH 

1006  FORMATC////,'  AVERAGE  HEAT  OVER  PCB  SURFACE  IS:  ’,F9.4,1 

+X,A13,  ) 

GOTO  1000 

ENDIF 

C 

C  CHOICE  #3 

C 

C 

ELSEIFC SELH.  EQ.  '3' )  THEN 

1200  WRITEC*,1201) 

1201  FORMATC////,’  YOU  HAVE  SELECTED  TO  ENTER  THE  HEAT  NODALLY' 

+,//,'  ENTER  THE  TOTAL  NUMBER  OF  NODES  DESIGNATED  FOR  HEAT  INP 

+UT:  ',2X,  ) 

READ  *,TOTNOD 
C 
C 

C  THIS  IS  DONT  NODE  BY  NODE.  GET  NUMBER  OF  ENTRIES  AND  THEN  LOOP  UNTIL 
C  ALL  ENTRIES  HAVE  BEEN  MADE. 

C 

C 

C  TELL  USER  MAXIMUM  NUMBER  OF  ENTRIES  POSSIBLE 
C 

IF  CTOTNOD.GT.NPL)  THEN 
WRITEC*,1202)  NPL 

1202  FORMATC///,'  THE  MAXIMUM  ENTRY  IS:  ’,14) 
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WRITE(*,2239) 

2239  FORMAK///,'  PLEASE  PRESS  <ENTER>  TO  CONTINUE  ') 

READ(*,2240)  DUMMY 

2240  FORMAT(Al) 

CALL  CLS 
GOTO  1200 
END  IF 

C 

C  MAKE  ENTRIES 

C 

C 

DO  1203  I=l,TOTNOD 
NC=I 

CALL  CLS 

WRITE(*,1205)  NC.TOTNOD 

FORMATC///,'  THIS  IS  NUMBER  ',13,'  OF  ',13,'  ENTRIES') 

WRITE(*,5345) 

FORMATC//,'  ENTER  THE  NODE  NUMBER  FOR  HEAT  INPUT:  ',2 

READ  (*,8032)  NN 
FORMATC 14) 

UMMY(I)=NN 

IF  (NN.  EQ.  O.OR.  NN.  GT.NPL)  THEN 
GOTO  1204 
ENDIF 

WRITE(*,1206)  UH 

FORMATC/,'  ENTER  THE  HEAT  INPUT  (',A6,'):  ',2X,  ) 

RFAD  *  fJHFAT 

HEAT(IH(NN),JH(NN))  =  NHEAT 
1203  CONTINUE 
C 

C  PROVIDE  OPPORTUNITY  TO  MAKE  CORRECTIONS 

C 

C 

1301  CALL  CLS 

WRITE(*,1302)  TOTNOD 

1.302  FORMATC///,'  YOU  HAVE  MADE  ',13,'  NODAL  ENTRIES:  ',/) 

DO  11203  IM=1, TOTNOD 

WRITEC*, 11204)  UMMY(IM),HEAT(IH(UMMY(IM)),JH(UMMYCIM))),UH 
11204  FORMATC IX,'  NODE  NUMBER  ' ,14, ' : ' , 1X,F9. 4, 1X,A6) 

11203  CONTINUE 

WRITEC*, 1303) 

1303  FORMATC/,'  DO  YOU  WISH  TO  MAIUi  ANY  MORE  ENTRIES  OR  CORRECTION 

+S?  (Y  OR  N):  ' ,2X,  ) 

READ(*,1304)  AHN 

1304  FORMATC Al) 

IFCAHN.EQ.  'Y')  THEN 

CALL  CLS 
GOTO  1200 

ELSEIFCAHN.  EQ.  'N')  THEN 
GOTO  1305 
ELSE 

GOTO  1301 

1305  ENDIF 

ELSEIFCSELH. EQ. '4' )  THEN 


1204 

1205 
5345 

+X,  ) 

8032 

1206 
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ooocjno  nnnoonnonoonorjonono 


GOTO  1306 
1306  ENDIF 
END 

SUBROUTINE  PCBS3(XEW, XEL,XCW,XCL,XUKE,XDELE,XUKC, SELECT, X1L1,X2L1, 
+X2L2,X3L1,X3L2,X3L3,X4L1,X4L2,X4L3,X4L4.XWIDE,XDEEP,XPL,XIH,XJH,X 
+E  AT ,  XOEF ,  X I BT ,  XPRT ,  XVvTlT ,  XRT ,  XLT ,  XFT ,  XBT ,  XN  ITS  ) 


TITLE: 

SUBROUTINE: 

AUTHOR: 

DATE: 

COMPILER: 

LINKER: 


MODEL  BUILDER 
PCBS3 

LT  STEVE  GLASER 
09  JUL  1991 

MICROSOFT  VERSION  4.  01 
MICROSOFT  VERSION  3. 55 


THIS  SUBROUTINE  IS  THE  "HEART"  OF  THE  MODEL  BUILDER.  THE 
NODE  COEFFICIENTS  ARE  CALCULATED  IN  THIS  SUBROUTINE  AND 
PLACED  IN  MATRICES.  PCBS3  IN  TURN  ALSO  CALLS  OTHER 
SUBROUTINES:  SI,  S2,  S3,  AND  S4.  THESE  SUBROUTINES  CALLED 
BY  PCBS3  GENERATE  THE  OUTPUT  DATA  FILE  FOR  THE  THERMAL 
ANALYZER. 


DEFINE  REAL  VARIABLES 


REAL  XEL,XEW,XDELE,XUKE,XCL,XCW,X1L1,X2L1,X2L2,X3L1,X3L2,X3L3,X4L1 
+ , X4L2 , X4L3 , X4L4 , XUKC , EZB 

REAL  EYLR,EFXB,EYY,EZC11,EZC21,EZC22,EZC31,EZC32,EZC33,EZC41,EZC42 
+,EZC43,EZC44,EXX,DELX,DELY 

REAL  CYLR 1 1 , CYLR2 1 , CYLR22 , CYLR3 1 , CYLR32 , CYLR33 , CYLR4 1 , CYLR42 , CYLR4 
+3,CYLR44 

REAL  CFXB 1 1 , CFXB2 1 , CFXB22 , CFXB31 ,CFXB32 , CFXB33 , CFXB41 , CFXB42 , CFXB4 
+3,CFXB44 

REAL  CYY 1 1 , CYY2 1 , CYY22 , CYY3 1 , CYY32 , CYY33 , CYY4 1 , CYY42 , CYY43 , CYY44 
REAL  CXX 1 1 , CXX2 1 , CXX22 , CXX3 1 , CXX32 , CXX3 3 , CXX4 1 , CXX42 , CXX43 , CXX44 
REAL  CZEll ,CZE21 ,CZE22 ,CZE31 ,CZE32 ,CZE33,CZE41 ,CZE42 ,CZE43 ,CZE44 
REAL  PC  1 1 , PC2 1 , PC22 , PC3 1 , PC 32 , PC33 , PC4 1 , PC42 , PC43 , PC44 
REAL  XEAT(100,100),XOEF(1000,12) 


DEFINE  CHARACTER  VARIABLES 


CHARACTER*!  SELECT, ANS ,XNITS 
DEFINE  INTEGER  VARIABLES 
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c 

c 

INTEGER  XVv’IDE  ,XDEEP  ,XPL,NUMC ,  I  ,N,  IB  ,USEL 

INTEGER  XIH( 1000), XJH( 1000) 

C 

C*****vr*****Vr*****-,v*****D£7ERNIj{v;£  COEFFICIENTS************************** 

C********************************* ************************************** 

C 

C 

C  DETERMINE  INCREMENTAL  MEASUREMENTS  IN  THE  X  AND  Y  DIRECTIONS 

C 

C 

C  XNITS  IS  A  CODE  USED  IN  THE  OUTPUT  DATA  FILE.  1  MEANS  ENGLISH  UNITS, 
C  AND  2  MEANS  SI  UNITS 
C 
C 

IF  (XNITS. EQ.  'E' )  THEN 
USEL=1 
ELSE 

USEL=2 

ENDIF 

C 

C  DETERMINE  INCREMENTS  IN  THE  ARRAY 
C 

DELX  =  XEL/XDEEP 
DELY  =  XEW/XVTDE 
nil  CONTINUE 
C 
C 

C  COPPER  COVERAGE,  ONT  COPPER  LAYER  CASE 

C 

C 

IF( SELECT. EQ. ' l' )  THEN 
931  CALL  CLS 
C 

C  ENTER  PERCENT  COPPER  COVERAGE 

C 

C 

WRITE(*,39A1) 

3941  FORMAT(////, '  ENTER  THE  PERCENT  COPPER  COVERAGE  (i. e.  67)  FOR 

+  LAYER  1:  ' ,2X,  ) 

READ  *,PC11 

IF  (PCll. GT. 100)  THEN 

GOTO  931 

ENDIF 

PC11=PC11/100 

C 

C 

C  COPPER  COVERAGE,  TW’O  COPPER  LAYERS  CASE 
C 
C 
C 

c 
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ELSEIFCSELECT. EQ. '2' )  THEN 

932  CALL  CLS 
VRITE(*,3942) 

3942  FORMATC////, '  ENTER  THE  PERCENT  COPPER  COVERAGE  (i. e.  67)  FOR 

+  LA^-ER  1:  '  ,2X,  ) 

READ  *,PC21 
IF  (PC21. GT. 100)  THEN 
GOTO  932 
END  IF 

PC21=PC21/100 

933  WRITE(*,3943) 

3943  FORMATC/,'  ENTER  THE  PERCENT  COPPER  COVERAGE  (i. e.  67)  FOR  LA 

+YER  2:  ' ,2X,  ) 

READ  *,PC22 
IF  (PC22. GT. 100)  THEN 
GOTO  933 
END  IF 

PC22=PC22/100 

C 

C 

C  COPPER  COVERAGE,  THRE^  -  "'PER  LAYERS  CASE 

C 

C 

ELSEIFCSELECT. EQ. ' 3' )  THEN 
942  CALL  CLS 

WRITEC*,4942) 

4942  FORMATC////,'  ENTER  THE  PERCENT  COPPER  COVERAGE  Ci- e.  67)  FOR 

+  LAYER  1:  ' ,2X,  ) 

READ  *,PC31 
IF  CPC31. GT. 100)  THEN 
GOTO  942 
END  IF 

PC31=PC31/100 
433  WRITEC*,4943) 

4943  FORMATC/,'  ENTER  THE  PERCENT  COPPER  COVERAGE  Ci. e.  67)  FOR  LA 

+YER  2:  ' ,2X,  ) 

READ  *,PC32 
IF  CPC32. GT. 100)  THEN 
GOTO  433 
END  IF 

PC32=PC32/100 
633  WRITEC*,8943) 

8943  FORMATC/,'  ENTER  THE  PERCENT  COPPER  COVERAGE  C i. e.  67)  FOR  LA 

+YER  3:  ' ,2X,  ) 

READ  *,PC33 
IF  CPC33. GT. 100)  THEN 
GOTO  633 
END  IF 

PC33=PC33/100 

C 

C 

C  COPPER  COVERAGE,  FOUR  COPPER  LAYERS  CASE 

C 

C 

ELSEIFCSELECT. EQ. '4' )  THEN 
242  CALL  CLS 
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V»'RITE(*,2942) 

2942  FORMAK////, '  ENTER  THE  PERCENT  COPPER  COVERAGE  ( i.  e.  67)  FOR 

+  LAYER  1:  ' ,2X,  ) 

READ  *,PC41 
IF  (PC41.GT.  100)  THEN 
GOTO  242 
ENDIF 

PC41=PC41/100 
233  WRITE(*,2943) 

2943  FORMATC/,'  ENTER  THE  PERCENT  COPPER  COVERAGE  (i. e.  67)  FOR  LA 

+YER  2:  ' ,2X,  ) 

READ  *,PC42 
IF  (PC42.GT.  100)  THEN 
GOTO  233 
ENDIF 

PC42=PC42/100 
133  WRITE(*,1943) 

1943  FORMATC/,'  ENTER  THE  PERCENT  COPPER  COVERAGE  (i. e.  67)  FOR  LA 

+YER  3:  ',2X,  ) 

READ  *,PC43 

IF  (PC43.  GT. 100)  THEN 

GOTO  133 

ENDIF 

PC43=PC43/100 
33  WRITEC*  43) 

43  FORMATC/,’  ENTER  THE  PERCENT  COPPER  COVERAGE  Ci.e.  67)  FOR  LA 

+YER  4:  ' ,2X,  ) 

READ  *,PC44 
IF  CPC44.  GT. 100)  THEN 
GOTO  33 
ENDIF 

PC44=PC44/100 
ENDIF 
CALL  CLS 
C 
C 
C 

C  VERIFY  THAT  THE  SELECTED  ENTRIES  ARE  CORRECT 
C 
C 
C 

IFC SELECT.  EQ.  ' l' )  THEN 
305  WRITEC*, 300)  PC11*100 

300  FORMATC/////'  YOU  HAVE  SELECTED  ',F6.2,'%  COVERAGE.',//,' 

+  IS  THIS  YOUR  DESIRED  ENTRY?  CY  OR  N):  ',1X,  ) 

READ  C*, 30 DANS 

301  FORMATC Al) 

IFCANS. EQ.  'Y' )  THEN 
GOTO  333 

ELSEIFCANS.EQ.  'N' )  THEN 

GOTO  1111 

ELSE 

CALL  CLS 
GOTO  305 
ENDIF 
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ELSEIFCSELECT.  EQ.  '2' )  THEN 

306  WRITE(*,307)  PC21*100 , PC22*100 

307  FORMATC/////’  YOU  HAVE  SELECTED;  ' ,F6.  2 COVERAGE ,  LAYER  1 ' 

+  ,/,25X,F6.  2, COVERAGE,  LAYER  2.'  Jj,'  IS  THIS  YOUR  DESIRED 

+  ENTRY?  (Y  OR  N):  ' ,1X,  ) 

READ  (*,308)ANS 

308  FORMATC Al) 

IFCANS. EQ. 'Y' )  THEN 
GOTO  333 

ELSEIFCANS.EQ. 'N' )  THEN 

GOTO  nil 

ELSE 

CALL  CLS 
GOTO  306 
ENDIF 

ELSEIFCSELECT.  EQ.  '3')  THEN 

310  WRITEC*,311)  PC31*100,PC32*100,PC33*100 

311  FORMATC/////’  YOU  HAVE  SELECTED:  ',F6.2,'%  COVERAGE,  LAYER  l' 

+  ,/,25X,F6.  2, ’%  COVERAGE,  LAYER  2 ' ,/ , 25X,F6.  2 , ' %  COVERAGE,  LAYER  3. 
+',//,’  IS  THIS  YOUR  DESIRED  ENTRY?  CY  OR  N):  ',1X,  ) 

READ  C*,312)ANS 

312  FORMATCAl) 

IFCANS. EQ.  'Y' )  THEN 
GOTO  333 

ELSEIFCANS.EQ.  'N' )  THEN 

GOTO  nil 

ELSE 

CALL  CLS 
GOTO  310 
ENDIF 

ELSEIFCSELECT.  EQ.  '4' )  THEN 

320  VRITEC*,321)  PC41*100,PC42*100,PC43*100,PC44*100 

321  FORMATC/////'  YOU  HAVE  SELECTED:  ',F6.2,’%  COVERAGE,  LAYER  l’ 

+  ,/,25X,F6.  2, ’%  COVERAGE,  LAYER  2' ,/,25X,F6.  2, '%  COVERAGE,  LAYER  3' 
+  ,/,25X,F6.  2, '%  COVERAGE,  LAYER  4. ',//,’  IS  THIS  YOUR  DESIRED 

+  ENTRY?  CY  OR  N):  ' ,1X,  ) 

READ  C*,322)ANS 

322  FORMATCAl) 

IFCANS. EQ. 'Y' )  THEN 
GOTO  333 

ELSEIFCANS.EQ.  'N' )  THEN 

GOTO  nil 

ELSE 

CALL  CLS 
GOTO  320 
ENDIF 
ENDIF 

333  CONTINUE 

C 

C 

c 

C  GENERATE  CONSTANTS  FOR  THE  EPOXY  LAYERS 
C 

C  LEFT  OF  RIGHT  EDGE  TO  OUTSIDE 

C 

C 
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nononnon  nononon  nonoooooooon  ooonnnoo  nnnnon  onoon 


EYLR  =  2*XUKE*DELX*XDELE/DELY 


FRONT  OR  BACK  TO  OUTSIDE 
EFXB  =  2*XUKE*DELY*XDELE/DELX 

INNER  MATRIX  MOVEMENT  IN  THE  Y  DIRECTION 
EYY  =  XUKE*DELX*XDELE/DELY 

INNTR  MATRIX  MOVEMENT  IN  THE  X  DIRECTION 

EXX  =  XUKE*XDELE*DELY/DELX 

EPOXY  TO  COPPER 

ONE  COPPER  LAYER  CASE 

IF( SELECT.  EQ.  ' 1' )  THEN 

EZC11=  2*DELX*DELY/((XDELE/XUKE)+(X1L1*PC11/XUKC)) 
TWO  COPPER  LAYERS  CASE 

ELSEIFCSELECT.  EQ.  '2' )  THEN 

EZC21  =  2*DELX*DELY/((XDELE/XUKE)+(X2L1*PC21/XUKC)) 
EZC22  =  2*DELX*DELY/((XDELE/XUKE)+(X2L2*PC22/XUKC)) 

THREE  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ. '3' )  THEN 
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nnnnonn  o  onoonnoooooooo  nnono  noonooo 


EZC31  =  2*DELX*DELY/((XDELE/XUKE)+(X3L1*PC31/XUKC)) 
EZC32  =  2*DELX*DELY/((XDELE/XUKE)+(X3L2*PC32/XUKC)) 
EZC33  =  2*DELX*DELY/((XDELE/XUKE)+(X3L3*PC33/XUKC)) 


FOUR  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ.  '4')  THEN 

EZC41  =  2*DELX*DELY/((XDELE/XUKE)+(X4L1*PC41/XUKC)) 
EZC42  =  2*DELX*DELY/((XDELE/XUKE)+(X4L2*PC42/XUKC)) 
EZC43  =  2*DELX*DELY/((XDELE/XUKE)+(X4L3*PC43/XUKC)) 
EZC44  =  2*DELX*DELY/((XDELE/XUKE)+(X4L4*PC44/XUKC)) 


ELSE 

CONTINUE 

ENDIF 


EPOXY  TO  BOTTOM  OUTER  EDGE 
EZB  =  2*XUKE*DELX*DELY/XDELE 


GENERATE  CONSTANTS  FOR  THE  COPPER  LAYERS 


CCCCCCCCCCCCCCC 


LEFT  OF  RIGHT  EDGE  TO  OUTSIDE 

ONE  COPPER  LAYER  CASE 


IF( SELECT. EQ.  ' r )  THEN 
CYLR11=  2*XUKC*X1L1*PC11/DELY 


TWO  COPPER  LAYERS  CASE 


ELSEIFCSELECT.  EQ.  '2' )  THEN 
CYLR21=  2*XUKC*X2L1*PC21/DELY 
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non  nnnnnn  nnnnonnnnn  nonnnnn  nnnonon 


CYLR22=  2*XUKC*X2L2*PC22/DELY 


THREE  COPPER  LAYERS  CASE 


ELSEIF(SELECT.EQ. '3' )  THEN 
CYLR31=  2*XUKC*X3L1*PC31/DELY 
CYLR32=  2*XUKC*X3L2*PC32/DELY 
CYLR33=  2*XUKC*X3L3*PC33/DELY 


FOUR  COPPER  LAYERS  CASE 


ELSEIFC SELECT. EQ. '4' )  THEN 
CYLR41=  2*XUKC*X4L1*PC41/DELY 
CYLR42=  2*XUKC*X4L2*PC42/DELY 
CYLR43=  2*XUKC*X4L3*PC43/DELY 
CYLR44=  2*XUKC*X4L4*PC44/DELY 

ELSE 
CONTINUE 
END  IF 


FRONT  OR  BACK  TO  OUTSIDE 


ONE  COPPER  LAYER  CASE 


IF(  SELECT.  EQ.  '  r  )  THEN 

'1FXB11=  2*XUKC*DELY*X1L1*PC11/DELX 


TWO  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ. '2' )  THEN 
CFXB21=  2*XUKC*DELY*X2L1*PC21/DELX 
CFXB22=  2*XUKC*DELY*X2L2*PC22/DELX 


THREE  COPPER  LAYERS  CASE 
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noonnno  nnoooon  oooooonn  oononoo  onoo 


ELSEIFCSELECT.  EQ.  '3' )  THEN 
CFXB31=  2*XUKC*DELY*X3L1*PC31/DELX 
CFXB32=  2*XUKC*DELY*X3L2*PC32/DELX 
CFXB33=  2*XUKC*DELY*X3L3*PC33/DELX 


FOUR  COPPER  LAYERS  CASE 


ELSEIFCSELECT.  EQ.  '4' )  THEN 
CFXB41=  2*XUKC*DELY*X4L1*PC41/DELX 
CFXB42=  2*XUKC*DELY*X4L2*PC42/DELX 
CFXB43=  2*XUKC*DELY*X4L3*PC43/DELX 
CFXB44=  2*XUKC*DELY*X4L4*PC44/DELX 

ELSE 

CONTINUE 

ENDIF 

INNER  MATRIX  MOVEMENT  IN  THE  Y  DIRECTION 


ONE  COPPER  LAYER  CASE 


IF( SELECT.  EQ.  ' l' )  THEN 
CYY11=  XUKC*DELX*X1L1*PC11/DELY 


TWO  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ. '2' )  THEN 
CYY21=  XUKC*DELX*X2L1*PC21/DELY 
CYY22=  XUKC*DELX*X2L2*PC22/DELY 


THREE  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ. '3' )  THEN 
CYY31=  XUT(C*DELX*X3L1*PC31/DELY 
CYY32=  XUKC*DETJ(*X3L2*PC32/DELY 
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ooo  ooonnnn  oonoDoo  onnoooonono  oooonoo 


CyY33=  XUKC*DELX*X3L3*PC33/DELY 


FOUR  COPPER  LAYERS  CASE 


ELSEIFCSELECT.  EQ.  '4' )  THEN 
CYY41=  XUKC*DELX*X4Ll*PC4l/DELY 
CYY42=  XUKC*DELX*X4L2*PC42/DELY 
CYY43=  XUKC*DELX*X4L3*PC43/DELY 
CYY44=  XUKC*DELX*X4L4*PC44/DELY 


ELSE 

CONTINUE 

ENDIF 


INNER  MATRIX  MOVEMENT  IN  THE  X  DIRECTION 


ONE  COPPER  LAYER  CASE 


IF( SELECT. EQ. ' 1' )  THEN 
CXX11=  XUKC*X1L1*PC11*DELY/DELX 


TWO  COPPER  LAYERS  CASE 


ELSEIFCSELECT.  EQ. '2' )  THEN 
CXX21=  XUKC*X2L1*PC21*DELY/DELX 
CXX22=  XUKC*X?L2*PC22*DELY/DELX 


THREE  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ. '3' )  THEN 
CXX31=  XUKC*X3L1*PC31*DELY/DELX 
CXX32=  XUKC*X3L2*PC32*DELY/DELX 
CXX33=  XUKC*X3L3*PC33*DELY/DELX 


FOUR  COPPER  LAYERS  CASE 
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onoooon  non  o  non  nooonoo  ocsooonnoo  oono 


ELSEIF(SELECT.  EQ.  '4')  THEN 
CXX41=  XUKC*X4L1*PC41*DELY/DELX 
CXX42=  XUKC*X4L2*PC42*DELY/DELX 
CXX43=  XUKC*X4L3*PC43*DELY/DELX 
CXX44=  XUKC*X4L4*PC44*DELY/DELX 


ELSE 

CONTINUE 

ENDIF 


COPPER  TO  EPOXY  (OR  AIR) 


ONE  COPPER  LAYER  CASE 


IF( SELECT.  EQ.  ' 1' )  THEN 

CZE11=  2*XUKC*DELX*DELY/(X1L1*PC11) 


TWO  COPPER  LAYERS  CASE 


ELSEIFCSELECT. EQ.  '2')  THEN 

CZE21  =  2*DELX*DELY/((XDELE/XUKE)+(X2L1*PC21/XUKC)) 
CZE22  =  2*XUKC*DELX*DELY/(X2L2*PC22) 


THREE  COPPER  LAYER  CASE 


ELSEIFCSELECT.  EQ.  '3' )  THEN 

CZE31  =  2*DELX*DELY/((XDELE/XUKE)+(X3L1*PC31/XUKC)) 
CZE32  =  2*DELX*DELY/((XDELE/XUKE)+(X3L2*PC32/XUKC)) 
CZE33  =  2*XUKC*DELX*DELY/(PC33*X3L3) 


FOUR  COPPER  LAYER  CASE 


ELSEIFCSELECT.  EQ.  '4' )  THEN 

CZE41  =  2*DELX*DELY/((XDELE/XUKE)+(X4L1*PC41/XUKC)) 
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o  o  o  n  n  n  o  o  n  o  n  n  n  o  o  o  o  n  rs  n  o  o  o  o  o  o  o  n  n  n  o  o  o 


CZE42  =  2*DELX*DELY/((XDELE/XUKE)+(X4L2*PC42/XUKC)) 
CZE43  =  2*DELX*DELY/((XDELE/XUKE)+(X4L3*PC43/XUKC)) 
CZE44  =  2*XUKC*DELX*DELY/(PC44*X4L4) 

ELSE 

CONTINUE 
END  IF 

CALL  THE  OUTPUT  DATA  FILE  SUBROUTINES 


ONE  COPPER  LAYER  CASE 


IF( SELECT.  EQ.  ' 1' )  THEN 

CALL  S1(EYLR,EFXB,EYY,EXX,EZC11,EZB,CYLR11,CFXB11,CYY11,CXX11,CZE1 
+1 , XPL , XIH , XJH , XEAT , XWIDE , XDEEP , XOEF ,XIBT , XPRT , XWRT , XRT , XLT , XFT , XBT 
+,USEL) 


TWO  COPPER  LAYER  CASE 


ELSEIF(SELECT.EQ.  '2' )  THEN 

CALL  S2 ( EYLR , EFXB , EYY , EXX , EZC2 1 ,EZC22 , EZB , CYLR2 1 , CYLR22 , CFXB2 1 , CFX 
+B 2 2 , C YY2 1 , C YY2 2 , CXX2 1 , CXX22 , CZE2 1 , CZE22 , XPL , XIH , XJH , XEAT , XWIDE , XDE 
+EP , XOEF , XI BT , XPRT , XWRT , XRT , XLT ,XFT , XBT , USEL) 


THREE  COPPER  LAYER  CASE 


ELSEIFCSELECT.  EQ.  '3' )  THEN 

CALL  S3(EYLR, EFXB, EYY, EXX, EZC31,EZC32,EZC33, EZB, CYLR31,CYLR32,CYLR 
+33,CFXB31,CFXB32,CFXB33,CYY31,CYY32,CYY33,CXX31,CXX32,CXX33,CZE31, 
+CZE32,CZE33, XPL, XIH, XJH, XEAT, XWIDE, XDEEP, XOEF, XIBT, XPRT, XWRT, XRT, X 
+LT,XFT, XBT, USEL) 


FOUR  COPPER  LAYER  CASE 


ELSEIFCSELECT.  EQ.  '4' )  THEN 

CALL  S4( EYLR , EFXB , EYY , EXX , EZC4 1 , EZC42 , EZC43 , EZC44 , EZB , CYLR4 1 , CYLR4 
+2 ,CYLR43 ,CYLR44 ,CFXB41 ,CFXB42 ,CFXB43 ,CFXB44 ,CYY41 ,CYY42 ,CYY43 ,CYY4 
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+4 , CXX4 1 , CXX42 , CXX43 , CXX44 , CZL4 1 , C2E42 , CZE43 , CZE44 , XPL , XIH , XJH , XEAT 
+ ,  XW  IDE ,  XDEEP ,  XOEF ,  XI BT ,  XPRT ,  XVv'RT ,  XRT ,  XLT ,  XFT ,  XBT ,  USEL) 

ELSE 

CONTINUE 
END  IF 
END 

SUBROUTINE  COPPER(EW, EL, ANSN, UK, SELECT, T1L1,T2L1,T2L2,T3L1,T3L2,T3 
+L3 , T4L1 , T4L2 , T4L3 , T4L4 , SPEVAL , UKC , CL , CW) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


TITLE: 

AUTHOR: 

DATE: 

COMPILER: 

LINKER: 


MODEL  BUILDER 
LT  STEVE  GLASER 
09  JUL  1991 

MICROSOFT  VERSION  4.  01 
MICROSOFT  VERSION  3.55 


DEFINE  REAL  VARIABLES 

SUBROUTINE  COPPER  GATHERS  ALL  THE  INFORMATION  REQUIRED  FOR  THE  COPPER 
LAYERS. 

REAL  EL,EW ,UN ,T1L1 ,T2L1 ,T2L2 ,T3L1 ,T3L2 ,T3L3 ,T4L1 ,T4L2 ,T4L3 ,T4L4 , VO 
+L 1 , V0L2 , V0L3 , V0L4 , TIW 1 , T2W 1 , T2W2 , T3W 1 , T3W2 , T3W3 , T4W 1 , T4W2 , T4W3 , T4W 
+4, UKC 


INTEGER  NWIDE,NDEEP,NPL 
C  DEFINE  ALL  ONE  CHARACTER  VARIABLES 
C 

CHARACTER* 1  ANSN , SELECT , WORL , SPEVAL , ANSE 1 , ANSTl , ANSKl 
C  CHARACTER  VARIABLES  OF  MORE  THAN  ONE  POSITION 


CHARACTER  L0C0*2 ,UK*10 ,UT*1 
C 


C  DEFINE  MATRICES 

C 

C 

C 

C 

C 

C  COPPER  LAYER  CHARACTERISTICS 

C 

C 

5458  CALL  CLS 
31171  WRITE(*, 18000) 

18000  FORMAT( III,' ****************************************************** 
+  ' *******iir**********cOPPER  LAYER  CHARACTERISTICS******************* 

+  '  *:fr*TfrTfr'i*nlnfrTV':fr'Jr^r^'******************’^***frAifrA*^*yr***********  ill  <i'i>i 

+******** ' , / , ) 

C 

C  PROVIDE  CORRECT  UNIT  ABBREVIATIONS 
C 


AREA=EW*EL 
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IFCANSN. EQ. 'S' )  THEN 

UN=12*2.  54*12*2.  54*.  00134374*2. 54 
WRITE(*. 15400) 

15400  FORMATC^  ALL  ENTRIES  ARE  IN  SI  NOTATION.  ' ,/) 

LOCO=’cm' 

C  UK='Watts/cm/C' 

C  UT='C' 

ELSEIFCANSN.  EQ.  'E' )  THEN 
UN=12*12*0. 00134374 
WRITE(*, 15401) 

15401  FORMATC  all  ENTRIES  ARE  IN  ENGLISH  NOTATION',/) 

LOCO=' in' 

C  UK='Btu/hr/F' 

C  UT='F' 

ENDIF 

C 

C 

c 

WRITE(*, 15402)  EL, LOCO 

15402  FORMATC/,'  THE  COPPER  LAYER  LENGTH  IS  THE  SAME  AS  THE  EPOXY  LAYER 
+:  ' ,F9.4,1X,A2,2X,  ) 

CL=EL 

C 

V/RITEC*,  15403)  EW.LOCO 

15403  FORMATC//,'  THE  COPPER  LAYER  WIDTH  IS  THE  SAME  AS  THE  EPOXY  LAYER 
+;  ' ,F9.4,1X,A2,2X,  ) 

CW=EW 

14337  WRTTEC*, 16403)  SELECT 

16403  FORMATC//,'  YOU  SELECTED  ',A1,'  COPPER  LAYERCS)  FOR  THE  PCB.  ',/,' 
+  YOU  NOW  WILL  BE  ASKED  TO  ENTER  THE  COPPER  LAYER  THICKNESS' ,2X,  ) 
99999  WRITEC*, 31112) 

31112  FORMATC//,'  DO  YOU  WANT  TO  SPECIFY  THICKNESS  BY  LENGTH  OR  WEIGHT? 
+',/,'  ENTER  L  IF  YOU  WANT  TO  ENTER  LENGTH,  W  IF  YOU  WANT  TO  ENTER 
+  WEIGHT:  ' ,2X,  ) 

READC*, 31113)  WORL 

31113  FORMATC Al) 

IFCWORL. EQ.  ' L' .  AND.  SELECT.  EQ.  ' l’ )  THEN 
CALL  CLS 

WRITEC*, 31114)  LOCO 

31114  FORMATC////,'  ENTER  THE  THICKNESS  FOR  LAYER  1  C’,A2,'):  ' ,2X,  ) 
READ  *,T1L1 

GOTO  31170 

ELSEIFCWORL.  EQ.  ' L' .  ANT.  SELECT.  EQ.  '2’)THEN 
CALL  CLS 

WRITEC*, 31115)  LOCO 

31115  FORMATC////,’  E.NTER  THE  THICKNESS  FOR  LAYER  1  C',A2,'):  '  ,2X,  ) 
READ  *,T2L1 

WRITEC*,31116)  LOCO 

31116  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  2  C',A2,'):  ' ,2X,  ) 
READ  *,T2L2 

GOTO  31170 

ELSEIFCWORL. EQ.  ' L' .  AND.  SELECT.  EQ.  ’ 3' )THEN 
CALL  CLS 

WRITEC*, 31117)  LOCO 

31117  FORMATC////,'  ENTER  THE  THICKNESS  FOR  LAYER  1  C',A2,’):  ' ,2X,  ) 
READ  *,T3L1 


WRITE(*. 31118)  LOCO 

31118  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  2  (’,A2,'):  ',2X,  ) 
READ  *,T3L2 

WRITEC*, 31119)  LOCO 

31119  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  3  (',A2,'):  ’ ,2X,  ) 
READ  *,T3L3 

GOTO  31170 

ELSEIFCWORL. EQ.  ' L' .  AND.  SELECT.  EQ.  ’4')THEN 
CALL  CLS 

VRITEC*. 31120)  LOCO 

31120  FORMATC////,'  ENTER  THE  THICKNESS  FOR  LAYER  1  C',A2,'):  ',2X,  ) 
READ  *,T4L1 

WRITEC*,31121)  LOCO 

31121  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  2  C',A2,'):  ' ,2X,  ) 
READ  *,T4L2 

WRITEC*, 31122)  LOCO 

31122  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  3  C ’ ,A2, ' ):  ' ,2X,  ) 
READ  *,T4L3 

WRITEC*, 31123)  LOCO 

31123  FORMATC/,'  ENTER  THE  THICKNESS  FOR  LAYER  4  C'.A2,'):  ' ,2X,  ) 
READ  *,T4L4 

GOTO  31170 

ELSEIFCWORL.  EQ.  'W' .  AND.  SELECT.  EQ.  '  l'  )  THEN 

CALL  CLS 

W'RITEC*,31142) 

31142  FORMATC////,'  ENTER  THE  WEIGHT  FOR  LAYER  1  Coz):  ',2X,  ) 

IFCUW. EQ. 'gm' )  THEN 

READ  *,T1W1 
T1W1=T1W1*62.  5 
ELSE 

READ  *,T1W1 
ENDIF 

V0L1=UN*T1W1 
T1L1=V0L1/AREA 
GOTO  31170 

ELSEIFCWORL. EQ.  'W' .  AND.  SELECT.  EQ.  '2')THEN 
CALL  CLS 
WRITEC*, 31143) 

31143  FORMATC////,'  ENTER  THE  WTIGHT  FOR  LAYER  1  Coz):  ' ,2X,  ) 

IFCUW. EQ. 'gm' )  THEN 

READ  *,T2W1 
T2W1=T2W1*62.  5 
ELSE 

READ  *,T2W1 
ENDIF 

V0L1=UN*T2W1 
T2L1=V0L1/AREA 
WRITEC *,31 144) 

31144  FORMATC/,'  ENTER  THE  WEIGHT  FOR  LAYER  2  Coz):  ',2X,  ) 

IFCUW. EQ. 'gm' )  THEN 

READ  *,T2W2 
T2W2=T2W2*62.  5 
ELSE 

READ  *,T2W2 
ENDIF 

VOL2=UN*T2W2 
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T2L2=V0L2/AREA 
GOTO  31170 

ELSEIF(WORL. EQ.  'W' .  AND.  SELECT. EQ.  '3')THEN 
CALL  CLS 
WRITE (*,31 145) 

31145  FOR.MATC////,  '  E.NTER  THE  WEIGHT  FOR  LAYER  1  (oz):  '  ,2X,  ) 
IF(UW'.  EQ.  'gm'  )  THEN 

READ  *,T3W1 
T3W1=T3W1*62,  5 
ELSE 

READ  *,T3W1 
ENDIF 

V0L1=UN*T3W1 
T3L1=V0L1/AREA 
WRITE(*, 31146) 

31146  FORMATC/,'  ENTER  THE  WEIGHT  FOR  LAYER  2  (oz):  ' ,2X,  ) 
IF(UW. EQ.  'gm' )  THEN 

READ  *,T3W2 
T3W2=T3W2*62.  5 

FT 

READ  *,T3W2 
ENDIF 

VOL2=UN*T3W2 

T3L2=VOL2/AREA 

WRITE(*,31147) 

31147  FORMATC/,'  ENTER  "niE  WEIGHT  FOR  LAYER  3  (oz):  ’ ,2X,  ) 
IF(UW. EQ. 'gm' )  THEN 

READ  *,T3W3 
T3W3=T3W3*62.  5 
ELSE 

READ  *,T3W3 
ENDIF 

VOL3=UN*T3W3 
T3L3=VOL3/AREA 
GOTO  31170 

ELSE IF(  WORL.  EQ.  '  W ' .  ANT).  SELECT.  EQ.  '  4 '  )THEN 
CALL  CLS 
WRITEC*, 31148) 

31148  FORMATC////,'  ENTER  THE  WEIGHT  FOR  LAYER  1  (oz):  ' ,2X,  ) 
IFCirW.  EQ.  'gm'  )  THEN 

READ  *,T4W1 
T4W1=T4W1*62.  5 
ELSE 

READ  *,T4W1 
ENDIF 

V0L1=UN*T4W1 
T4L1=V0L1/AREA 
WRITEC*, 31149) 

31149  FORMATC/,'  ENTER  THE  WEIGHT  FOR  LAYER  2  (oz):  ' ,2X,  ) 
IFCUW.  EQ.  'gm' )  THEN 

READ  *,T4W2 
T4W2=T4W2*62.  5 
ELSE 

READ  *,T4W2 
ENDIF 

VOL2=UN*T4W2 


T4L2=V0L2/AREA 
WRITEC*, 31150) 

31150  FORMATC/,'  ENTER  THE  WEIGHT  FOR  LAYER  3  (oz):  ',2X,  ) 

IF(UV. EQ. 'gm' )  THEN 

READ  *,T4W3 
T4W3=T4W3*62.  5 
ELSE 

READ  *,T4V3 
ENDIF 

V0L3=UN*T4W3 
T4L3=V0L3/AREA 
WRITE(*, 31151) 

31151  FORMATC/,'  ENTER  THE  WEIGHT  FOR  LAYER  4  (oz):  ' ,2X,  ) 

IFCUW. EQ. 'gm')  THEN 

READ  *,T4W4 
T4W4=T4W4*62.  5 
ELSE 

READ  *,T4W4 
ENDIF 

V0L4=UN*T4W4 
T4L4=V0L4/AREA 
GOTO  31170 
CALL  CLS 
GOTO  31170 
ELSE 

CALL  CLS 
GOTO  31171 
ENDIF 

31170  IFCSPEVAL.EQ. 'B')  THEN 
GOTO  41186 
ELSE 

CONTINUE 

ENDIF 

WRITEC*, 31172)  UK 

31172  FORMATC/,'  ENTER  COPPER  LAYER  THERMAL  CONDUCTIVITY  (',A10,'): 

+2X,  ) 

READ  *,UKC 
C 

C  MAKE  CHANGES  OR  CORRECTIONS  TO  COPPER  LAYER  ENTRIES 
C 

31193  IF( SELECT.  EQ.  ' 1' )THEN 
CALL  CLS 

WRITEC*, 31173)  CL, LOCO 

31173  FORMATC///,'  YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  COPPE 
+R  LAYER(S).  ' ,///, 

+  '  1.)  LENGTH:  ' ,F9. 4, 1X,A2,2X,  ) 
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WRITEC*, 73)  CW.LOCO 

FORMATC/,'  2.)  WIDTH: 

’ ,F9.4,1X,A2,2X,  ) 

74 

WRITEC*, 74)  TlLl.LOCO 

FORMATC/,'  3.)  THICKNESS  LAYER 

1;  ' ,F9.  4,1X,A2,2X,  ) 

75 

WRITEC*, 75)  UKC,UK 

FORMATC/,'  4.  )  k: 

' ,F9. 4,1X,A10,2X,  ) 

311 

ELSEIFC  SELECT.  EQ.  ' 2 ' )THEN 

CALL  CLS 

WRITEC*, 311)  CL, LOCO 

FORMATC///,'  YOU  HAVE  MADE  THE 

FOLLOWING  ENTRIES  FOR  THE  COPPE 
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+R  LAYER(S). ' ,///, 

+'  1.)  LENGTH:  ' ,F9. 4, 1X,A2 ,2X,  ) 

WRITE(*,312)  CW.LOCO 

312  FORMATC/,'  2.  )  WIDTH:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITE(*,314)  T2Ll,LOCO 

314  FORMATC/,'  3.)  THICKNESS  LAYER  1:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITE(*,315)  T2L2,LOCO 

315  FORMATC/,'  THICKNEbS  LAYER  2:  ' ,F9. 4, IX, A2 , 2X,  ) 

WRITEC*,316)  UKC,UK 

316  FORMATC/,'  4.)  k:  ' ,F9. 4, 1X,A10,2X,  ) 

ELSEIFC  SELECT.  EQ.  ' 3 ' )THEN 

CALL  CLS 

WRITEC*,317)  CL, LOCO 

317  FORMATC///,’  YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  COPPE 
+R  LAYERCS).  ',///, 

+'  1.)  LENGTH:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,318)  CW,LOCO 

318  FORMATC/,'  2.  )  WIDTH:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,319)  T3Ll,LOCO 

319  FORMATC/,'  3.)  THICKNESS  LAYER  1:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,320)  T3L2,LOCO 

320  FORMATC/,'  THICKNESS  LAYER  2:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,321)  T3L3,LOCO 

321  FORMATC/,'  THICKNESS  LAYER  3:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,322)  UKC,UK 

322  FORMATC/,'  4.  )  k:  ' ,F9. 4, 1X,A10,2X,  ) 

ELSEIFC SELECT. EQ.  '4')THEN 

CALL  CLS 

TOITEC*,323)  CL, LOCO 

323  FORMATC///,'  YOU  HAVE  MADE  THE  FOLLOWING  ENTRIES  FOR  THE  COPPE 
+R  LAYERCS).  ' ,///, 

+'  1.)  LENGTH:  ’ ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,324)  CW,LOCO 

324  FORMATC/,'  2.  )  WIDTH:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,325)  T4Ll,LrC0 

325  FORMATC/,'  3.)  THICKNESS  LAYER  1:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,326)  T4L2,L0C0 

326  FORMATC/,'  THICKNESS  LAYER  2:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,327)  T4L3,L0C0 

327  FORMATC/,'  THICKNESS  LAYER  3:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,328)  T4L4,L0C0 

328  FORMATC/,'  THICKNESS  LAYER  4:  ' ,F9. 4, 1X,A2,2X,  ) 

WRITEC*,329)  UKC,UK 

329  FORMATC/,'  4.  )  k:  ' ,F9. 4, 1X,A10,2X,  ) 

ENDIF 

C 

31194  WRITEC*, 71180) 

71180  FORMATC///,'  DO  YOU  WISH  TO  MAKE  ANY  CHANGES?  SELECT  Y  FOR  Y 

+ES  AND  N  FOR  NO:  ',2X,  ) 

READC*,31181)ANSE1 
31181  FORMATC Al) 

IFCANSEl.EQ. 'Y' )  THEN 
56342  WRITEC*, 31184) 

31184  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  THICKNESS?  CY  OR  N 
+):  ',2X,  ) 


READ(*,31185)ANST1 
31185  FORMAT(Al) 

PRINT  * 

IF(AN3T1.EQ. 'Y' )  THEN 
SPEVAL='B' 

CALL  CL= 

WRITE(*,!J9998) 

99998  FORMATC////) 

GOTO  99999 

ELSEIF(ANST1.EQ. 'N')  THEN 

CONTINUE 

ELSE 

CALL  CLS 
WRITE(*, 56343) 

56343  FORMATC////) 

GOTO  56342 
ENDIF 


C 

41186  CALL  CLS 

WRITEC*, 31187) 

31187  FORMATC////) 

PRINT  THE  CURRENT  ENTRY  FOR  THERMAL  CONDUCTIVITY  IS  ' ,UK 

+C  '  '  UK 

32186  ’  WRITEC*, 31188) 

31188  FORMATC/,'  WOULD  YOU  LIKE  TO  CHANGE  THE  THERMAL  CONDUCTIVIT 

+Y?  CY  OR  N):  ' ,2X,  ) 

READC*,31189)ANSK1 

31189  FORMATCAl) 

PRINT  * 

IFCANSKl.EQ.  'Y')  THEN 
WRITEC*, 31190)  UK 

31190  FORMATC/,'  ENTER  THE  THERMAL  CONDUCTIVITY  C ’ ,A10, ' ):  ' ,2X,  ) 

READ  *,UKC 

ELSEIFCANSKl.EQ.  'N')  THEN 

GOTO  31191 

ELSE 

CALL  CLS 
WRITEC*, 56344) 

56344  FORMATC////) 

GOTO  32186 
ENDIF 

31191  GOTO  31193 
C 

C 

ELSEIFCANSEl.EQ.  'N' )  THEN 

GOTO  31180 

ELSE 

CALL  CLS 
GOTO  31194 
ENDIF 
31180  END 


C 

C  TITLE:  MODEL  BUILDER 

C  SUBROUTINE:  SI 

C  DATE:  09  JUL  91 
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AUTHOR:  LT  STEVE  GLASER 

COMPILER:  MICROSOFT  VERSION  4. 01 

LINKER:  MICROSOFT  VERSION  3.55 

CALLED  FROM  SUBROUTINE  PCBS3.  THIS  SUBROUTINE  GENERATES  THE  OUTPUT 

DATA  FILE  FOR  THE  ONE  COPPER  LAYER  CASE 


SUBROUTINE  SI  (GYLR.GFXB ,GYY,GXX,GZC11 ,GZB .GYLRll .GFXBll ,GYY11 ,GXX 
+11,GZE11,GPL,GIH,GJH,GEAT,GWIDE,GDEEP,G0EF,GIBT,GPRT,GWRT,GRT,GLT, 
+GFT,GBT,GSEL) 

REAL  GYLR,GFXB,GYY,GXX,GZC11,GZB,GYLR11,GFXB11,GYY11,GXX11,GZE11,A 
-J-CC ,  DAMP ,  CONFAC ,  G I BT ,  GPRT ,  GWRT ,  GRT ,  GLT ,  GFT ,  GBT 
REAL  GEAT(100,100),GOEF(1000,12) 


CHARACTER*!  SELECT, ANS 
CHARACTER  DATAF*79 ,NAME*6 

INTEGER  GPL , I , N , IB , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS,Dl,D2,D3,D4,D5,D6,D7,MAXIT,LOCVAR 

INTEGER  GIH( 1000) ,GJH( 1000) ,GCON( 1000,12) 

C0NTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0.  8 
C0UNT=2*GPL 
LOCVAR=COUNT 


898  CALL  CLS 


WRITE(*,803) 

803  FORMATC///,’  THIS  PROGRAM  CREATES  AN  OUTPUT  DATA  FILE  FOR  ENTRY 
+  INTO  THE ' , / , 

+'  EXISTING  THERMAL  ANALYZER,  FURTHERMORE,  THIS  PROGRAM  DOES',/, 

+'  NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE.  THEREFORE’,/, 

+'  THE  USER  WILL  NAME  THE  DATA  FILE  FOR  EACH  RUN  OF  THIS  ' ,/, 

+  '  PROGRAM.  THE  FILE  NA.ME  IS  LIMITED  TO  SIX  CHARACTERS,  AND',/, 

+'  SHOULD  NOT  HAVE  ANY  SPACES.  ',///, 
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+’  PLEASE  ENTER  THE  DESIRED  DATA  FILE  NAME:  ',2X  ) 

READ(*,804)  NAME 

804  FORMAT(A6) 

552  WRITE(*,910)NAME 

910  FORMATC /////,'  YOU  SELECTED  ’,A6,'  FOR  YOUR  DATA  FILE  NAME 
+  ’  ,//) 

911  WRITE(*,811) 

811  FORMATC'  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N 

+FOR  NO:  ' ,2X,  ) 

READ(*,812)  ANS 

812  FORMAT(Al) 

C 

IFCANS.EQ.  'N' )  THEN 
GOTO  898 
ELSE 

CONTINUE 

ENDIF 

IFCANS.EQ.  'Y' )  THEN 
GOTO  897 
ELSE 

CALL  CLS 
GOTO  552 
897  ENDIF 

C  ALLOW  THE  USER  TO  PROVIDE  TITLE  LINE  FOR  DATA  FILE 

833  CALL  CLS 
WRITEC*,805) 

805  FORMATC////,'  ENTER  THE  DESIRED  TITLE  TO  BE  PLACED  ON  LINE 

+’,/,'  NUMBER  ONE  OF  THE  OUTPUT  DATA  FILE: ',///, 

+'  '  2X  ) 

READC*,806)’DAfAF 

806  FORMATC A79) 

835  WRITEC*,831) 

831  FORMATC////,'  DO  YOU  WISH  TO  CHANGE  THE  TITLE  OF  YOUR  OUTPU 

+T  DATA  FILE?' ,/,'  ENTER  Y  FOR  YES  AND  N  FOR  NO:  ’ ,2X,  ) 

READC*,832)  ANS 

832  FORMATC Al) 

C 

IFCANS. EQ.  'Y' )  THEN 
GOTO  833 
ELSE 

CONTINUE 

ENDIF 

IFCANS.EQ. 'N')  THEN 
GOTO  834 
ELSE 

GOTO  835 

834  ENDIF 


C**************c0EFFICIENTS  for  epoxy  and  COPPER  LAYERS***************** 

C*********************************************************************** 


C 


DO  90  1=1, GPL 
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N=1 

IB=GPL+I 

C 

C 

C*********************************C0RNERS******************************** 

C 

IF  ((GIH(I).EQ.  l.OR.  GIH(I).EQ.  GDEEP).AND.  (GJH(I).EQ.  1.  OR.  GJH(  I).  EQ 
+. GWIDE))  THEN 

DETERMINE  COEFFICIENTS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).EQ.  0.0)  THEN 
GC0N(I,N)  =  6 

ELSE 

GCONd.N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
LEFT  EDGE 

IF  (GJH(I). EQ.  1)  THEN 
LEFT  COEFFICIENT 
G0EF(I,N)  =  GYLRll 
G0EF(IB,N)  =  GYLR 
N=N+1 

GCONd.N)  =  7551 
GCONCIB.N)  =  7551 

RIGHT  COEFFICIENT 
GOEFd.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCONd.N)  =  10  *  (1+1)  +  1 
GCONCIB.N)  =  10*  (IB+l)  +  1 

RIGHT  EDGE 

ELSE  IF  (GJHd).EQ.  GWIDE)  THEN 

LEFT  COEFFICIENT 

GOEFd.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCON(I.N)  =  10*(I-1)+1 
GCONCIB.N)  =  10*(IB-1)+1 

RIGHT  COEFFICIENT 

GOEF(I.N)  =  GYLRll 
GOEFCIB.N)  =  GYLR 
N=N+1 

GCON(I.N)  =  7541 
GCONCIB.N)  =  7541 
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END  IF 


FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
FRONT  EDGE 

IF  (GIH(I). EQ.  1)  THEN 
FRONT  COEFFICIENT 
G0EF(I,N)=GFXB11 
GOEF(IB,N)=GFXB 
N=N+1 

GCONd.N)  =7521 
GC0N(IB,N)  =  7521 

BACK  COEFFICIENT 
G0EF(I,N)=GXX11 
GOEF(IB,N)=GXX 
N=N+1 

GC0N(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 

BACK  EDGE 

ELSEIF  (GIH(I).EQ. GDEEP)  THEN 
FRONT  COEFFICIENT 
G0EF(I,N)=GXX11 
GOEF(IB,N)=GXX 
N=N+1 

GCONd.N)  =10*d-GWIDE)+l 
GCONCIB.N)  =  10*(IB-GWIDE)+1 

BACK  COEFFICIENT 
GOEFd,N)=GFXBll 
GOEFdB,N)=GFXB 
N=N+1 

GCONd.N)  =  7511 
GCONCIB.N)  =  7511 
END  IF 

TOP  COEFFICIENT 
GOEFd.N)=GZEll 
GOEFdB.N)=GZCll 
N=N+1 

GCONd.N)  =  7511 
GCONCIB.N)  =  10*1+1 

BOTTOM  COEFFICIENT 

GOEFCI.N)  =  GZCll 
GOEFCIB.N)  =  GZB 
N=N+1 

GCONCI.N)  =10*CI+GPL)+1 
GCONCIB.N)  =  10*CIB+GPL)+1 

HEAT  INPUT 

IFCGEATCGIHCD  ,GJHCI)).NE.  0.  0)  THEN 
GOEFCI.N)  =  GEATCGIHCD.GJHCD) 

N=N+1 


95 


ooooo  no  noon  no  on 


GCON(I,N)=9991 

ENDIF 


ELSEIF((GIH(I).EQ.  1.  OR.  GIH(  I).  EQ.  GDEEP).  AND.  (GJH(I).NE.  l.OR.GJH(I) 
+. NE.GWIDE))  THEN 

DETERMINE  NUMBER  OF  CONNECTIONS  FOR  COPPER  AND  EPOXY  LAYERS 
IF  (GEAT(GIH(I),GJH(I)).EQ.  0.0)  THEN 
GC0N(I,N)  =  6 

ELSE 

GCONd.N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS 
LEFT  COEFFICIENT 
G0EF(I,N)  =  GYYll 
G0EF(IB,N)  =  GYY 
N=N+1 

GC0N(I,N)  =  10*(I-1)+1 
GC0N(IB,N)  =  10*(IB-1)+1 

RIGHT  COEFFICIENT 
GOEFd.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCONd.N)  =  10  *  (1+1)  +  1 
GC0NdB,N)  =  10*  (IB+l)  +  1 


FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
FRONT  EDGE 

IF  (GIHCI). EQ. 1)  THEN 
C  FRONT  COEFFICIENT 
GOEFd,N)=GFXBll 
GOEFdB,N)=GFXB 
N=N+1 

GCONd.N)  =7521 
GCONCIB.N)  =  7521 

C  BACK  COEFFTCIENT 

GOEFd,N)=GXXll 
GOEFdB,N)=GXX 
N=N+1 

GCONd.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 

C  BACK  EDGE 

ELSEIF  (GIHCI). EQ. GDEEP)  THEN 
C  FRONT  COEFFICIENT 

GOEFd.N)=GXXll 
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GOEF(IB,N)=GXX 

N=N+1 

GC0N(I,N)  =10*(I-GWIDE)+1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 

C  BACK  COEFFICIENT 
G0EF(I,N)=GFXB11 
GOEF(IB,N)=GFXB 
N=N+1 

GCON(I,N)  =  7511 
GC0N(IB,N)  =  7511 
END  IF 

TOP  COEFFICIENT 
G0EF(I,N)=GZE11 
G0EF(IB,N)=GZC11 
N=N+1 

GC0N(I,N)  =  7511 
GCON(IB,N)  =  10*1+1 

BOTTOM  COEFFICIENT 

G0EF(I,N)  =  GZCll 
G0EF(IB,N)  =  GZB 
N=N+1 

GC0N(I,N)  =10*(I+GPL)+1 
GC0N(IB,N)  =  10*(IB+GPL)+1 

HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0. 0)  THEN 
GOEFd.N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 

GCON(I,N)=9991 
END  IF 


LEFT  AND  RIGHT  EDGES  EXCLUDING  CORNERS 

ELSEIF((GIH(I).  EQ.  l.OR.  GIH(I).EQ.  GDEEP).AND.  (GJH(I).EQ.  1.  OR.  GJH(I) 
+. EQ. GWIDE))  THEN 

DETERMINE  COEFFICIENTS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).NE. 0.0)  THEN 
GC0N(I,N)  =  7 

ELSE 

GCONd.N)  =  6 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GCONCIB.N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
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C  LEFT  EDGE 

IF  (GJH(I). EQ. 1)  THEN 
C  LEFT  COEFFICIENT 

GOEF(I,N)  =  GYLRll 
GOEF(IB,N)  =  GYLR 
N=N+1 

GCONd.N)  =  7551 
GCON(IB,N)  =  7551 
C 

C  RIGHT  COEFFICIENT 

GOEFd.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCONd.N)  =  10  *  (1+1)  +  1 
GCONCIB.N)  =  10*  (IB+l)  +  1 
C 

C  RIGHT  EDGE 

ELSE IF  (GJH(I). EQ. GWIDE)  THEN 
C 

C  LEFT  COEFFICIENT 

C 

GOEFCI.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCONCI.N)  =  10*d-l)+l 
GCONCIB.N)  =  10*dB-l)+l 
C 

C  RIGHT  COEFFICIENT 

C 

GOEFd.N)  =  GYLRll 
GOEFCIB.N)  =  GYLR 
N=N+1 

GCONd.N)  =  7541 
GCONCIB.N)  =  7541 
END  IF 
C 
C 

C  FRONT  COEFFICIENT 

G0EF(I.N)=GXX11 
GOEF(IB.N)=GXX 
N=N+1 

GCONCI.N)  =10*(I-GWIDE)+1 
GCONCIB.N)  =  10*(IB-GWIDE)+1 

C  BACK  COEFFICIENT 

G0EF(I.N)=GXX11 
GOEF(IB.N)=GXX 
N=N+1 

GCONd.N)  =10*(I+GWIDE)+1 
GCONCIB.N)  =  10*(IB+GWIDE)+1 

C 

C  TOP  COEFFICIENT 

G0EF(I,N)=GZE11 
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G0EF(IB,N)=GZC11 

N=N+1 

GCONd.K)  =  7511 
GCON(IB,N)  =  10*1+1 

BOTTOM  COEFFICIENT 

G0EF(I,N)  =  GZCll 
G0EF(IB,N)  =  GZB 
N=N+1 

GC0N(I,N)  =10*(I+GPL)+1 
GC0N(IB,N)  =  10*(IB+GPL)+1 

HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE. 0. 0)  THEN 
GOEF(I.N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 

GC0N(I,N)=9991 

ENDIF 


****DETERMINE  COEFFICIENTS  FOR  ALL  NODES  NOT  TOUCHING  AN  EDGE******* 


ELSEIFC ( GIH( I) .  NE.  1.  OR.  GIH( I ) .  NE.  GDEEP) .  AND.  ( GJH( I ) . NE.  1.  OR.  GJH( I ) 
+.  NE.GWIDE))  THEN 

DETERMINE  CONNECTIONS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I) ,GJH(I)). NE. 0. 0)  THEN 
GCONd.N)  =  7 

ELSE 

GCONCI.N)  =  6 
ENDIF 


CONNECTIONS  FOR  EPOXY  LAYER 
GCONCIB.N)  =  6 


LEFT  COEFFICIENT 
GOEFCI.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCON(I,N)  =  10*d-l)  +  l 
GCONCIB.N)  =  10*dB-l)  +  l 


RIGHT  COEFFICIENT 
GOEFCI.N)  =  GYYll 
GOEFCIB.N)  =  GYY 
N=N+1 

GCONCI.N)  =  10  *  d  +  1)  +  1 
GCONCIB.N)  =  10*  CIB+1)  +  1 
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C  FRONT  COEFFICIENT 
GOEF(I,N)=GXXll 
GOEF(IB,N)=GXX 
N=N+1 

GCONd.N)  =10*(I-GWIDE)  +  1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 

C  BACK  COEFFICIENT 

GOEF(I,N)=GXXll 
GOEF(IB,N)=GXX 
N=N+1 

GCON(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 

C 

C  TOP  COEFFICIENT 

GOEF(I,N)=GZEll 
GOEF(IB,N)=GZCll 
N=N+1 

GCON(I,N)  =  7511 
GCON(IB,N)  =  10*1+1 
C 

C  BOTTOM  COEFFICIENT 

C 

GOEF(I,N)  =  GZCll 
GOEF(IB,N)  =  GZB 
N=N+1 

GCONd.N)  =10*d+GPL)+l 
GCONCIB.N)  =  10*dB+GPL)+l 
C 

C  HEAT  INPUT 

IF(GEAT(GIHd)  .GJHCI)).  NE.  0.  0)  THEN 
GOEFd.N)  =  GEAT(GIHd) , GJHCI)) 

N=N+1 

GCONCI ,N)=9991 

ENDIF 

ENDIF 

90  CONTINUE 

C 

C  GENERATE  DATA  FILE 

C 

OPEN  ( 3, F I LE=NAME,FORM=' FORMATTED' ,ACCESS=' DIRECT' ,RECL=1 08, STATUS 
+='NEW' ) 

WRITE(3,909)  DATAF 
909  FORMATC 1X,A79) 

WRITEC  3 , 908 )  COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908  FORMATC 2X, 9(13, 5X)) 

WRITE(3,907)  ZER, ZER, ZER 
907  FORMATC 2X, 3(13, 5X)) 

WRITEC 3,9081)  NMAX ,TMAX,HTRS ,D1 ,D2 ,D3 ,D4 ,D5 ,D6 ,D7 
9081  FORMATC 2X, 9(1 3, 5X)) 

WR I TE  C  3 , 9  0  5 )  ACC , D AMP , MAX I T , CONFAC , G I BT 

905  F0RMAT(1X,2(F9. 7 , IX) , 14X , 12 , 1X,F9.  7 , 1X,F9.  5) 

WR I TE ( 3 , 9  0  6 )  GPRT , GFT , GBT , GRT , GLT , GWRT 

906  FORMATC 1X,6(F12.  4, IX)) 
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DO  112  I=l,LOCVAR 

WRITE(3,9100)  GCON(I,1),'JCON(I,2),GCON(I,3),UCON(I,4) 

+ON( 1 , 6 ) , GCON( 1,7), GCON( i  ,  8 ) 

9100  FORMAT(I4,3X,7(I5,7X)) 

WRITEC  3,9110)  GOEF( 1,1), GOEF( 1,2), GOEFC 1,3), GOEF( 1,4) 
+EF(I,6),GOEF(I,7) 

9110  FORMAT(7(F9.  3,3X)) 

112  CONTINUE 
CLOSE  (3) 

CALL  CLS 

WRITE(*,999)  NAME 

999  FORMATC////, '  THE  OUTPUT  DATA  HAS  BEEN  PLACED  IN  A  FILE 

+NAMED  ',A6, /////,'  <PRESS  ENTER  TO  CONTINUE>') 

READ(*,5912)  ANS 
5912  FORMAT(Al) 


,GC0NC',5),GC 
,G0EF(I,5) ,G0 


ENT) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


TITLE: 

SUBROUTINE: 

DATE: 

AUTHOR: 

COMPILER: 

LINKER: 


MODEL  BUILDER 
S2 

09  JUL  91 
LT  STEVE  GLASER 
MIC^^OSOFT  VERSION  4.01 
MICROSOFT  VERSION  3.  55 


S2  IS  CALLED  FROM  PCBS3.  S2  IS  THE  OUTPUT  DATA  FILE  GENERATOR  FOR  THE 
THE  COPPER  LAYER  CASE. 


SUBROUTINE  S2  ( GYLR,GFXB ,GYY,GXX,GZC21 ,GZC22 ,GZB ,GYLR21 ,GYLR22 ,GFX 
+B21,GFXB22,GYY21,GYY22,GXX21,GXX22,GZE21,GZE22,GPL,GIH,GJH,GEAT,GW 
+IDE , GDEEP , GOEF , GI BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GSEL) 


REAL  GYLR , GFXB , GYY , GXX , GZC2 1 , GZB , GYLR2 1 , GFXB2 1 , GYY2 1 , GXX2 1 , GZE2 1 , A 
+CC , DAMP , CONFAC , GI BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GYLR2  2 , GFXB  22 , GYY2  2 , 
+GXX22,GZE22,GZC22 


REAL  GEAT(100,100),G0EF(1000,12) 


CHARACTER*!  SELECT, ANS 
CHARACTER  DATAF*79 ,NAME*6 

INTEGER  GPL , I , N , I B , GW I DE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS,Dl,D2,D3,D4,D5,D6,D7,MAXIT,LOCVAR,IC,ID 

INTEGER  GIH( 1000) ,GJH( 1000) ,GCON( 1000,12) 

C0NTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 
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D5=0 

D6=0 

D7=0 

ACC=0. 05 

DAMP=0.  666667 

MAXIT=12 

CONFAC=0.  8 

C0UNT=4*GPL 

L0CVAR=C0UNT 


C 

C 

898  CALL  CLS 

WRITE(*,803) 

803  FORMATC///,'  THIS  PROGRAM  CREATES  AN  OUTPUT  DATA  FILE  FOR  ENTRY 
+  INTO  THE ' , / , 

+'  EXISTING  THERMAL  ANALYZER,  FURTHERMORE,  THIS  PROGRAM  DOES',/, 

+'  NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE.  THEREFORE',/, 

+'  THE  USER  WILL  NAME  THE  DATA  FILE  FOR  EACH  RUN  OF  THIS  ' ,/, 

+'  PROGRAM.  THE  FILE  NAME  IS  LIMITED  TO  SIX  CHARACTERS,  AND',/, 

+'  SHOULD  NOT  HAVE  ANY  SPACES.  ',///, 

+'  PLEASE  ENTER  THE  DESIRED  DATA  FILE  NAME:  ' ,2X  ) 


READC*,804)  NAME 

804 

FORMATC A6) 

552 

WRITEC*,910)NAME 

FORMATC/////,'  YOU  SELECTED  ',A6,'  FOR  YOUR  DATA 

910 

FILE  NAME 

+'.//) 

911 

WRITEC*,811) 

FORMATC’  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y 

811 

FOR  YES  AND  N 

+FOR  NO:  ' ,2X,  ) 

READC*,812)  ANS 

812 

FORMATC Al) 

C 

IFCANS.EQ. 'N')  THEN 
GOTO  898 
ELSE 

CONTINUE 
END  IF 

IFCANS.EQ. 'Y')  THEN 
GOTO  897 
ELSE 

CALL  CLS 
GOTO  552 
897  ENDIF 

C  ALLOW  THE  USER  TO  PROVIDE  TITLE  LINE  FOR  DATA  FILE 

833  CALL  CLS 

WRITE(*,805) 

805  FORMATC////,'  ENTER  THE  DESIRED  TITLE  TO  BE  PLACED  ON  LINE 

+  ',/,'  NUMBER  ONE  OF  THE  OUTPUT  DATA  FILE:  ',///, 

+'  '  2X  ) 

READC*,806)’dATAF 

806  FORMATC A79) 
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835  WRITE(*,831) 

831  FORMATC////, '  DO  YOU  WISH  TO  CHANGE  THE  TITLE  OF  YOUR  OUTPU 

+T  DATA  FILE?' ENTER  Y  FOR  YES  AND  N  FOR  NO:  ' ,2X,  ) 

READ(*,832)  ANS 

832  FORMAT(Al) 

C 

IF(ANS. EQ. 'Y' )  THEN 
GOTO  833 
ELSE 

CONTINUE 

ENTDIF 

IFCANS.EQ.  'N')  THEN 
GOTO  834 
ELSE 

GOTO  835 
834  ENDIF 

C*******5V*************************************************************** 

C**************cOEFFICIENTS  FOR  EPOXY  AND  COPPER  LAYERS***************** 

c 

DO  90  1=1, GPL 
N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

C 

CVf*'***i<r'*****yr*********************QORNERS******************************** 

C 

IF  ((GIH(I).EQ.  l.OR.  GIH(I).EQ.  GDEEP).AND.  (GJH(I).EQ.  1.  OR.  GJH(  I).  EQ 
+.  GWIDE))  THEN 
C 

C  DETERMINE  COEFFICIENTS  FOR  TOP  COPPER 
IF  (GEAT(GIH(I) ,GJH(I)).EQ.  0.  0)  THEN 
GCONd.N)  =  6 

ELSE 

GCON(I,N)  =  7 
ENDIF 
C 

C  CONNECTIONS  FOR  EPOXY  LAYER 

GCON(IB,N)  =  6 
GCON(ID,N)  =  6 
GCON(IC,N)  =  6 
C 
C 

C  LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

C 

C  LEFT  EDGE 

IF  (GJH(I).EQ.  1)  THEN 
C  LEFT  COEFFICIENT 

GOEF(I,N)  =  GYLR21 
GOEF(IC,N)  =  GYLR22 
GOEF(IB,N)  =  GYLR 
GOEF(ID,N)  =  GYLR 
N=N+1 


oooo  nn o  ooo  no  no 


GCON(I,N)  =  7551 
GCON(IB,N)  =  7551 

GCON(IC,N)  =  7551 

GCON(ID,N)  =  7551 

RIGHT  COEFFICIENT 
GOEF(I,N)  =  GYY21 
GOEF(IC,N)  =  GYY22 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10  *  (I+l)  +  1 
GC0N(IB,N)  =  10*  (IB+1)  +  1 
GCONdC.N)  =  10  *  (IC+1)  +  1 

GC0N(ID,N)  =  10  *  (ID+1)  +  1 

RIGHT  EDGE 

ELSEIF  (GJH(I).EQ. GWIDE)  THEN 

LEFT  COEFFICIENT 

GOEFd.N)  =  GYY21 
GOEFdC.N)  =  GYY22 
GOEFCIB.N)  =  GYY 
GOEFCID.N)  =  GYY 


N=N+1 

GCONCI.N)  =  10*d-l)+l 
GCONCIB.N)  =  10*dB-l)+l 
GCONCIC.N)  =  10*dC-l)+l 
GCONCID.N)  =  10*dD-l)+l 


RIGHT  COEFFICIENT 

GOEFCI.N)  =  GYLR21 
GOEFCIC.N)  =  GYLR22 
GOEFCIB.N)  =  GYLR 
GOEFCID.N)  =  GYLR 
N=N+1 

GCONCI.N)  =  7541 
GCONCIB.N)  =  7541 
GCONCIC.N)  =  7541 
GCONCID.N)  =  7541 
END  IF 

FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
FRONT  EDGE 

IF  CGIHCI).EQ. 1)  THEN 
C  FRONT  COEFFICIENT 
G0EFCI,N)=GFXB21 
GOEFCIC.N)=GFXB22 
GOEFCIB.N)=GFXB 
GOEFCID.N)=GFXB 
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N=N+1 

GCONd.N)  =7521 
CtCONCIB.N)  =  7521 
GCON(IC,N)  =  7521 
GCON(ID,N)  =  7521 

C  BACK  COEFFICIENT 
GOEF(I,N)=GXX21 
G0EF(IC,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCONd.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 
GCONCIC.N)  =  10*dC+GWIDE)+l 
GCONCID.N)  =  10*dD+GWIDE)+l 


C  BACK  EDGE 

ELSEIF  (GIH(I). EQ. GDEEP)  THEN 
C  FRONT  COEFFICIENT 

GOEFd.N)=GXX21 
GOEFdC,N)=GXX22 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
N=N+1 

GCONd.N)  =10*d-GWIDE)+l 
GCONCIB.N)  =  10*dB-GWIDE)+l 
GCONCIC.N)  =  10*CIC-GWIDE)+1 
GCONCID.N)  =  10*CID-GWIDE)+1 

C  BACK  COEFFICIENT 

GOEFCI,N)=GFXB21 
GOEFCIC,N)=GFXB22 
GOEFCIB.N)=GFXB 
GOEFCID.N)=GFXB 
N=N+1 

GCONCI.N)  =  7511 
GCONCIB.N)  =  7511 
GCONCIC.N)  =  7511 
GCONCID.N)  =  7511 

ENDIF 

TOP  COEFFICIENT 
GOEFCI.N)=GZE21 
GOEFCIC.N)=GZE22 
GOEFCIB.N)=GZC21 
GOEFCID.N)=GZC22 
N=N+1 

GCONCI.N)  =  7511 
GCONCIC.N)  =  7511 
GCONCIB.N)  =  10*1+1 
GCONCID.N)  =  10*IC+1 
C 
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BOTTOM  COEFFICIENT 

G0EF(I,N)  =  GZC21 
G0EF(IC,N)  =  GZC22 
GOEF(IB,N)  =  GZE22 
G0EF(ID,N)  =  GZB 
N=N+1 

GCONd.N)  =10*(I+GPL)+1 
GCON(IB,N)  =  10*(IB+GPL)+1 
GC0N(IC,N)  =  10*(IC+GPL)+1 
GCON(ID,N)  =  10*(ID+GPL)-fl 


HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0.0)  THEN 
G0EF(I,N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 

GCON(I,N)=9991 

ENDIF 


ELSEIF((GIH(I).EQ.  1.  OR.  GIH(I).  EQ.  GDEEP).  AND.  (GJH(I).NE.  1.  OR.  GJH(I) 
+.NE.  GWIDE))  THEN 

DETERMINE  NUMBER  OF  CONNECTIONS  FOR  COPPER  AND  EPOXY  LAYERS 
IF  (GEAT(GIH(I),GJH(I)).EQ.  0.  0)  THEN 
GC0N(I,N)  =  6 

ELSE 

GC0N(I,N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 
GC0N(ID,N)  =  6 
GC0N(IC,N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS 
LEFT  COEFFICIENT 
G0EF(I,N)  =  GYY21 
G0EF(IC,N)  =  GYY22 
G0EF(IB,N)  =  GYY 
G0EF(ID,N)  =  GYY 
N=N+1 

GC0N(I,N)  =  10*(I-1)+1 
GCON(IB,N)  =  10*(IB-1)+1 
GC0N(IC,N)  =  10*(IC-1)+1 
GC0N(ID,N)  =  10*(ID-1)+1 


RIGHT  COEFFICIENT 
GOEF(I,N)  =  GyY21 
G0EF(IC,N)  =  GYY22 
G0EF(IB,N)  =  GYY 
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GOEF(ID,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10  *  (I+l)  +  1 
GCONCIB.N)  =  10*  (IB+1)  +  1 
GC0N(IC,N)  =  10  *  (IC+1)  +  1 
GCON(ID,N)  =  10  *  (ID+1)  +  1 


FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
FRONT  EDGE 

IF  (GIH(I).EQ.  1)  THEN 
C  FRONT  COEFFICIENT 

GOEF(I,N)=GFXB21 
GOEF(IC,N)=GFXB22 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
N=N+1 

GCON(I,N)  =7521 
GCON(IB,N)  =  7521 
GCON(IC,N)  =  7521 
GCON(ID,N)  =  7521 

C  BACK  COEFFICIENT 

GOEF(I,N)=GXX21 
GOEF(IC,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCON(I,N)  =10*(I+GWIDE)+1 
GCONCIB.N)  =  10*(IB+GWIDE)+1 
GCONCIC.N)  =  10*(IC+GWIDE)+1 
GCONCID.N)  =  10*(ID+GWIDE)+1 

C  BACK  EDGE 

ELSEIF  (GIH(I). EQ. GDEEP)  THEN 
C  FRONT  COEFFICIENT 

GOEF(I,N)=GXX21 
GOEF(IC,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCONCI.N)  =10*(I-GWIDE)+1 
GCONCIB.N)  =  10*CIB-GWIDE)+1 
GCONCIC.N)  =  10*CIC-GWIDE)+1 
GCONCID.N)  =  10*CID-GWIDE)+1 

C  BACK  COEFFICIENT 

GOEFCI,N)=GFXB21 
G0EFCIC.N)=GFXB22 
GOEFCIB.N)=GFXB 
G0EFCID.N)=GFXB 
N=N+1 

GCONCI.N)  =  7511 
GCONCIB.N)  =  7511 
GCONCIC.N)  =  7511 
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GCON(ID,N)  =7511 
ENDIF 
C 

C  TOP  COEFFICIENT 

GOEF(I,N)=GZE21 
GOEF(IC,N)=GZE22 
GOEF(IB,N)=GZC21 
GOEF(ID,N)=GZC22 
N=N+1 

GCONd.N)  =  7511 
GCON(IC,N)  =  7511 
GCON(IB,N)  =  10*1+1 
GCON(ID,N)  =  10*IC+1 
C 

C  BOTTOM  COEFFICIENT 

C 

GOEFd.N)  =  GZC21 
GOEFCIC.N)  =  GZC22 
GOEFCIB.N)  =  GZB 
GOEFCID.N)  =  GZB 
N=N+1 

GCONCI.N)  =10*d+GPL)+l 
GCONCIB.N)  =  10*dB+GPL)  +  l 
GCONCIC.N)  =  10*dC+GPL)+l 
GCONdD.N)  =  10*dD+GPL)+l 
C 

C  HEAT  INPUT 

IF(GEAT(GIHd),GJHd)).NE.  0.  0)  THEN 
GOEFCI.N)  =  GEAT(GIHd),GJHd)) 

N=N+1 

GCONd,N)=9991 

ENDIF 

cxcxcxc 

cxcxcxc 

c 

c 

C  LEFT  AND  RIGHT  EDGES  EXCLUDING  CORNERS 

C 

ELSEIF((GIHd).  EQ.  1.  OR.  GIHd ).  EQ.  GDEEP).AND.  (GJHd).EQ.  1.  OR.  GJHd) 
+. EQ. GWIDE))  THEN 
C 

C  DETERMINE  COEFFICIENTS  FOR  TOP  LAYER 

IF  (GEAT(GIHd)  ,GJHd)).  NE.  0.  0)  THEN 
GCONCI.N)  =  7 

ELSE 

GCONCI.N)  =  6 
ENDIF 
C 

C  CONNECTIONS  FOR  EPOXY  LAYER 

GCONCIB.N)  =  6 
GCONCIC.N)  =  6 
GCONCID.N)  =  6 


C 

c 
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non  a  n  Ci  on  no  no 


C  LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

C 

C  LEFT  EDGE 

IF  (GJH(I).  EQ. 1)  THEN 
C  LEFT  COEFFICIENT 

GOEF(I,N)  =  GYLR21 
GOEF(IC,N)  =  GYLR22 
G0EF(IB,N)  =  GYLR 
GOEF(ID,N)  =  GYLR 
N=N+1 

GCONd.N)  =  7551 
GC0N(IB,N)  =  7551 
GC0N(IC,N)  =  7551 
GC0N(ID,N)  =  7551 

RIGHT  COEFFICIENT 
GOEF(I,N)  =  GYY21 
GOEF(IC,N)  =  GYY22 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
N=N+1 

GC0N(I,N)  =  10  *  (I+l)  +  1 
GC0N(IB,N)  =  10*  (IB+1)  +  1 
GC0N(IC,N)  =  10  *  (IC+1)  +  1 
GC0N(ID,N)  =  10  *  (ID+1)  +  1 

RIGHT  EDGE 

ELSEIF  (GJH(I).EQ. GWIDE)  THEN 

LEFT  COEFFICIENT 
GOEFd.N)  =  GYY21 
GOEFCIC.N)  =  GYY22 
GOEFCIB.N)  =  GYY 
GOEFCID.N)  =  GYY 


N=N+1 

GCONd.N)  =  10*d-l)  +  l 
GCONCIB.N)  =  10*dB-l)  +  l 
GCONCIC.N)  =  10*dC-l)+l 
GCONCID.N)  =  10*dD-l)+l 

RIGHT  COEFFICIENT 

GOEFCI.N)  =  GYLR21 
GOEFCIC.N)  =  GYLR22 
GOEFCIB.N)  =  GYLR 
GOEFCID.N)  =  GYLR 
N=N+1 

GCONCI.N)  =  7541 
GCONCIB.N)  =  7541 
GCONCIC.N)  =  7541 
GCONCID.N)  =  7541 
ENDIF 


FRONT  COEFFICIENT 
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oon  no  non  no 


G0EF(I,N)=GXX21 

GOEF(IC,N)=GXX22 

GOEF(IB,N)=GXX 

GOEF(ID,N)=GXX 

N=N+1 

GCONd.N)  =10*(I-GWIDE)+1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 
GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCON(ID,N)  =  10*(ID-GWIDE)+1 

C  BACK  COEFFICIENT 
G0EF(I,N)=GXX21 
G0EF(IC,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GC0N(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GC0N(IC,N)  =  10*(IC+GWIDE)+1 
GC0N(ID,N)  =  10*(ID+GWIDE)+1 


TOP  COEFFICIENT 
G0EF(I,N)=GZE21 
GOEF(IC,N)=GZE22 
G0EF(IB,N)=GZC21 
G0EF(ID,N)=GZC22 
N=N+1 

GC0N(I,N)  =7511 
GC0N(IC,N)  =  7511 
GC0N(IB,N)  =  10*1+1 
GC0N(ID,N)  =  10*IC+1 

BOTTOM  COEFFICIENT 

G0EF(I,N)  =  GZC21 
G0EF(IC,N)  =  GZC22 
G0EF(IB,N)  =  GZE22 
G0EF(ID,N)  =  GZB 
N=N+1 

GC0N(I,N)  =10*f I+GPL)+1 
GC0N(IB,N)  =  10*(IB+GPL)+1 
GC0N(IC,N)  =  10*(IC+GPL)+1 
GC0N(ID,N)  =  10*(ID+GPL)+1 

HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0.0)  THEN 
GOEFd.N)  =  GEAT(GIHd),GJHd)) 
N=N+1 

GCONd,N)=9991 
END  IF 


C****DETERMINE  COEFFICIENTS  FOR  ALL  NODES  NOT  TOUCHING  AN  EDGE******* 
C 


no 


cjooo  o  onon  on  on 


ELSEIF((GIH(I).NE.  1.  OR.  GIH(  I).  NE.  GDEEP).  AND.  (GJH(I).NE.  1.  OR.  GJH(I) 
+.  NE.GWIDE))  THEN 

DETERMINE  CONNECTIONS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).NE.  0.0)  THEN 
GC0N(I,N)  =  7 

ELSE 

GCONd.N)  =  6 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 
GC0N(IC,N)  =  6 
GC0N(ID,N)  =  6 


LEFT  COEFFICIENT 
G0EF(I,N)  =  GYY21 
G0EF(IC,N)  =  GYY22 

G0EF(IB,N)  =  GYY 

G0EF(ID,N)  =  GYY 

N=N+1 

GC0N(I,N)  =  10  *  (I-l)  +  1 
GC0N(IB,N)  =  10*  (IB-1)  +  1 
GC0N(IC,N)  =  10  *  (IC-1)  +  1 

GC0N(ID,N)  =  10  *  (ID-1)  +  1 


RIGHT  COEFFICIENT 
G0EF(I,N)  =  GYY21 
G0EF(IC,N)  =  GYY22 

G0EF(IB,N)  =  GYY 

G0EF(ID,N)  =  GYY 

N=N+1 

GC0N(I,N)  =  10  *  (I+i)  +  1 
GC0N(IB,N)  =  10*  (IB+1)  +  1 
GC0N(IC,N)  =  10  *  (IC+1)  +  1 

GC0N(ID,N)  =  10  *  (ID+1)  +  1 


FRONT  COEFFICIENT 
G0EF(I,N)=GXX21 
GOEF(IC,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GC0N(I,N)  =10*(I-GWIDE)+1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 
GC0N(IC,N)  =  10*(IC-GWIDE)+1 
GC0N(ID,N)  =  10*(ID-GWIDE)+1 


C  BACK  COEFFICIENT 
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nnovo  on  non  on 


G0EF(I,N)=GXX21 

G0EF(IC,N)=GXX22 

GOEF(IB,N)=GXX 

GOEF(ID,N)=GXX 

N=N+1 

GCON(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GCON(IC,N)  =  10*( IC+GWIDE)+1 
GCCN(ID.N)  =  10*(ID+GWIDE)+1 


TOP  COEFFICIENT 
G0EF(I.N)=GZE21 
G0EF(IC,N)=GZE22 
G0EF(IB,N)=GZC21 
GOEF(ID,N)=GZC22 
N=N+1 

GCONd.N)  =  7511 
GC0N(IB,N)  =  10*1+1 
GC0N(IC,N)  =  7511 
GC0N(ID,N)  =  10*IC+1 


BOTTOM  COEFFICIENT 

GOEFd.N)  =  GZC21 
GOEFdC.N)  =  GZC22 
GOEFCIB.N)  =  GZB 
GOEFCID.N)  =  GZE22 
N=N+1 

GCONd.N)  =10*d+GPL)+l 
GCONCIB.N)  =  10*dB+GPL)+l 
GCONCIC.N)  =  10*dC+GPL)+l 
GCONCID.N)  =  10*dD+GPL)+l 


HEAT  INPUT 

IF(GEAT(GIHd),GJHd)).NE.  0.  0)  THEN 
GOEFd.N)  =  GEAT(GIHd).GJHd)) 

N=N+1 

GC0Nd.N)=9991 
END  IF 
END  IF 
CONTINUE 

GENERATE  DATA  FILE 

OPEN  ( 3. F I LE=NAME.FORM=' FORMATTED' ,ACCESS=' DIRECT' .RECL=108 . STATUS 
+='NEW' ) 

U'RITE(3.909)  DATAF 
909  FORMATdX.A79) 

WR ITE ( 3 . 9 08 )  COUNT . CONTEMP , ZER . ZER . ZER . ZER . ZER . ZER , GSEL 
908  F0RMAT(2X.9d3,5X)) 

WRITE(3,907)  ZER, ZER. ZER 
907  FORMAT(2X,3d3,5X)) 


112 


WRITEC 3,9081)  NMAX ,TMAX,HTRS ,D1 ,D2 ,D3 ,D4,D5 ,D6 ,D7 
9081  F0RMAT(2X,9(I3,5X)) 

WR ITE ( 3 , 905 )  ACC , DAMP , MAXIT , CONFAC , GI BT 

905  F0RMAT(1X,2(F9.  7 , IX) , 14X, 12 , 1X,F9.  7,1X,F9.5) 

WRITEC 3,906)  GPRT,GFT,GBT,GRT,GLT,GWRT 

906  FORMATC 1X,6(F12. 3,1X)) 

DO  112  I=l,LOCVAR 

WRITEC  3,9100)  GCONC 1 , 1) ,GCONC 1 ,2) ,GCONC 1 , 3) ,GCONC 1 ,4) ,GCONC 1 ,5) ,GC 
+0NCI,6),GC0NCI,7),GC0NCI,8) 

9100  FORMATC 14, 3X,7C 15, 7X)) 

WRITEC 3 ,9110)  GOEFC 1 , 1) ,GOEFC 1,2) ,GOEFC 1,3) ,GOEFC 1,4) ,GOEFC 1 ,5) ,GO 
+EFCI,6),GOEFCI,7) 

9110  FORMATC 7CF9.  3, 3X)) 

112  CONTINUE 
CLOSE  C3) 

CALL  CLS 

WRITEC *,999)  NAME 

999  FORMATC////,'  THE  OUTPUT  DATA  HAS  BEEN  PLACED  IN  A  FILE 

+NAMED  ',A6, /////,'  <PRESS  ENTER  TO  CONTINUE>') 

kEADC",5912)  ANS 
5912  FORMATC Al) 

ENT) 

C 

C 

C  TITLE:  MODEL  BUILDER 

C  SUBROUTINE:  S3 

C  DATE:  09  JUL  91 

C  AUTHOR:  LT  STEVE  GLASER 

C  COMPILER:  MICROSOFT  VERSION  4.01 

C  LINKER:  MICROSOFT  VERSION  3.55 

C 

C  SUBROUTINE  S3  IS  CALLED  FROM  SUBROUTINE  PCBS3.  S3  GENERATES  THE 
C  OUTPUT  DATA  FILE  FOR  THE  THREE  COPPER  LAYER  CASE. 

C 

SUBROUTINE  S3  CGYLR ,GFXB ,GYY,GXX,GZC31 ,GZC32 ,GZC33 ,GZB ,GYLR31 ,GYLR 
+32,GYLR33,GFX331,GFXB32,GFXB33,GYY31,GYY32,GYY33,GXX31,GXX32,GXX33 
+,GZE31,GZE32,GZE33,GPL,GIH,GJH,GEAT,GWIDE,GDEEP,GOEF,GIBT,GPRT,GWR 
+T,GRT,GLT,GFT,GBT,GSEL) 

REAL  GYLR , GFXB , GYY , GXX , GZC3 1 , GZB , GYLR3 1 , GFXB3 1 , GYY3 1 , GXX3 1 , GZE3 1 , A 
+CC , DAMP , CONFAC , GI BT , GPRT , GWRT , GRT , GLT, GFT , GBT , GYLR32 , GFXB32 , GYY32 , 
+GXX32 ,GZE32 ,GZC32 ,GYLR33 ,GFXB33 ,GYY33,GXX33 ,GZE33 ,GZC33 

REAL  GEATC 100, 100) , GOEFC 1000,12) 


CHARACTER*!  SELECT, ANS 
CHARACTER  DATAF*79 ,NAME*6 

INTEGER  GPL , I , N , I B , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS,D1,D2,D3,D4,D5,D6,D7,MAXIT,L0CVAR,IC,ID,IE,IF 

INTEGER  GIHC 1000) ,GJHC 1000) ,GCON( 1000, 12) 

CONTEMP=6 

ZER=0 
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NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0.  8 
COUNT=6*GPL 
LOCVAR=COUNT 


C 

C 

898  CALL  CLS 


803 


804 

552 

910 

911 
811 


812 

C 


897 


WRITE(*,803) 
FORMATC///, ' 

+  INTO  THE’ ,/, 


THIS  PROGRAM  CREATES  AN  OUTPUT  DATA  FILE  FOR  ENTRY 


+  ' 
+  ' 
+  ’ 
+  ' 

+  ’ 


EXISTING  THERMAL  ANALYZER,  FURTHERMORE,  THIS  PROGRAM  DOES',/, 
NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE.  THEREFORE',/, 
THE  USER  WILL  NAME  THE  DATA  FILE  FOR  EACH 
PROGRAM.  THE  FILE  NAME  IS  LIMITED  TO  SIX 
SHOULD  NOT  HAVE  ANY  SPACES. 

PLEASE  ENTER  THE  DESIRED  DATA  FILE  NAME: 

READ(*,804)  NAME 
FORMAT(A6) 

WRITE(*,910)NAME 

FORMATC/////,'  YOU  SELECTED  ',A6,’  FOR  YOUR  DATA  FILE  nAME 

+' .//) 


RUN  OF  THIS  ' ,/, 
CHARACTERS,  A.ND’,/, 
'  ,///, 

’.2X  ) 


WRITEC*  811) 

FORMATC'  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N 
+FOR  NO:  ' ,2X,  ) 

READC*,812)  ANS 
FORMATC Al) 


IFCANS. EQ. 'N' )  THEN 
GOTO  898 
ELSE 

CONTINUE 

ENDIF 

IFCANS. EQ.  'Y' )  THEN 
GOTO  897 
ELSE 

CALL  CLS 
GOTO  552 
ENDIF 


C  ALLOW  THE  USER  TO  PROVIDE  TITLE  LINE  FOR  DATA  FILE 


114 


833  CALL  CLS 
WRITE(*,805) 

805  FORMATC////, ’  ENTER  THE  DESIRED  TITLE  TO  BE  PLACED  ON  LINE 

+  NUMBER  ONE  OF  THE  OUTPUT  DATA  FILE:  ',///, 

+'  ' , 2X,  ) 

READ(*,806)’DAfAF 

806  FORMATC A79) 

835  WRITE(*,831) 

831  FORMATC////,'  DO  YOU  WISH  TO  CHANGE  THE  TITLE  OF  YOUR  OUTPU 

+T  DATA  FILE?',/,'  ENTER  Y  FOR  YES  AND  N  FOR  NO:  ',2X,  ) 

READC*,832)  ANS 

832  FORMATC Al) 

C 

IFCANS.EQ. 'Y')  THEN 
GOTO  833 
ELSE 

CONTINUE 

ENDIF 

IFCANS. EQ. 'N' )  THEN 
GOTO  834 
ELSE 

GOTO  835 

834  ENDIF 
C 

c*********************************************************************** 
C**************coEFFICIENTS  FOR  EPOXY  AND  COPPER  LAYERS***************** 

C********Vc*****************************Vr******************************** 

c 

DO  90  1=1, GPL 
N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

IE=4*GPL+I 

IF=5*GPL+I 


(]*********************************QQjy^£J^g******************************** 

c 

IF  CCGIHCD.EQ.  l.OR.  GIHCI).EQ.GDEEP).AND.  CGJHCI).EQ.  1.  OR.  GJHCI).  EQ 
+. GWIDE))  THEN 
C 

C  DETERMINE  COEFFICIENTS  FOR  TOP  COPPER 

IF  CGEATCGIHCI), GJHCI)). EQ. 0.0)  THEN 
GCONCI.N)  =  6 

ELSE 

GCONCI,N)  =  7 
ENDIF 
C 

C  CONNECTIONS  FOR  EPOXY  LAYER 

GCONCIB,N)  =  6 
GCONCID,N)  =  6 
GCONCIC,N)  =  6 
GCONCID,N)  =  6 
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GCON(IE,N)  =  6 
GC0N(IF,N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
LEFT  EDGE 

IF  (GJH(I).EQ. 1)  THEN 
LEFT  COEFFICIENT 
GOEF(I,N)  =  GYLR31 
G0EF(IC,N)  =  GYLR32 
G0EF(IE,N)  =  GYLR33 
GOEF(IB,N)  =  GYLR 
GOEF(ID,N)  =  GYLR 
GOEF(IF,N)  =  GYLR 
N=N+1 

GCON(I,N)  =  7551 
GCON(IB,N)  =  7551 
GCON(IC,N)  =  7551 
GC0N(ID,N)  =  7551 
GCON(IE,N)  =  7551 
GCON(IF,N)  =  7551 


RIGHT  COEFFICIENT 
GOEF(I,N)  =  GYY31 
G0EF(IC,N)  =  GYY32 
G0EF(IE,N)  =  GYY33 
GOEF(IB,N)  =  GYY 
G0EF(ID,N)  =  GYY 
G0EF(IF,N)  =  GYY 


N=N+1 


GCON(I,N)  =  10  *  (I+l)  +  1 
GC0N(IB,N)  =  10*  (IB+l)  +  1 
GC0N(IC,N)  =  10  *  (IC+1)  +  1 

GC0N(ID,N)  =  10  *  (ID+1)  +  1 

GC0N(IE,N)  =  10  *  (IE+1)  +  1 

GC0N(IF,N)  =  10  *  (IF+1)  +  1 


RIGHT  EDGE 

ELSEIF  (GJH(I). EQ. GWIDE)  THEN  % 

LEFT  COEFFICIENT 

G0EF(I,N)  =  GYY31 
G0EF(IC,N)  =  GYY32 
GOEF(IE,N)  =  GYY33 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
G0EF(IF,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10*(I-1)+1 
GC0N(IB,N)  =  10*(IB-1)+1 
GC0N(IC,N)  =  10*(IC-1)+1 
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noon  o  o  o 


GCON(ID,N)  =  10*(ID-1)+1 
GCON(IE,N)  =  10*(IE-1)+1 
GCON(IF,N)  =  10*(IF-1)+1 

RIGHT  COEFFICIENT 

GOEF(I,N)  =  GYLR31 
GOEF(IC,N)  =  GYLR32 
G0EF(IE,N)  =  GYLR33 
G0EF(IB,N)  =  GYLR 
G0EF(ID,N)  =  GYLR 
G0EF(IF,N)  =  GYLR 
N=N+1 

GCON(I,N)  =  7541 
GCON(IB,N)  =  7541 
GC0N(IC,N)  =  7541 
GC0N(ID,N)  =  7541 
GCON(IE,N)  =  7541 
GC0N(IF,N)  =  7541 

ENDIF 

FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
FRONT  EDGE 

IF  (GIH(I).EQ.  1)  THEN 
C  FRONT  COEFFICIENT 

GOEF(I,N)=GFXB31 
G0EF(IC,N)=GFXB32 
G0EF(IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(EF,N)=GFXB 
N=N+1 

GCONd.N)  =  7521 
GCON(IB,N)  =  7521 
GCON(IC,N)  =  7521 
GCON(ID,N)  =  7521 
GCON(IE,N)  =  7521 
GCON(IF,N)  =  7521 
C  BACK  COEFFICIENT 

GOEF(I,N)=GXX31 
GOEF(IC,N)=GXX32 
G0EF(IE,N)=GXX33 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
N=N+1 

GCON(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GCON(IC,N)  =  10*(IC+GWIDE)+1 
GCON(ID,N)  =  10*(ID+GWIDE)+1 
GCON(IE,N)  =  10*(IE+GWIDE)+1 
GCON(IF,N)  =  10*(IF+GWIDE)+1 
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o  n  n  o  o 


C 


BACK  EDGE 

ELSEIF  (GIH(I).EQ.  GDEEP)  THEN 
C  FRONT  COEFFICIENT 

GOEF(I,N)=GXX31 
GOEF(IC,N)=GXX32 
G0EF(IE,N)=GXX33 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
N=N+1 

GCON(I,N)  =  10*(I-GWIDE)+1 
GCON(I3,N)  =  10*(IB-GWIDE)+1 
GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCON(ID,N)  =  10*(ID-GWIDE)+1 
GCON(IE,N)  =  10*(IE-GWIDE)+1 
GCON(IF,N)  =  10*(IF-GWIDE)+1 

C  BACK  COEFFICIENT 

GOEF(I,N)=GFXB31 
GOEF(IC,N)=GFXB32 
G0EF(IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF,N)=GFXB 
N=N+1 

GCONd.N)  =  7511 
GCON(IB,N)  =  7511 
GCONdC.N)  =  7511 
GCONCID.N)  =  7511 
GCONCIE.N)  =  7511 
GCONCIF.N)  =  7511 

ENDIF 

TOP  COEFFICIENT 
GOEFd,N)=GZE31 
GOEFdC,N)=GZE32 
GOEFdE,N)=GZE33 
GOEFdB,N)=GZC31 
GOEFdD,N)=GZC32 
GOEFdF,N)=GZC33 
N=N+1 

GCONCI.N)  =  7511 
GCONCIC.N)  =  7511 
GCONCIE.N)  =  7511 
GCONCIB.N)  =  10*1+1 
GCONCID.N)  =  10*IC+1 
GCONdF.N)  =  10*IE+1 

BOTTOM  COEFFICIENT 

GOEFCI.N)  =  GZE31 
GOEKIC.N)  =  GZE32 
GOEFCIE.N)  =  GZE33 
GOEFCIB.N)  =  GZC31 
GOEFdD.N)  =  GZC32 
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no  o n o o  on  no 


G0EF(IF,N)  =  GZB 
N=N+1 

GC0N(I,N)  =10*(I+GPL)+1 
GCON(IB,N)  =  10*(IB+GPL)+1 
GC0N(IC,N)  =  10*(IC+GPL)+1 
GC0N(ID,N)  =  10*(ID+GPL)+1 


C 

C  HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE. 0. 0)  THEN 
GOEFd.N)  =  GEAT(GIH(I),GJH(I)) 

N=N+1 

GCON(I,N)=9991 

ENDIF 

C 

ELSEIF((GIH(I).EQ.  1.  OR.  GIH( I). EQ.  GDEEP) . AND. (GJH(I).NE.  1.  OR.  GJH(I) 
+.NE.GWIDE))  THEN 

DETERMINE  NUMBER  OF  CONNECTIONS  FOR  COPPER  AND  EPOXY  LAYERS 
IF  (GEAT(GIH(I) ,GJH(I)).EQ.  0. 0)  THEN 
GCONd.N)  =  6 

ELSE 

GCONd.N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GCONCIB.N)  =  6 
GCONCID.N)  =  6 
GCONCIC.N)  =  6 
GCONdE.N)  =  6 
GCONCIF.N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS 

LEFT  COEFFICIENT 
GOEFd.N)  =  GYY31 
GOEFCIC.N)  =  GYY32 
GOEFCIE.N)  =  GYY33 
GOEFCIB.N)  =  GYY 
GOEFCID.N)  =  GYY 
GOEFCIF.N)  =  GYY 
N=N+1 

GCONd.N)  =  10*d-l)+l 
GCONCIB.N)  =  10*dB-l)+l 
GCONCIC.N)  =  10*CIC-1)+1 
GCONCID.N)  =  10*CID-1)+1 
GCONCIE.N)  =  10*CIE-1)+1 
GCONCIF.N)  =  10*CIF-1)+1 

RIGHT  COEFFICIENT 
GOEFCI.N)  =  GYY31 
GOEFCIC.N)  =  GYY32 
GOEFCIE.N)  =  GYY33 
GOEFCIB.N)  =  GYY 
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GOEF(ID,N)  =  GYY 

GOEF(IF,N)  =  GYY 

N=N+1 

GCONd.N)  =  10  *  (I+l)  +  1 
GCON(IB,N)  =  10*  (IB+1)  +  1 
GCON(IC,N)  =  10  *  (IC+1)  +  i 

GCON(ID,N)  =  10  *  (ID+1)  +  1 

GCON(IE,N)  =  10  *  (IE+1)  +  1 

GCON(IF,N)  =  10  *  (IF+1)  +  1 

FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

FRONT  EDGE 

IF  (GIH(I).EQ. 1)  THEN 

C  FRONT  COEFFICIENT  * 

G0EF(I,N)=GFXB31 
G0EF(IC,N)=GFXB32 
G0EF(IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF,N)=GFXB 
N=N+1 

GCONd.N)  =7521 
GCONdB.N)  =  7521 
GCONCIC.N)  =  7521 
GCONdD.N)  =  7521 
GCONdE.N)  =  7521 
GCONCIF.N)  =  7521 

C  BACK  COEFFICIENT 

G0EFd,N)=GXX31 
G0EFdC,N)=GXX32 
GOEFdE,N)=GXX33 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF,N)=GXX 
N=N+1 

GCONd.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 
GCONCIC.N)  =  10*dC+GWIDE)+l 
GCONCID.N)  =  10*dD+GWIDE)+l 
GCONCIE.N)  =  10*dE+GWIDE)+l 

GCONCIF.N)  =  10*CIF+GWIDE)+1  1 

C  BACK  EDGE 

ELSEIF  CGIHCD.EQ.  GDEEP)  THEN 
C  FRONT  COEFFICIENT 

G0EFCI,N)=GXX31 
G0EFCIC,N)=GXX32 
G0EFCIE,N)=GXX33 
GOEFCIB,N)=GXX 
GOEFCID,N)=GXX 
GOEFCIF.N)=GXX 
N=N+1 

GCONCI.N)  =10*CI-GWIDE)+1 
GCONCIB.N)  =  10*CIB-GWIDE)+1 
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GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCONCID.N)  =  10*(ID-GWIDE)+1 
GCON’(IE,N)  =  10*(IE-GWIDE)+1 
GCON(IF,N)  =  10*(IF-GWIDE)+1 

C  BACK  COEFFICIENT 
GOEF(I,N)=GFXB31 
GOEF(IC,N)=GFXB32 
GOEF(IE,N)=GFXB33 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF,N)=GFXB 
N=N+1 

GCONd.N)  =  7511 
GCON(IB,N)  =7511 
GCON(IC,N)  =  7511 
GCON(ID,N)  =  7511 
GCON(IE,N)  =7511 
GCON(IF,N)  =7511 
END  IF 
C 

C  TOP  COEFFICIENT 
GOEF(I.N)=GZE31 
GOEF(IC,N)=GZE32 
GOEF(IE,N)=GZE33 
GOEF(IB,N)=GZC31 
GOEF(ID,N)=GZC32 
G0EF(IF,N)=GZC33 

N=N+1 

GCONd.N)  =  7511 
GCONCIC.N)  =  7511 
GCONCIB.N)  =  10*1+1 
GCONCID.N)  =  10*IC+1 
GCONCIE.N)  =7511 
GCONCIF.N)  =  10*IE+1 


C 

C  BOTTOM  COEFFICIENT 
C 

GOEFCI.N)  =  GZE31 
GOEFCIC.N)  =  GZE32 
GOEFCIE.N)  =  GZE33 
GOEFCIB.N)  =  GZC31 
GOEFCID.N)  =  GZC32 
GOEFCIF.N)  =  GZB 


N=N+1 


GCONd.N)  =10*d+GPL)+l 
GCONCIB.N)  =  10*CIB+GPL)+1 
GCONCIC.N)  =  10*CIC+GPL)+1 
GCONCID.N)  =  10*CID+GPL)+1 
GCONCIE.N)  =  10*CIE+GPL)+1 
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GCON(IF,N)  =  10*(IF+GPL)+1 


C  HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE. 0.0)  THEN 
GOEFd.N)  =  GEAT(GIH(I)  ,GJH(I)) 
N=N+1 

GCON(I,N)=9991 
END  IF 

cxcxcxc 

cxcxcxc 


LEFT  AND  RIGHT  EDGES  EXCLUDING  CORNERS 

ELSEIFC ( GIH( I ) .  EQ.  1. OR.  GIH( I ) . EQ. GDEEP) . AND. ( GJH( I) .  EQ. 
+. EQ. GWIDE))  THEN 

DETERMINE  COEFFICIENTS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).NE.O.O)  THEN 
GCON(I,N)  =  7 

ELSE 

GCON(I,N)  =  6 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 
GCONdC.N)  =  6 
GCONCID.N)  =  6 
GCONCIE.N)  =  6 
GCONCIF.N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
LEFT  EDGE 

IF  (GJHd).  EQ.  1)  THEN 
LEFT  COEFFICIENT 
GOEFCI.N)  =  GYLR31 
GOEFdC.N)  =  GYLR32 
GOEFCIE.N)  =  GYLR33 
GOEFCIB.N)  =  GYLR 
GOEFdD.N)  =  GYLR 
GOEFCIF.N)  =  GYLR 
N=N+1 

GCONCI.N)  =  7551 
GCONCIB.N)  =  7551 
GCONdC.N)  =  7551 
GCONCID.N)  =  7551 
GCONCIE.N)  =  7551 
GCONCIF.N)  =  7551 

RIGHT  COEFFICIENT 
GOEFCI.N)  =  GYY31 
GOEFCIC.N)  =  GYY32 
GOEFCIE.N)  =  GYY33 


.OR.  GJHCI) 
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c 

c 

c 

c 


c 

c 

c 


c 

c 


GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
N=N+1 


GCON(I,N)  =  10  *  (I+l)  +  1 
GCON(IB,N)  =  10*  (IB+1)  +  1 
GCON(IC,N)  =  10  *  (IC+1)  +  1 

GCON(ID,N)  =  10  *  (ID+1)  +  1 

GC0N(IE,N)  =  10  *  (IE+1)  +  1 

GCON(IF,N)  =  10  *  (IF+1)  +  1 


RIGHT  EDGE 

ELSEIF  (GJH(I).EQ. GWIDE)  THEN 


LEFT  COEFFICIENT 
GOEFd.N)  =  GYY31 
GOEF(IC,N)  =  GYY32 
GOEF(IE,N)  =  GYY33 
GOEF(IB.N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10*(I-1)+1 
GCON(IB,N)  =  10*(IB-1)+1 
GCON(IC,N)  =  10*(IC-1)+1 
GCON(ID,N)  =  10*(ID-1)+1 
GCON(IE,N)  =  10*(IE-1)+1 
GCON(IF,N)  =  10*(IF-1)+1 


RIGHT  COEFFICIENT 


GOEF(I,N) 

GOEF(IC,N) 

GOEF(IE,N) 

GOEF(IB,N) 

GOEF(ID,N) 

GOEF(IF,N) 

N=N+1 

GCONd.N) 

GCONCIB.N) 

GCONCIC.N) 

GCONCID.N) 

GCONCIE.N) 

GCONCIF.N) 

ENDIF 


GYLR31 

GYLR32 

GYLR33 

GYLR 

GYLR 

GYLR 

7541 

7541 

7541 

7541 

7541 

7541 


FRONT  COEFFICIENT 
GOEFd.N)  =GXX31 
GOEFdC,N)=GXX32 
GOEFdE,N)=GXX33 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF,N)=GXX 
N=N+1 

GCONd.N)  =10*d-GWIDE)+l 
GCONCIB.N)  =  10*dB-GWIDE)+l 
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GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCON(ID,N)  =  10*(ID-GWIDE)+1 
GCON(IE,N)  =  10*(IE-GWIDE)+1 
GCON(IF,N)  =  10*(IF-GWIDE)+1 
C  BACK  COEFFICIENT 
G0EF(I,N)=GXX31 
GOEF(IC,N)=GXX32 
GOEF(IE,N)=GXX33 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
N=N+1 

GCONd.N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GCON(IC,N)  =  10*(IC+GWIDE)+1 
GCON(ID,N)  =  10*(ID+GWIDE)+1 
GCON(IE,N)  =  10*(IE+GWIDE)+1 
GCON(IF,N)  =  10*(IF+GWIDE)+1 

TOP  COEFFICIENT 
GOEF(I,N)=GZE31 
GOEF(IC,N)=GZE32 
GOEF(ID,N)=GZE33 
GOEF(IB,N)=GZC31 
GOEF(ID,N)=GZC32 
GOEF(IF,N)=GZC33 
N=N+1 

GCONCI.N)  =  7511 
GCON(IC,N)  =  7511 
GCON(IB,N)  =  10*1+1 
GCONCID.N)  =  10*IC+1 
GCONdE.N)  =  7511 
GCONCIF.N)  =  10*IE+1 

BOTTOM  COEFFICIENT 


GOEFCI.N)  =  GZE31 
GOEFCIC.N)  =  GZE32 
GOEFCIE.N)  =  GZE33 


GOEFCIB.N) 

GOEFdD,N) 

GOEFCIF.N) 

N=N+1 

GCONCI.N) 

GCONCIB.N) 

GCONCIC.N) 

GCONCID.N) 

GCONCIE.N) 

GCONCIF.N) 


GZC31 

GZC32 

EZB 

10*CI+GPL)+1 

10*CIB+GPL)+1 

10*CIC+GPL)+1 

10*CID+GPL)+1 

10*CIE+GPL)+1 

10*CIF+GPL)+1 


« 


HEAT  INPUT 

IFCGEATCGIHCI),GJHCI)).NE. 0.0)  THEN 
GOEFCI.N)  =  GEATCGIHCD.GJHCD) 
N=N+1 
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GC0N(I,N)=9991 

ENDIF 


****DETERMINE  COEFFICIENTS  FOR  ALL  NODES  NOT  TOUCHING  AN  EDGE******* 


ELSEIF((GIH(I).  NE.  1.  OR.  GIH(I).  NE.  GDEEP). AND.  (GJH(I).NE.  1.  OR.  GJH( I) 
+. NE.GWIDE))  THEN 

DETERMINE  CONNECTIONS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).NE.O. 0)  THEN 
GCON(I,N)  =  7 

ELSE 

GCON(I,N)  =  6 
ENDIF 


CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 
GCON(IC,N)  =  6 
GCON(ID,N)  =  6 
GC0N(IE,N)  =  6 
GCON(IF,N)  =  6 


LEFT  COEFFir.r 
G0EF(I,N)  viYYSl 
G0EF(IC..>,)  =  GYY32 
GOEFCL.N)  =  GYY33 
GOEFCIB.N)  =  GYY 
GOLF(ID,N)  =  GYY 
G0EF(IF,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10  *  (I-l)  +  1 
GC0N(IB,N)  =  10*  (IB-1)  +  1 
GC0N(IC,N)  =  10  *  (IC-1)  +  1 

GCON(ID,N)  =  10  *  (ID-1)  +  1 

GCON(IE,N)  =  10  *  (IE-1)  +  1 

GC0N(IF,N)  =  10  *  (IF-1)  +  1 

C  RIGHT  COEFFICIENT 
GOEF(I,N)  =  GYY31 
GOEF(IC,N)  =  GYY32 
G0EF(IE,N)  =  GYY33 
G0EF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
N=N+1 

GC0N(I,N)  =  10  *  (I+l)  +  1 
GCON(IB,N)  =  10*  (IB+1)  +  1 
GCONdC.N)  =  10  *  (IC+1)  +  1 

GC0N(ID,N)  =  10  *  (ID+1)  +  1 

GCON(IE,N)  =  10  *  (IE+1)  +  1 
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GCON(IF,N)  =  10  *  (IF+1)  +  1 


FRONT  COEFFICIENT 
G0EF(I,N)=GXX31 
GOEF(IC,N)=GXX32 
G0EF(IE,N)=GXX33 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
N=N+l 

GCONCI.N)  =10*(I-GWIDE)+1 
GC0N(IB,N)  =  10*(IB-GWIDE)+1 
GCONdC.N)  =  10*(IC-GWIDE)+1 
GCONdD.N)  =  10*dD-GWIDE)+l 
GCONCIE.N)  =  10*dE-GWIDE)  +  l 
GCONCIF.N)  =  10*dF-GWIDE)+l 
BACK  COEFFICIENT 
GOEFd,N)=GXX31 
GOEFdC,N)=GXX32 
GOEFdE,N)=GXX33 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF.N)=GXX 
N=N+1 

GCONCI.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 
GCONCIC.N)  =  10*dC+GWIDE)+l 
GCONCID.N)  =  10*dD+GWIDE)  +  l 
GCONCIE.N)  =  10*CIE+GWIDE)+1 
GCONCIF.N)  =  10*CIF+GWIDE)+1 
TOP  COEFFICIENT 
GOEFCI,N)=GZE31 
GOEFCIC.N)=GZE32 
GOEFCIE.N)=GZE33 
GOEFCIB,N)=GZC31 
GOEFCID,N)=GZC32 
GOEFCIF.N)=GZC33 
N=N+1 

GCONCI.N)  =  7511 
GCONCIB.N)  =  10*1+1 
GCONCIC.N)  =  7511 
GCONCID.N)  =  10*IC+1 
GCONCIE.N)  =  7511 
GCONCIF.N)  =  10*IE+1 


BOTTOM  COEFFICIENT 

GOEFCI.N)  =  GZE31 
GOEFCIC.N)  =  GZE32 
GOEFCIE.N)  =  G2E33 
GOEFCIB.N)  =  GZC31 
GOEFCID.N)  =  GZC32 
GOEFCIF.N)  =  GZC33 
N=N+1 

GCONCI.N)  =10*CI+GPL)+1 
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GCON(IB,N)  =  10*(IB+GPL)+1 
GCON(IC,N)  =  10*(IC+GPL)+1 
GCON(ID,N)  =  10*(ID+GPL)+1 
GCON(IE,N)  =  10*(IE+GPL)+1 
GCONCIF,N)  =  10*(IF+GPL)+1 
C 

C  HEAT  INPUT 

IF(GEAT(GIH(I) ,GJH(I)).NE. 0.  0)  THEN 
GOEFd.N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 


GCON(I,N)=9991 

ENDIF 

END  IF 

90  CONTINUE 

C 

C  GENERATE  DATA  FILE 

C 

OPEN  ( 3, FILE=NAME,F0RM=' FORMATTED' ,ACCESS=’ DIRECT* ,RECL=108, STATUS 
+='NEW' ) 

WRITE(3,909)  DATAF 
909  F0RMAT(1X,A79) 

WRITE( 3 , 908 )  COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908  F0RMAT(2X,9(I3,5X)) 

WRITE(3,907)  ZER, ZER, ZER 
907  F0RMAT(2X,3(I3,5X)) 

WRITE(3,9081)  NMAX,TMAX,HTRS,D1,D2,D3,D4,D5,D6,D7 
9081  FORMAT(2X,9(I3,5X)) 

WRITE(3,905)  ACC , DAMP ,MAXIT,CONFAC ,GIBT 

905  F0RMAT(1X,2(F9. 7,1X),14X,I2,1X,F9. 7,1X,F9.5) 

WR I TE  (  3 , 9  0  6  )  GPRT ,  GFT ,  GBT ,  GRT ,  GLT ,  GV,TIT 

906  F0RMAT(1X,6(F12. 3,1X)) 

DO  112  I=1,L0CVAR 

WRITE( 3,9100)  GC0N(I,1),GC0N(I,2),GC0N(I,3),GC0N(I,4),GC0N(I,5),GC 
+0N(I,6) ,GCON(I,7),GCON(I,8) 

9100  FORMAT(I4,3X,7(I5,7X)) 

WRITE ( 3 , 9 1 1 0 )  GOEF( 1,1), GOEF( 1,2), GOEF( 1,3), GOEF( 1,4), GOEF( I , 5 ) , GO 
+EF(I,6),GOEF(I,7) 

9110  FORMAT(7(F9. 3,3X)) 

112  CONTINUE 
CLOSE  (3) 

CALL  CLS 

WRITE(*,999)  NAME 

999  FORMATC ////,'  THE  OUTPUT  DATA  HAS  BEEN  PLACED  IN  A  FILE 

+NAMED  ',A6, /////,'  <PRESS  ENTER  TO  CONTINUE>’) 

READ(*,5912)  ANS 
5912  FORMATC Al) 

END 

C********************************************** ************************ 


C 

c 

c 

c 

c 

c 

c 

c 


TITLE:  MODEL  BUILDER 

SUBROUTINE:  S4 

DATE:  09  JUL  91 

AUTHOR:  LT  STEVE  GLASER 

COMPILER:  MICROSOFT  VERSION  4.  01 

LINKER:  MICROSOFT  VERSION  3.  55 

SUBROUTINE  S4  IS  CALLED  FROM  SUBROUTINE  PCBS3. 


S4  GENERATES  THE 
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OUTPUT  DATA  FILE  FOR  THE  FOUR  COPPER  LAYER  CASE. 


SUBROUTINE  S4  (GYLR ,GFXB ,GYY,GXX,GZC41 ,GZC42 ,GZC43 ,GZC44 ,GZB .GYLR4 
+1 , GYLR42 , GYLR43 , GYLR44 , GFXB4 1 , GFXB42 , GFXB43 , GFXB44 , GYY4 1 , GYY42 , GYY 
+43 , GYY44 , GXX4 1 , GXX42 , GXX43 , GXX44 , GZE4 1 , GZE42 , GZE43 , GZE44 , GPL , GIH , G 
+ JH , GE AT , GW I DE , GDEE  P , GOEF , G I BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GSEL ) 

REAL  GYLR , GFXB , GYY , GXX , GZC4 1 , GZB , GYLR4 1 , GFXB4 1 , GYY4 1 , GXX4 1 , GZE4 1 , A 
+CC , DAMP , CONFAC , G I BT , GPRT , GTOT , GRT , GLT , GFT , GBT , GYLR42 , GFXB42 , GYY42 , 
+GXX42 , GZE42 , GZC42 , GYLR43 , GFXB43 , GYY43 , GXX43 , GZE43 , GZC43 , GYLR44 , GFX 
+B44 , GYY44 , GXX44 , GZE44 , GZC44 

REAL  GEAT(100,100) ,GOEF(1000,12) 


CHARACTER*!  SELECT, ANS 
CHARACTER  DATAF*79 ,NAME*6 

INTEGER  GPL , I , N , I B , GW I DE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS,Dl,D2,D3,D4,D5,D6,D7,MAXIT,LOCVAR,IC,ID,IE,IF,IG,IH 

INTEGER  GIH( 1000) ,GJH( 1000) ,GCON( 1000,12) 

CONTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0.  666667 
MAXIT=12 
CONFAC=0.  8 
COUNT=8*GPL 
LOCVAR=COUNT 


898  CALL  CLS 


WRITE(*,803) 

803  FORMATC///,'  THIS  PROGRAM  CREATES  AN  OUTPUT  DATA  FILE  FOR  ENTRY 
+  INTO  THE ' , / , 

+'  EXISTING  THERMAL  ANALYZER,  FURTHERMORE,  THIS  PROGRAM  DOES',/, 

+'  NOT  ERASE  OR  WRITE  OVER  THE  EXISTING  DATA  FILE.  THEREFORE',/, 

+  '  THE  USER  WILL  NAME  THE  DATA  FILE  FOR  EACH  RUN  OF  THIS  '  ,/, 

+'  PROGRAM.  THE  FILE  NAME  IS  LIMITED  TO  SIX  CHARACTERS,  AND',/, 

+'  SHOULD  NOT  HAVE  ANY  SPACES.  ',///, 

+'  PLEASE  ENTER  THE  DESIRED  DATA  FILE  NAME:  ' ,2X  ) 


128 


1 


x. 


READ(*,804)  NAME 
804  F0RMAT(A6) 

552  WRITE(*,910)NAME 

910  F0RMAT( /////, '  YOU  SELECTED  ',A6,'  FOR  YOUR  DATA  FILE  NAME 

+’,//) 

911  VRITE(*,811) 

811  FORMATC'  IS  THIS  THE  DESIRED  SELECTION?  ENTER  Y  FOR  YES  AND  N 

+FOR  NO:  ' ,2X,  ) 

READ(*,812)  ANS 

812  FORMATC Al) 

C 

IFCANS. EQ.  'N' )  THEN 
GOTO  898 
ELSE 

CONTINUE 

ENDIF 

IFCANS.  EQ.  'Y')  THEN 
GOTO  897 
ELSE 

CALL  CLS 
GOTO  552 
897  ENDIF 


C  ALLOW  THE  USER  TO  PROVIDE  TITLE  LINE  FOR  DATA  FILE 

833  CALL  CLS 

WRITEC*,805) 

805  FORMATC////,'  ENTER  THE  DESIRED  TITLE  TO  BE  PLACED  ON  LINE 

+’,/,'  NUMBER  ONE  OF  THE  OUTPUT  DATA  FILE: ',///, 

+'  ’.2X.  ) 


READC*,806)  DATAF 

806 

FORMATC A7 9) 

835 

WRITEC*,831) 

831 

FORMATC////,' 

DO  YOU  WISH 

TO  CHANGE  THE 

TITLE  OF 

YOUR  OUTPU 

+T  DATA  FILE?' ,/,' 

ENTER  Y 

FOR  YES  AND  N 

FOR  NO: 

',2X,  ) 

READC*,832)  ANS 

832 

c 

FORMATC Al) 

\j 

IFCANS.  EQ.  'Y')  THEN 

GOTO  833 
ELSE 


CONTINUE 

ENDIF 

IFCANS.  EQ.  'N' )  THEN 
GOTO  834 
ELSE 

GOTO  835 
834  ENDIF 
C 


C**************coeFFICIENTS  for  epoxy  and  copper  layers***************** 

Q************************  ************************************************ 


c 


DO  90  1=1, GPL 
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N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

IE=4*gPL+I 

IF=5*GPL+I 

IG=6*GPL+I 

IH=7*GPL+I 

C 

C********Vf************************cORNERS******************************** 

C 

IF  ((GIH(I).EQ.  l.OR.  GIH(I).EQ.  GDEEP).AND.  (GJH(I).EQ.  1.  OR.  GJH(I).  EQ 
+. GWIDE))  THEN 

DETERMINE  COEFFICIENTS  FOR  TOP  COPPER 
IF  (GEAT(GIH(I),GJH(I)).Eq.  0. 0)  THEN 
GC0N(I,N)  =  6 

ELSE  ^ 

GC0N(I,N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GC0N(IB,N)  =  6 
GCON(ID.N)  =  6 
GC0N(IC,N)  =  6 
GC0N(ID,N)  =  6 
GC0N(IE,N)  =  6 
GC0N(IF,N)  =  6 
GC0N(IG,N)  =  6 
GC0N(IH,N)  =  6 


LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 
LEFT  EDGE 

IF  (GJH(I). EQ. 1)  THEN 
C  LEFT  COEFFICIENT 

G0EF(I,N)  =  GYLR41 
GOEF(IC,N)  =  GYLR42 
G0EF(IE,N)  =  GYLR43 

G0EF(IG,N)  =  GYLR44  ( 

GOEF(IB,N)  =  GYLR 

G0EF(ID,N)  =  GYLR 

GOEF(IF,N)  =  GYLR 

G0EF(IG,N)  =  GYLR 

G0EF(IH,N)  =  GYLR 

N=N+1 

GCONd.N)  =  7551 
GC0N(IB,N)  =  7551 
GC0N(IC,N)  =  7551 
GCON(ID,N)  =  7551 
GC0N(IE,N)  =  7551 
GC0N(IF,N)  =  7551 
GC0N(IG,N)  =  7551 
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GCON(IH,N)  =  7551 


C  RIGHT  COEFFICIENT 
G0EF(I,N)  =  GYY41 
G0EF(IC,N)  =  GYY42 
GO£F(IE,N)  =  GYY43 
GOEF(IG,N)  =  GYY44 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
GOEF(IH,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10  *  (I+l)  +  1 
f  GC0N(IB,N)  =  10*  (IB+1)  +  1 

GC0N(IC,N)  =  10  *  (IC+1)  +  1 

GC0N(ID,N)  =  10  *  (ID+1)  +  1 

^  GC0N(IE,N)  =  10  *  (IE+1)  +  1 

GC0N(IF,N)  =  10  *  (IF+1)  +  1 

GC0N(IG,N)  =  10  *  (IG+1)  +  1 

GC0N(IH,N)  =  10  *  (IH+1)  +  1 


RIGHT  EDGE 

ELSEIF  (GJH(I).EQ. GWIDE)  THEN 

LEFT  COEFFICIENT 

GOEFd.N)  =  GYY41 
G0EF(IC,N)  =  GYY42 
G0EF(IE,N)  =  GYY43 
G0EF(IG,N)  =  GYY44 
G0EF(IB,N)  =  GYY 
G0EF(ID,N)  =  GYY 
G0EF(IF,N)  =  GYY 
G0EF(IH,N)  =  GYY 
N=N+1 

GC0N(I,N)  =  10*(I-1)+1 
GC0N(IB,N)  =  10*(IB-1)+1 
GC0N(IC,N)  =  10*(IC-1)+1 
GC0N(ID,N)  =  10*(ID-1)+1 
GC0N(IE,N)  =  10*(IE-1)+1 
GC0N(IF,N)  =  10*(IF-1)+1 
GC0N(IG,N)  =  10*(IG-1)+1 
GC0N(IH,N)  =  10*(IH-1)+1 

RIGHT  COEFFICIENT 

GOEFd.N)  =  GYLR41 
GOEFCIC.N)  =  GYLR42 
GOEFCIE.N)  =  GYLR43 
GOEFCIG.N)  =  GYLR44 
GOEFCIB.N)  =  GYLR 
GOEFCID.N)  =  GYLR 
GOEFCIF.N)  =  GYLR 
GOEFCIH.N)  =  GYLR 
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GCONd.N)  =  7541 
GCON(IB,N)  =  7541 
GCON(IC,N)  =  7541 
GCON(ID,N)  =  7541 
GC0N(IE,N)  =  7541 
GCON(IF,N)  =  7541 
GCON(IG,N)  =  7541 
GCON(IH,N)  =  7541 
ENDIF 

C 

C  FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

C 

C  FRONT  EDGE 

IF  (GIH(I).EQ. 1)  THEN 

C  FRONT  COEFFICIENT 

G0EF(I,N)=GFXB41 
G0EF(IC,N)=GFXB42 
G0EF(IE,N)=GFXB43 
G0EF(IG,N)=GFXB44 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF,N)=GFXB 
GOEF(IH,N)=GFXB 
N=N+1 

GCONd.N)  =  7521 
GCONCIB.N)  =  7521 
GCONCIC.N)  =  7521 
GCONCID.N)  =  7521 
GCONCIE.N)  =  7521 
GCONCIF.N)  =  7521 
GCONCIG.N)  =  7521 
GCONCIH.N)  =  7521 

C  BACK  COEFFICIENT 
GOEFd,N)=GXX41 
G0EFdC,N)=GXX42 
GOEFdE,N)=GXX43 
G0EFdG,N)=GXX44 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF,N)=GXX 
GOEFdH,N)=GXX 
N=N+1 

GCONd.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 
GCONdC.N)  =  10*dC+GWIDE)+l 
GCONCID.N)  =  10*CID+GWIDE)+1 
GCONCIE.N)  =  10*CIE+GWIDE)+1 
GCONCIF.N)  =  10*CIF+GWIDE)+1 


C  BACK  EDGE 

ELSEIF  CGIHCI).EQ.GDEEP)  THEN 
C  FRONT  COEFFICIENT 

G0EFCI.N)=GXX41 
GOEFCIC.N)=GXX42 
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G0EF(IE,N)=GXX43 

G0EF(IG,N)=GXX44 

GOEF(IB,N)=GXX 

GOEF(ID,N)=GXX 

GOEF(IF,N)=GXX 

GOEF(IH,N)=GXX 

N=N+1 

GCON(I,N)  =  10*(I-GWIDE)+1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 
GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCON(ID,N)  =  10*(ID-GWIDE)+1 
GCON(IE,N)  =  10*(IE-GWIDE)+1 
GCON(IF,N)  =  10*(IF-GWIDE)+1 
GCON(IG,N)  =  10*(IG-GWIDE)+1 
GCON(IH,N)  =  10*(IH-GWIDE)+1 
C  BACK  COEFFICIENT 
GOEF(I,N)=GFXB41 
G0EF(IC,N)=GFXB42 
GOEF(IE,N)=GFXB43 
G0EF(IG,N)=GFXB44 
GOEF(IB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF.N)=GFXB 
GOEF(IH,N)=GFXB 
N=N+1 

GCON(I,N)  =  7511 
GCON(IB,N)  =  7511 
GCON(IC,N)  =  7511 
GCON(ID,N)  =  7511 
GCON(IE,N)  =  7511 
GCON(IF,N)  =  7511 
GCON(IG,N)  =  7511 
GCON(IH,N)  =  7511 
ENDIF 

TOP  COEFFICIENT 
G0EF(I,N)=GZE41 
G0EF(IC,N)=GZE42 
G0EF(IE,N)=GZE43 
G0EF(IG,N)=GZE44 
G0EF(IB,N)=GZC41 
GOEF(ID,N)=GZC42 
GOEF(IF,N)=GZC43 
G0EF(IH,N)=GZC44 
N=N+1 

GCON(I,N)  =  7511 
GCON(IC,N)  =  7511 
GCON(IE,N)  =  7511 
GCON(IB,N)  =  10*1+1 
GCON(ID,N)  =  10*IC+1 
GCON(IF,N)  =  10*IE+1 
GCON(IG,N)  =  7511 
GCON(IH,N)  =  10*IF+1 

BOTTOM  COEFFICIENT 
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GOEFd.N)  =  GZE41 
GOEFdC.N)  =  GZE42 
GOEFCIE.N)  =  GZE43 
GOEFCIG.N)  =  GZE44 
GOEFCIB.N)  =  GZC41 
GOEFCID.N)  =  GZC42 
GOEFCIF.N)  =  GZC43 
GOEFdH.N)  =  GZB 
N=N+1 

GCONCI.N)  =10*d+GPL)+l 
GCONCIB.N)  =  10*dB+GPL)+l 
GCONCIC.N)  =  10*(IC+GPL)+1 
GCONCID.N)  =  10*dD+GPL)+l 
GCONCIE.N)  =  10*(IE+GPL)+1 
GCON(IF,N)  =  10*dF+GPL)+l 
GCONCIG.N)  =  10*dG+GPL)+l 
GCONdH.N)  =  10*dH+GPL)+l 
C  HEAT  INPUT 

IF(GEAT(GIHd),GJHd)).NE.  0.  0)  THEN 
GOEFd.N)  =  GEAT(GIHd),GJHd)) 
N=N+1 

GC0Nd,N)=9991 

ENDIF 


C 

ELSEIF((GIHd).EQ.  1.  OR.  GIH(  I).  EQ.  GDEEP).AND.  (GJHd).NE.  l.OR.  GJHfl) 
+.NE.GWIDE))  THEN 

DETERMINE  NUMBER  OF  CONNECTIONS  FOR  COPPER  AND  EPOXY  LAYERS 
IF  (GEAT(GIHd),GJHd)).EQ.  0.0)  THEN 
GCONd.N)  =  6 

ELSE 

GCONCI.N)  =  7 
ENDIF 

CONNECTIONS  FOR  EPOXY  LAYER 
GCONCIB.N)  =  6 
GCONCID.N)  =  6 
GCONCIC.N)  =  6 
GCONCIE.N)  =  6 
GCONCIF.N)  =  6 
GCONCIC.N)  =  6 
GCONCIH.N)  =  6 
C 

C  LEFT  AND  RIGHT  COEFFICIENTS 

C  LEFT  COEFFICIENT 

GOEFCI.N)  =  GYY41 
GOEFCIC.N)  =  GYY42 
GOEFCIE.N)  =  GYY43 
GOEFCIG.N)  =  GYY44 
GOEFCIB.N)  =  GYY 
GOEFCID.N)  =  GYY 
GOEFCIF.N)  =  GYY 
GOEFCIH.N)  =  GYY 
N=N+l 

GCONCI.N)  =  10*CI-1)+1 
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GC0N(IB,N)  =  10*(IB-1)+1 
GCON(IC,N)  =  10*(IC-1)+1 
GC0N(ID,N)  =  10*(ID-1)+1 
GC0N(IE,N)  =  10*(IE-1)+1 
GC0N(IF,N)  =  10*(IF-1)+1 
GC0N(IG,N)  =  10*(IG-1)+1 
GCON(IH,N)  =  10*(IH-1)+1 

RIGHT  COEFFICIENT 
GOEF(I,N)  =  GYY41 
GOEF(IC,N)  =  GYY42 
G0EF(IE,N)  =  GYY43 
GOEF(IG,N)  =  GYY44 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
GOEF(IH,N)  =  GYY 
N=N+1 

GCONd.N)  =  10  *  (I+l)  +  1 
GCON(IB,N)  =  10*  (IB+1)  +  1 
GCON(IC,N)  =  10  *  (IC+1)  +  1 
GC0N(ID,N)  =  10  *  (ID+1)  +  1 
GCON(IE,N)  =  10  *  (IE+1)  +  1 
GC0N(IF,N)  =  10  *  (IF+1)  +  1 
GCONCIG.N)  =  10  *  (IG+1)  +  1 
GCONdH.N)  =  10  *  (IH+l)  +  1 

FRONT  AND  BACK  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

FRONT  EDGE 

IF  (GIHd).EQ.  1)  THEN 
C  FRONT  COEFFICIENT 
G0EFd,N)=GFXB41 
G0EFdC,N)=GFXB42 
GOEFdE,N)=GFXB43 
G0EFdG,N)=GFXB44 
GOEFdB,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEFdF,N)=GFXB 
GOEFdH,N)=GFXB 
N=N+1 

GCONd.N)  =7521 
GCONCIB.N)  =  7521 
GCONdC.N)  =  7521 
GCONCID.N)  =  7521 
GCONCIE.N)  =  7521 
GCONCIF.N)  =  7521 
GCONCIG.N)  =  7521 
GCONCIH.N)  =  7521 

C  BACK  COEFFICIENT 
G0EFd.N)=GXX41 
G0EFdC.N)=GXX42 
G0EFdE.N)=GXX43 
G0EFdG.N)=GXX44 
GOEFdB.N)=GXX 
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GOEF(ID,N)=GXX 

GOEF(IF,N)=GXX 

GOEF(IH,N)=GXX 

N=N+1 

GCONd.N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GCON(IC,N)  =  10*(IC+GWIDE)+1 
GCON(ID,N)  =  10*(ID+GVIDE)+1 
GCON(IE,N)  =  10*(IE+GWIDE)+1 
GCON(IF,N)  =  10*(IF+GWIDE)+1 
GCONCIG.N)  =  10*(IG+GWIDE)+1 
GCON(IH,N)  =  10*(IH+GWIDE)+1 
C  BACK  EDGE 

ELSEIF  (GIHCI).EQ. GDEEP)  THEN 
C  FRONT  COEFFICIENT 

G0EF(I,N)=GXX41 
G0EF(IC,N)=GXX42 
GOEF(IE,N)=GXX43 
G0EF(IG,N)=GXX44 
GOEF(IB,N)=GXX 
GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
GOEF(IH,N)=GXX 
N=N+1 

GCONd.N)  =10*d-GWIDE)+l 
GCONCIB.N)  =  10*dB-GWIDE)+l 
GCONCIC.N)  =  10*dC-GWIDE)+l 
GCONCID.N)  =  10*dD-GVIDE)+l 
GCONCIE.N)  =  10*dE-GWIDE)+l 
GCONCIF.N)  =  10*dF-GWIDE)+l 
GCONCIG.N)  =  10*dE-GWIDE)+l 
GCONCIH.N)  =  10*dF-GWIDE)+l 
C  BACK  COEFFICIENT 
G0EFd,N)=GFXB41 
GOEFdC,N)=GFXB42 
G0EFdE,N)=GFXB43 
G0EFdG,N)=GFXB44 
GOEFdB,N)=GFXB 
GOEFdD,N)=GFXB 
GOEFdF,N)=GFXB 
GOEFdH,N)=GFXB 
N=N+1 

GCONCI.N)  =  7511 
GCONdB.N)  =  7511 
GCONCIC.N)  =  7511 
GCONCID.N)  =  7511 
GCONCIE.N)  =  7511 
GCONCIF.N)  =  7511 
GCONCIG.N)  =7511 
GCONCIH.N)  =  7511 
ENDIF 

TOP  COEFFICIENT 
G0EFCI,N)=GZE41 
GOEFCIC.N)=GZE42 
G0EFCIE.N)=GZE43 
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G0EF(IG,N)=GZE44 

G0EF(IB,N)=GZC41 

G0EF(ID,N)=GZC42 

G0EF(IF,N)=GZC43 

G0EF(IH,N)=GZC44 

N=N+1 

GCON(I,N)  =  7511 
GC0N(IC,N)  =  7511 
GCON(IB,N)  =  10*1+1 
GC0N(ID,N)  =  10*IC+1 
GC0N(IE,N)  =  7511 
GC0N(IF,N)  =  10*IE+1 
GC0N(IG,N)  =  7511 
GC0N(IH,N)  =  10*IG+1 

BOTTOM  COEFFICIENT 

G0EF(I,N)  =  GZE41 
G0EF(IC,N)  =  GZE42 
G0EF(IE,N)  =  GZE43 
G0EF(IG,N)  =  GZE44 
G0EF(IB,N)  =  GZC41 
G0EF(ID,N)  =  GZC42 
G0EF(IF,N)  =  GZC43 
G0EF(IH,N)  =  GZB 
N=N+1 

GC0N(I,N)  =10*(I+GPL)+1 
GC0N(IB,N)  =  10*(IB+GPL)+1 
GC0N(IC,N)  =  10*(IC+GPL)+1 
GC0N(ID,N)  =  10*(ID+GPL)+1 
GC0N(IE,N)  =  10*(IE+GPL)+1 
GC0N(IF,N)  =  10*(IF+GPL)+1 
GC0N(IG,N)  =  10*(IG+GPL)+1 
GC0N(IH,N)  =  10*(IH+GPL)+1 

HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0.  0)  THEN 
GOEFd.N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 

GC0N(I,N)=9991 
END  IF 


LEFT  AND  RIGHT  EDGES  EXCLUDING  CORNERS 

ELSEIF((GIH(I).EQ.  1.  OR.  GIH( I ).  EQ.  GDEEP).  AND.  (GJH(I).EQ.  l.OR.GJH(I) 
+. EQ. GWIDE))  THEN 

DETERMINE  COEFFICIENTS  FOR  TOP  LAYER 
IF  (GEAT(GIHa'>,GJH(I)).NE.  0.  0)  THEN 
GCON(I,N)  ==  7 

ELSE 

GCON(I,N)  =  6 
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ENDIF 


CONNECTIONS  FOR  EPOXY  LAYER 
GCON(IB,N)  =  6 
GCON(IC,N)  =  6 
GCON(ID,N)  =  6 
GCON(IE,N)  =  6 
GCON(IF,N)  =  6 
GCON(IG,N)  =  6 
GC0N(IH,N)  =  6 

LEFT  AND  RIGHT  COEFFICIENTS  DEPENDING  ON  WHICH  EDGE 

LEFT  EDGE 

IF  (GJH(I).EQ.  1)  THEN 

LEFT  COEFFICIENT 
GOEFfl.N)  =  GYLR41 
GOEF(IC,N)  =  GYLRA2 
GOEF(IE,N)  =  GYLR43 
GOEF(IG,N)  =  GYLR44 
GOEF(IB,N)  =  GYLR 
GOEF(ID,N)  =  GYLR 
G0EF(IF,N)  =  GYLR 
GOEF(IH,N)  =  GYLR 
N=N+1 

GCON(I,N)  =  7551 
GCON(IB,N)  =  7551 
GCON(IC,N)  =  7551 
GC0N(ID,N)  =  7551 
GCON(IE,N)  =  7551 
GCON(IF,N)  =  7551 
GCON(IG,N)  =  7551 
GCON(IH,N)  =  7551 

RIGHT  COEFFICIENT 
G0EF(I,N)  =  GYY41 
G0EF(IC,N)  =  GYY42 
GOEF(IE,N)  =  GYY43 
GOEF(IG,N)  =  GYY44 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
GOEF(IH,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10  *  (I+l)  +  1 
GC0N(IB,N)  =  10*  (IB+1)  +  1 
GCONdC.N)  =  10  *  (IC+1)  +  1 

GCON(ID,N)  =  10  *  (ID+1)  +  1 

GCON(IE,N)  =  10  *  (IE+1)  +  1 

GC0N(IF,N)  =  10  *  (IF+1)  +  1 

GC0N(IG,N)  =  10  *  (IG+1)  +  1 

GC0N(IH,N)  =  10  *  (IH+1)  +  1 

RIGHT  EDGE 

ELSEIF  (GJH(I).EQ.  GWIDE)  THEN 

LEFT  COEFFICIENT 
GOEF(I,N)  =  GYY41 
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GOEF(IC,N)  =  GYY42 
GOEF(IE,N)  =  GYY43 
G0EF(IG,N)  =  GYY44 
GOEF(IB,N)  =  GYY 
GOEF(ID,N)  =  GYY 
GOEF(IF,N)  =  GYY 
GOEF(IH,N)  =  GYY 
N=N+1 

GCON(I,N)  =  10*(I-1)+1 
GCON(IB,N)  =  10*(IB-1)+1 
GCON(IC,N)  =  10*(IC-1)+1 
GCON(ID,N)  =  10*(ID-1)+1 
GCON(IE,N)  =  10*(IE-1)+1 
GCON(IF.N)  =  10*(IF-1)+1 
GCON(IG,N)  =  10*(IG-1)+1 
GCON(IH,N)  =  10*(IH-1)+1 
RIGHT  COEFFICIENT 

GOEF(I,N)  =  GYLR41 
GOEF(IC,N)  =  GYLR42 
GOEFdE.N)  =  GYLR43 
GOEF(IG,N)  =  GYLR44 
G0EF(IB,N)  =  GYLR 
G0EF(ID,N)  =  GYLR 
GOEF(IF,N)  =  GYLR 
GOEF(IH,N)  =  GYLR 
N=N+1 

GCONd.N)  =  7541 
GCONCIB.N)  =  7541 
GCONdC,N)  =  7541 
GCONCID.N)  =  7541 
GCONCIE.N)  =  7541 
GCONCIF.N)  =  7541 
GCONCIG.N)  =  7541 
GCONCIH.N)  =  7541 
END  IF 

FRONT  COEFFICIENT 
GOEFCI.N)  =GXX41 
G0EFdC,N)=GXX42 
GOEFdE,N)=GXX43 
G0EFdG,N)=GXX44 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF,N)=GXX 
GOEFdH,N)=GXX 
N=N+1 

GCONCI.N)  =10*d-GWIDE)+l 
GCONCIB.N)  =  10*dB-GWIDE)+l 
GCONCIC.N)  =  10*dC-GWIDE)+l 
GCONCID.N)  =  10*dD-GWIDE)+l 
GCONCIE.N)  =  10*CIE-GWIDE)+1 
GCONCIF.N)  =  10*CIF-GWIDE)+1 
GCONCIG.N)  =  10*CIG-GWIDE)+1 
GCONCIH.N)  =  10*CIH-GWIDE)+1 
C  BACK  COEFFICIENT 
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non 


G0EF(I,N)=GXX41 

G0EF(IC,N)=GXX42 

G0EF(IE,N)=GXX43 

G0EF(IG,N)=GXX44 

GOEF(IB,N)=GXX 

GOEF(ID,N)=GXX 

GOEF(IF,N)=GXX 

GOEF(IH,N)=GXX 

N=N+1 

GCON(I,N)  =10*(I+GWIDE)+1 
GCON(IB,N)  =  10*(IB+GWIDE)+1 
GCON(IC,N)  =  10*(IC+GWIDE)+1 
GCON(ID,N)  =  10*(ID+GWIDE)+1 
GCON(IE,N)  =  10*(IE+GWIDE)+1 
GCON(IF,N)  =  10*(IF+GWIDE)+1 
GC0N(IG,N)  =  10*(IG+GWIDE)+1 
GCON(IH,N)  =  10*(IH+GWIDE)+1 
C  TOP  COEFFICIENT 
G0EF(I,N)=GZE41 
G0EF(IC,N)=GZE42 
G0EF(IE,N)=GZE43 
G0EF(IG,N)=GZE44 
G0EF(IB,N)=GZC41 
G0EF(ID,N)=GZC42 
G0EF(IF,N)=GZC43 
G0EF(IH,N)=GZC44 
N=N+1 

GCON(I,N)  =7511 
GCON(IC,N)  =  7511 
GCON(IB,N)  =  10*1+1 
GC0N(ID,N)  =  10*IC+1 
GC0N(IE,N)  =  7511 
GC0N(IF,N)  =  10*IE+1 
GCCN(IG,N)  =  7511 
GC0N(IH,N)  =  10*IG+1 

BOTTOM  COEFFICIENT 


G0EF(I,N) 

G0EF(IC,N) 

G0EF(IE,N) 

G0EF(IG,N) 

G0EF(IB,N) 

G0EF(ID,N) 

G0EF(IF,N) 

G0EF(IH,N) 

N=N+1 

GCONd.N) 

GCON(IB,N) 

GCON(IC,N) 

GCON(ID,N) 

GC0N(IE,N) 

GCON(IF,N) 

GCON(IG,N) 


GZE41 

GZE42 

GZE43 

GZE44 

GZC43 

GZC43 

GZC43 

EZB 

10*(I+GPL)+1 

10*(IB+GPL)+1 

10*(IC+GPL)+1 

10*(ID+GPL)+1 

10*(IE+GPL)+1 

10*(IF+GPL)+1 

10*(IG+GPL)+1 


140 


o  Cl  o  on  oo  oon  no 


GCON(IH,N)  =  10*(IH+GPL)+1 
HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0.0)  THEN 
GOEFd.N)  =  GEAT(GIH(I),GJH(I)) 
N=N+1 

GC0N(I,N)=9991 

ENDIF 


****dETERMINE  COEFFICIENTS  FOR  ALL  NODES  NOT  TOUCHING  AN  EDGE******* 


ELSEIF((GIH(I).NE.  1.  OR.  GIH( I).  NE.  GDEEP).  AND.  (GJH(I).NE.  l.OR.  GJH(I) 
+. NE. GWIDE))  THEN 

DETERMINE  CONNECTIONS  FOR  TOP  LAYER 
IF  (GEAT(GIH(I),GJH(I)).NE.O.  0)  THEN 
GCONd.N)  =  7 

ELSE 

GC0Nd,N)  =  6 
ENDIF 


CONNECTIONS  FOR  EPOXY  LAYER 
GCONCIB.N)  =  6 
GCONCIC.N)  =  6 
GCONCID.N)  =  6 
GCONCIE.N)  =  6 
GCONCIF.N)  =  6 
GCONCIG.N)  =  6 
GCONCIH.N)  =  6 


C 


LEFT  COEFFICIENT 
GOEFCI.N)  =  GYY41 
GOEFCIC.N)  =  GYY42 
GOEFCIE.N)  =  GYY43 
GOEFCIG.N)  =  GYY44 
GOEFCIB.N)  =  GYY 
GOEFCID.N)  =  GYY 
GOEFCIF.N)  =  GYY 
GOEFCIH.N)  =  GYY 
N=N+1 


GCON( 

GCON( 

GCON( 

GC0N( 

GCON( 

GCON( 

GCON( 

GCON( 

RIGHT 

GOEF( 

G0EF( 

GOEF( 

G0EF( 

G0EF( 

G0EF( 


I,N) 

IB, N) 

IC, N) 

ID, N) 
lE.N) 

IF. N) 

IG. N) 

IH. N) 


10  *  (I-II  +  1 
10*  (IB-l)  +  1 
10  *  (IC-l)  +  1 
(ID-l)  +  1 
(lE-l)  +  1 
(IF-l)  +  1 
(IG-l)  +  1 
(IH-l)  +  1 


10 

10 

10 

10 

10 


if 

if 

if 

if 

if 


COEFFICIENT 
I.N)  =  GYY41 

IC, N)  =  GYY42 
IE,N)  =  GYY43 
IG,N)  =  GYY44 
IB,N)  =  GYY 

ID, N)  =  GYY 
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o  o 


GOEF(IF,N)  =  GYY 
GOEF(IH,N)  =  GYY 
N=N+1 

GCONd.N)  =  10  *  (I+l)  +  1 
GCON(IB,N)  =  10*  (IB+1)  +  1 
GCON(IC,N)  =  10  *  (IC+1)  +  1 
GCON(ID,N)  =  10  *  (ID+1)  +  1 
GCON(IE,N)  =  10  *  (IE+1)  +  1 
GCON(IF,N)  =  10  *  (IF+1)  +  1 
GCON(IG,N)  =  10  *  (IG+1)  +  1 
GCON(IH,N)  =  10  *  (IH+1)  +  1 

FRONT  COEFFICIENT 
G0EF(I,N)=GXX41 

G0EF(IC,N)=GXX42  , 

G0EF(IE,N)=GXX43  ^ 

G0EFCIG,N)=GXX44 

GOEF(IB,N)=GXX  < 

GOEF(ID,N)=GXX 
GOEF(IF,N)=GXX 
GOEF(IH,N)=GXX 
N=N+1 

GCON(I,N)  =10*(I-GWIDE)+1 
GCON(IB,N)  =  10*(IB-GWIDE)+1 
GCON(IC,N)  =  10*(IC-GWIDE)+1 
GCON(ID,N)  =  10*(ID-GWIDE)+1 
GCONdE.N)  =  10*dE-GWIDE)+l 
GCONCIF.N)  =  10*dF-GWIDE)  +  l 
GCONCIG.N)  =  10*dG-GWIDE)+l 
GCONCIH.N)  =  10*dH-GWIDE)+l 
C  BACK  COEFFICIENT 
G0EFd,N)=GXX41 
G0EFdC,N)=GXX42 
G0EFdE,N)=GXX43 
G0EFdH,N)=GXX44 
GOEFdB,N)=GXX 
GOEFdD,N)=GXX 
GOEFdF,N)=GXX 
GOEFdH,N)=GXX 
N=N+1 

GCONCI.N)  =10*d+GWIDE)+l 
GCONCIB.N)  =  10*dB+GWIDE)+l 

GCONCIC.N)  =  10*dC+GWIDE)+l  , 

GCONCID.N)  =  10*dD+GWIDE)+l 

GCONCIE.N)  =  10*dE+GWIDE)+l 

GCONdF,N)  =  10*dF+GWIDE)+l  v 

GCONCIG.N)  =  10*dG+GWIDE)+l 

GCONCIH.N)  =  10*dH+GWIDE)+l 

C  TOP  COEFFICIENT 
G0EFd,N)=GZE41 
G0EFdC,N)=GZE42 
G0EFdE,N)=GZE43 
G0EF(IG,N)=GZE44 
G0EFdB,N)=GZC41 
G0EFdD,N)=GZC42 
G0EFdF,N)=GZC43 
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G0EF(IH,N)=GZC44 

N=N+1 

GCON(I,N)  =7511 
GCON(IB,N)  =  10*1+1 
GC0N(IC,N)  =  7511 
GC0N(ID,N)  =  10*IC+1 
GCON(IE,N)  =  7511 
GCON(IF,N)  =  10*IE+1 
GC0N(IG,N)  =  7511 
GCON(IH.N)  =  10*IG+1 
C 

C  BOTTOM  COEFFICIENT 

C 

GOEF(I,N)  =  GZE41 
GOEF(IC,N)  =  GZE42 
GOEF(IE,N)  =  GZE43 
GOEF(IG,N)  =  GZE44 
GOEF(IB,N)  =  GZC41 
GOEF(ID,N)  =  GZC42 
GOEF(IF,N)  =  GZC43 
GOEF(IH,N)  =  EZB 
N=N+1 

GCON(I,N)  =10*(I+GPL)+1 
GCON(IB,N)  =  10*(IB+GPL)+1 
GCON(IC,N)  =  10*(IC+GPL)+1 
GCON(ID,N)  =  10*(ID+GPL)+1 
GCON(IE,N)  =  10*(IE+GPL)+1 
GCONdF.N)  =  10*(IF+GPL)+1 
GCON(IG,N)  =  10*(IG+GPL)+1 
GCON(IH,N)  =  10*(IH+GPL)+1 
C  HEAT  INPUT 

IF(GEAT(GIH(I),GJH(I)).NE.  0.0)  THEN 
GOEF(I,N)  =  GEAT(GIH(I) ,GJH(I)) 

N=N+1 

GC0N(I,N)=9991 

END  IF 

ENDIF 

90  CONTINUE 

C 

C  GENERATE  DATA  FILE 

C 

OPEN  ( 3, F I LE=NAME,FORM=' FORMATTED' ,ACCESS=' DIRECT' ,RECL=108 , STATUS 
+='NEW' ) 

WRITE(3,909)  DATAF 
909  F0RMAT(1X,A79) 

WRITE( 3 , 908 )  COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908  FORMAT(2X,9(I3,5X)) 

WRITE(3,907)  ZER, ZER. ZER 
907  F0RMAT(2X,3(I3,5X)) 

WRITE(3,9081)  NMAX,TMAX,HTRS,D1,D2,D3,D4,D5,D6,D7 
9081  FORMAT(2X,9(I3,5X)) 

WRITE(3,905)ACC,DAMP,MAXIT,CONFAC,GIBT 

905  F0RMAT(1X,2(F9. 7 , IX) , 14X , 12 , 1X,F9.  7 , 1X,F9.  5) 

WRITE ( 3 , 906 )  GPRT , GFT , GBT , GRT , GLT , GWRT 

906  F0RMAT(1X,6(F12. 3,1X)) 

DO  112  I=1,L0CVAR 


WRITE( 3,9100)  GC0N(I,1) ,GCON( I ,2) ,GCON( I , 3) ,GCON( I ,4) ,GCON( I ,5) ,GC 
+ON(I,6),GCON(I,7) ,GCON(I,8) 

9100  FORMAT(I4,3X,7(I5,7X)) 

WRITE(3,9110)  G0EF(I,1),G0EF(I,2),G0EF(I,3),G0EF(I,4),G0EF(I,5),G0 
+EF(I,6),GOEF(I,7) 

911^  FORMAT(7(F9.  3,3X)) 

112  CONTINUE 
CLOSE  (3) 

CALL  CLS 

WRITE(*,999)  NAME 

999  FORMATC////, '  THE  OUTPUT  DATA  HAS  BEEN  PLACED  IN  A  FILE 

+NAMED  ',A6, /////,'  <PRESS  ENTER  TO  CONTINUE>') 

READ(*,5912)  ANS 
5912  FORMAT(Al) 

END 

C  ********************************************************************** 
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